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SUMMARY 


Bends of large width occasionally exhibit cracks at the center 
of the width when minimum bend radii successfully employed on 
small widths are used. An analysis of the stress condition exist- 
ing at the center of a large width revealed an unfavorable bi- 
axial tension stress that reduces the ultimate elongation from 
that in simple tension. An experimental phctogrid strain survey 
on three aluminum alloys confirmed the results of the analysis and 
indicated the extent of bi-axiality of stress as the width 
decreased. 

Elongation date from standard tensile coupons were corrected 
for bi-axiality, and the results were then applied in determining 
theoretical minimum bend radii for various widths of bend. 


INTRODUCTION 


rag A RECENT INVESTIGATION at Republic on 
the formability properties of certain new alumi- 
num alloys, attention was given to the basic forming 


problem of minimum bend radii. In evaluating the 
bend radius test results of this investigation, an attempt 
was made to correlate the data on Alclad 24S-O (in- 
cluded in the program to provide a basis for comparison 
with the new materials) with the minimum bend radii 
in current use at various aircraft companies. 

The investigation indicated a marked disagreement 
on the minimum bend radius for various sheet gages of 
Alclad 24S-O, although apparently the stated value is 
the minimum that can be consistently attained in the 
shops of the various companies. A large share of this 
disagreement could probably be traced to the different 
handling techniques employed in each shop, although 
substantially the same equipment is employed every- 
where. Another source of discrepancy may be the con- 
sistency of results desired—i.e., some companies were 
willing to tolerate a small number of rejects in order to 
realize a smaller bend radius, while others were not. 
The number of rejects is probably caused by the vari- 
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ations in material properties from one lot of sheet to 
another. 

In spite of these possible causes for disagreement, it 
was felt that some fundamental factor in establishing 
minimum bend radii was responsible for the dis- 
crepancy. A probable answer was indicated in refer- 
ence 1, in which test data were presented which showed 
a marked increase in minimum bend radius as the speci- 
men width (the length of bend normal to the bend cross 
section) was increased. Since each company probably 
used different bend widths in establishing minimum 
bend radii for Alclad 24S-O, this effect as given in 
reference 1 was serious enough to account for the 
discrepancy. 

Aside from providing a reason for the disagreement 
on the Alclad 24S-O bend radii, these data provided an 
explanation of why sections of large bend width, such 
as stringers or longerons, occasionally crack in central 
regions when bend radii successfully employed on short 


Characteristic cracking at center of bend width caused 


Fic. 1. 
by bi-axial tension stress condition. 
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widths are used. This phenomenon is illustrated 
in Fig. 1 in which a specimen of 0.064-in. R301-T 
aluminum alloy was formed to a 90° bend. This in- 
vestigation was initiated, therefore, to determine the 
seriousness of the width effect on minimum bend radii. 


SYMBOLS 


= stress, lbs. per sq. in. 

= strain, in. per in. 

= coordinate in direction of bend 

= coordinate in width direction 

= coordinate in thickness direction 

= Poisson’s ratio, 0.50 in the plastic range 

= bend width in 

= bend radius to neutral axis, in. 

= inside bend radius, in. 

= thickness, in. 

= ratio of distance between qutside fiber of bend and 
neutral axis to total thickness 

= ratio of cladding thickness per face to total thickness 

= strain in bend direction 

= strain in width direction 

= strain in bend direction under simple tension, i.e., plane 

stress 
éuit. corrected for cladding thickness 
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Stress CONDITION IN WIDE BENDS 


Before an analysis of the stress condition in a bend is 
presented, it would be well to review some basic con- 
cepts of the response of metals to combined stress. 

In general, we are only concerned with an element 
in a two-dimensional stress field, since the stress condi- 
tion in a bend can be reduced to one of this nature. 
In Fig. 2, the three principal stress conditions that can 
be obtained on a two-dimensional element are illus- 
trated and a stress-strain curve for each condition is 
shown. Each of the combined stress conditions in- 
fluences the stress-strain curve in a manner such as to 
ease or aggravate the forming operation. Condition 
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Fic. 2. Various stress states on a two-dimensional element and 
the effect upon the stress-strain curve. 


(b) approaches a pure shear condition that is helpfy 
in forming because of the high elongations that can by 
obtained, whereas condition (c) approaches the highly 
undesirable bi-axial tension condition in which elonga. 
tion properties are considerably reduced from those ip 
simple tension (a). 

The principal forming load is often applied in one 
direction only but the geometry of the part is such that 
stresses are induced in directions other than that of the 
applied load. An element elongated in one direction 
tends to contract in the other two, and if these induced 
contractions are restrained, stresses are induced in the 
latter directions. The ratio of this induced latera 
deformation to applied axial deformation under an axial 
load is Poisson’s ratio. For aluminum alloys, this 
ratio has a value of 0.3 in the elastic range and is as. 
sumed equal to 0.5 in the plastic range of the material, 

Fig. 3 shows the deformation of an element under an 
applied stress. If the lateral deformation is com. 
pletely restricted—i.e., ¢, = O—a stress will be induced 
in this direction. If no lateral restraint is imposed 
then «, = —ve, and the stress in the lateral direction 
will be zero. 
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Fic. 3. Poisson’s ratio, 


Anticlastic Bending 


In the bending process, the outer surface of the bend 
is elongated and the inner surface compressed. Thus, 
tension and compression forming loads are applied 
to the outside and inside surfaces, respectively. These 
forming loads induce strains by Poisson’s effect in the 
width (y) and thickness (z) directions (see Fig. 4) of a 
magnitude equal to one-half the applied strain, since 
forming is accomplished in the plastic range. The 
z-direction strain is not restrained by the geometry of 
the formed part, and, consequently, in the bending 
process some thinning can occur. 

The strain in the width direction, however, is af- 
fected by the bend géometry. The outer fiber strain 
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Fic. 4. Bend parameters. 

induces a lateral contraction in the width direction, and, 
similarly, the inner fiber strain induces a lateral elonga- 
tion. Unless the bend curves in the width direction 
as Shown in Fig. 5, these induced strains will be re- 
strained. The geometry of the bent part is such as to 
completely restrain this anticlastic curvature beeamse- 
of the relatively large moment of inertia of the flange 
and the high degree of curvature of the material in the 
bend region. 

Photographs of sections through bend widths are 
illustrated in Fig. 6. Except for a highly localized 
curvature at the edge of the large bend widths and the 
complete anticlastic curvature in the smallest bend 
widths, the center of the width is straight and, conse- 
quently, the width direction strain must be zero. 





Fic.5. Anticlastic curvature of bend if no restraint were imposed 
in the width direction. 
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Therefore, it can be concluded that the induced stress 
along the bend width arises from the restriction upon 
latera] strain. This induced stress corresponds to a 
strain equal to one-half the strain in the direction of 
loading. 

At the center of the tension side of the bend, there 
now exists a bi-axial tension condition that is highly 
unfavorable and reduces the formability of the ma- 
terial. Thus, for bend widths in the usual range em- 
ployed in aircraft parts, an unfavorable bi-axial ten- 
sion condition does exist in the central region of the 
width. The bi-axial stress state reduces to a simple 
tension at the edges where the anticlastic curvature is 
not restrained, and, consequently, cracking may be con- 
fined only to the central region of the bend width where 
lowest elongation exists. 








Fic. 6. Section through bend showing (a) anticlastic curvature 
confined to edge of large width and (b) complete anticlastic 
curvature of small widths (0.4 to 0.1 in.). 


BEND Rapius TESTS ON VARIOUS WIDTHS 


Test Specimens 


The bending of large widths—e.g., stringers or 
longerons—is usually accomplished in a power brake, 
whereas the hydropress and hand brake are used largely 
in bending sections of small width. To obtain a com- 
plete range of bend widths in this investigation, the 
froming of all bends was performed on the power brake. 
Thus, if the effect of bend width upon minumum bend 
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TABLE 1 
Width of Bend Specimens 





Material Bend Width, In. 
0.064 Alclad 24S-T Specimens A1/,* 1/,* 1* -2% 5 10 25 50 
0.064 R301-T Specimens B 100 
0.064 Alclad 75S-T Specimens C 0.1° 0:2" 0.3 0.4 0.5 





* Two widths of this diménsion were included on each bend specimen. 


radius were serious, it would be manifest on the power 
brake where the widest range of widths are attainable. 
To limit the scope of the investigation, the test 
program was restricted to the -T temper of aluminum 
alloys R301, Alclad 24S, and Alclad 75S. In addition, 
the specimens were of 0.064-in. sheet only, and no at- 
tempt was made to investigate the effect of thickness. 
The test specimens used in the investigation were of 
two general types: bend strips of various widths, 
which were formed on the power brake, and standard 
tensile coupons of each material cut from the same sheet 
as the width specimens. The bend specimens—strips 
of the same dimensions except for the width, which was 
varied as listed in Table 1—were designed to form 
simultaneously as many different widths as possible and 
thus eliminate experimental errors due to: variations in 
bend angle which were anticipated #f the. variotis bend Fic. 7. Photogrid at center of large width. 
widths were formed individually, <*~. » =m 
The general experiniental approach was “to obtain 
data on shop bend tests of specimens A and B to indi- 
cate the severity of tie.bend width efféct: \ Int addi- : nsorr / S/f y ?/ 1/4) 
tion, elongation measurements” in-theawidth direction 7 ‘/ L, / 7 J / / i / / “y 
on C-type specimens were made for an analysis of the 
width effect. The tensile coupons were used to obtain 
the elongation properties of the sheet material used. 
A convenient method for the determination: of plastic 
strains and strain distributions exists in the photogrid 






/ 
ALCLAD ve ; S L/ /, 


process. This process, amply described in references 2 J ///; /f// 
2 and 3, consists of photographing a fine grid on the 2 A) /) 1 fs 
sensitized surface of the specimen as illustrated in | 
Figs. 7, 10, and 12. After the specimen has’ been 3 
tested, elongations between the grid lines are measured ~ | 
with a cathetometer. An application of this tech- OY 
nique in determining minimum bend radii is described ‘ T 7S// // 
in reference 4. UA AAA ] 
As a preliminary step, a 2-in. grid was printed at the ] Bye ano 
center of each of the widths on one A-type specimen of | Rebepconi a 
Alclad 24S-T and the specimen was bent to a 1/¢-in. oo os ' 2 5 ee o 
radius. The bent specimen with the grid at the center _ — WIOTH we SUS. 
of a width is shown in Fig. 7. After bending, a strain !* ® eke with width fos 


survey on each width indicated that elongations along 
the width of the bend were restricted to specimens less was blanketed by the various widths tested. The 
than 1 in. Consequently, all subsequent elongation formed specimens were visually inspected for cracking, 
surveys were confined to C-type specimens. which was readily discernible in all cases. To discover 
The bend specimens.of type A and B were bent 90° any cracks that may have developed under the clad- 
to a number of different radii in the region of ‘the ding, bend width specimens that formed successfully 
minimum radius so as to produce fracture in some immediately after a cracked specimen of larger width 
widths and satisfactory bends near the minimum in were sectioned at the center of the width and deep 
others. Thus the region of the minimum bend radius_ etched. 
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TABLE 2 
90° Power-Brake Bend Tests on Specimens A and B—0.064-In.’Sheet 
Inside Bend Width, In. -_—_ 
Material Radius, In. R;/t 1/, 1/4 1 2 5 10 25 50 100 
Alclad 24S-T 1/16 0.97 fr fr fr fr fr fr fr fr fr 
3/39 1.47 ok ok ok ok* cr cr cr cr fr 
1/s 1.95 ok ok ok ok ok ok ok ok ok 
Alclad 75S-T 5/30 2.44 ok ok ok* fr fr fr fr fr 
3/16 2.94 ok ok ok ok ok ok* cr cr a 
7/30 3.42 ok ok ok ok ok ok ok ok ok 
1/, 3.90 a - a es de i ry: as ok 
R301-T 1/s 1.95 fr fr fr fr fr fr fr fr iN 
5/30 2.44 ok ok ok ok ok* er fr cr cr 
3/16 2.94 ok ok ok ok ok ok ok ok ok 





* Specimens sectioned at center of width and deep-etched. 
Note: ok-—specimen formed successfully; cr—specimen cracked in center region of bend; fr—specimen fractured along entire 
bend line. 


Power-brake equipment for making bends in the It was interesting to observe the progressive severity 
shop varies in steps of 1/z-in. radius. Consequently, of cracking along the bend line as the specimen width 
it was impossible to determine exactly the minimum increased. A typical example is illustrated in Fig. 9. 
bend radius for each width specimen. At best, the 


minimum radius could be bracketed, and thus the ELONGATION SURVEYS 
radius was approached within a maximum tolerance of : 
i: Type C Specimens 
#!/s9 in. 
The shop bend tests are tabulated in Table 2, and As previously stated, the elongation survey was con- 
these data are plotted on Fig. 8. fined to type C specimens. Photogrids were printed 
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Fic. 9. Progressive severity of anticlastic cracking as bend width increases. 
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widths. 





Fic. 12. Gridded tensile coupon showing enlarged fracture 
region. 


on one specimen for each of the three alloys as shown 
in Fig. 10, and the specimens were bent to a 3/,-in. 


radius. The strain along the width, e,, was then deter- 
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mined by measuring the elongation between successive 
pairs of grid lines. 

The elongation data from these small width speci. 
mens are plotted on Fig. 11, with ¢,/e,%. as the ordinate 
and the specimen width as the abscissa. The ratio 
€,/€uit. is the strain in the width direction relative to the 
simple tensile strain in the bend direction. The latter 
strain, €,1,, is determined theoretically for the particular 
sheet thickness and bend radius according to the 
analysis presented in reference 4. 

Because of rather large scatter in the test data, it 
was thought desirable to group all data for the three 
alloys on one curve as shown in Fig. 11. On this 
curve, the elongation data obtained from the grid 
survey on an A-type specimen is also plotted to com- 
plete the range of widths tested. 

A curve was fitted to the test data as shown by the 
heavy line on Fig. 11. The equation of this curve ex. 
pressing the relationship between the lateral strain 
ratio and the width is 


éy/éut. = —0.02w-*/s (1) 


Theoretically, the limiting value of the ratio ¢,/e,, is 
0.5, which is Poisson’s ratio in the plastic range. 


Tensile Coupons 


Standard tensile coupons were cut cross-grain to 
duplicate the grain direction of the bend strips, viz., 
with the grain direction parallel to the bend line. Two 
coupons of each material were used. 

By use of grids photographed on the surface of the 
tensile coupons as shown in Fig. 12, the ultimate 
strain in various gage lengths was determined‘ for the 
three alloys. The strain results are plotted in Fig. 13 
and the data are used in the subsequent analysis. 

Along the heavy numbered lines of the grid (Fig. 12), 
the lateral contraction was determined for each tensile 
specimen, and applying the correction to the strain in 
the direction of applied load detailed in Appendix A, 
the ratio of lateral contraction e, to elongation €,),, was 
found. The results are plotted on the ordinate of 
Fig. 11 in the form of a distribution curve. The mean 
of the lateral contraction data was 0.35, indicating that 
in a tensile coupon the sides are restrained from com- 
plete contraction. Contraction would have been com- 
plete if the ratio were equal to 0.5, Poisson’s ratio in the 
plastic range. 


BEND WipTH ANALYSIS 


The strain in the y-direction which induces a stress 
in this direction is the difference between the lateral 
contraction arising by Poisson’s effect and the actual 
strain that does exist. This is illustrated on Fig. 14 
This stress, in turn, tends to produce a contraction i 
the x-direction so that the maximum tensile strail 
that can be attained is 
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€z = €ult. — vile, 6m (—veuie.) ] 
simplifying 
Ge = €uit.(1 r7 v?) ey (2) 


The relation between «¢,/e,;, and width has been 
empirically determined from lateral contraction meas- 
wrements on C-type specimens. From Eq. (1), 


€y = €uit.(—0.02w-*/*) 
Combining the latter with Eq. (2) leads to 
és = €uit.(L—v?) + veuie.(0.02w-'/*) 
Simplifying, 
és = €ut.[1 — v(v — 0.02w-*)] (3) 


It has been shown in reference 4 that the ultimate 
strain in the gage length of the bend can be predicted 
from ultimate strain-gage length data such as pre- 
snted in Fig. 13. These data are for a stress condi- 
tion approaching simple tension and, consequently, 
are inapplicable for predicting minimum bend radii 
for wide. bend specimens. 

By applying Eq. (3) to these data, however, the ulti- 
mate strain-gage length curves can be modified to ac- 
count for the effect of the width. For a given width, 
w, this equation reduced to the form 


Ss > éuit.(1 a K) 


where K is a constant for each width. Thus, the basic 
strain-gage length curve derived from tensile coupons 
can be corrected for the width effect, and the corrected 
curve can be used in solving for the minimum bend 
radius, 


Tensile Coupon Stress Condition 


Lateral contraction measurements on the tensile 
coupons showed the mean lateral contraction to be 
035 of the strain in the direction of applied load. 
Since the lateral contraction was less than 0.5, some 
restraint must have existed in the y-direction (Appendix 
A) similar to that in the width direction of a bend. 
This restraint in the width direction of the tensile 
specimen was probably caused by the geometry of the 
coupon in this plane. Lateral contraction measure- 
ments in the thickness direction indicated no restraint 
upon necking in this direction. 

In addition to the width direction stress on the tensile 
coupon, some restraint was imposed in the lateral direc- 
tions when necking occurred. However, since the 
tensile specimens were of -T condition material, necking 
is usually insignificant. Consequently, the restraint 
due to necking as outlined in reference 5 is small, and 
since scatter of the lateral contraction data was rela- 
lively large, the correction was neglected. 

On Fig. 11, the mean lateral contraction of 0.35 for 
the tensile specimens corresponded to a bend width 
0.1 in. Thus, the bend width for which minimum 
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Fic. 13. Ultimate strain in various gage lengths. Data obtained 
from gridded tensile coupons. 
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Fic. 14. State of strain in a wide bend. 


bend radii are theoretically determined was estab- 
lished. 


Summary of Analysis and Theoretical Results 


To determine the minimum bend radius for a bend of 
a given width, the following calculations are performed 
with the ultimate strain-gage length data of Fig. 13: 
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(a) Since these tests were conducted on clad TABLE 3 
material, correct the curve for cladding thickness as Cladding and Width Corrections for Ultimate Strain-Gage 
outlined in Appendix B. Length Data 
€ult., = €y/(1 — 2c) (A-3) w = 0.1 in. 
; Gage fult.. w=lin. w = 100in, 
(b) Correct the curve for tensile specimen lateral Length, eu. Eq. (A — fult.. fult., 
restriction to a bend width in which anticlastic curva- In. (Fig. 13) 3a, b, c) Eq. (6b) Eq. (6c) 
ture is complete and the forming stresses are simple Alclad 24S-T 
tension stresses, Eq. (3). 0.02 0.360 0.379 0.312 0.307 
0.06 0.330 0.347 0.285 0.281 
Eult.,=9 = €uit../[1 — 0.5(0.50 — 0.35) 0.10 0.306 0.322 0.264 0.261 
on “ 0.14 0.301 0.317 0.260 0.257 
fult.y=o = éult.,/0.925 (4) 0.22 0.274 0.288 0.236 0.233 
5 2 
It is convenient to combine Eqs. (A-3) and (4). 3 eo — 4 po fi ‘a 
Eult.y mo = €uit,/0.925(1 — 2c) (5) Alclad 75S-T 
. . i ; : 0.334 .27 .27 
(c) Correct for the width effect (including correc- ; a 5 moe 0 pour ; sas : it 
tions (a) and (b)), 0.10 0.219 0.238 0.196 0.193 
/ 0.14 0.194 0.211 0.173 0.171 
€xw = [ewte,/0.925(1 — 2c)][1 — v(v — 0.02/w**)] (6) 0.22 0.170 0.185 0.152 0. 150 
or 0.30 0.149 0.162 0.133 0.131 
0.50 0.130 0.141 0.116 0.114 
few = (Exit.-/0.925) [1 — v(v — 0.02/w*/*) ] (6a) R301-T 
From Fig. 13, the minimum ultimate strain in each 0.02 0.344 0.382 0.322 0.318 
0.06 0.285 0.317 0.260 0.257 
gage length for the three alloys was corrected for the 0.10 0.256 0.284 0.233 0.230 
cladding thickness following the method of Appendix B. 0.14 0.238 0.264 0.217 0.214 
The cladding-thickness ratios in the 0.064-in. gage are: 0.22 0.197 0.219 0.180 0.177 
= 0.30 0.170 0.189 0.155 0.153 
Alclad 24S-T c = 0.025 0.50 0.135 0.150 0.123 0.121 
Alclad 75S-T 0.040 Ky 
R301-T 0.050 
TABLE 4 
Using Eq. (A-3), for Theoretical Minimum Bend Radii for Various Widths 
Alclad 24S-T, €utt.. = €4,./0.95 (A-3a) w, In. R;, In. Rint 
bh Alclad 24S-T 0.1 0.092 1.44 th 
Alclad 75S-T, €utt., = €44,/0.92 (A-3b) 1 0.101 1.57 to 
100 0.121 1.90 
o € = / = at 
ROUT Qrttre =. Sein. 0-00 fAate) Alclad 75S-T 0.1 0.182 2.85 be 
The calculations for each material are presented in 1 0.252 3.93 fe 
‘ 100 0.272 4.24 
Table 3. : th 
Since the tensile test data without correction for Ret ~g Z joe oa is 
width applied directly to a bend of 0.1-in. width, these 100 0.295 3 59 
data were used directly in calculating the minimum in 
bend cap for a —" of 0.1 in. after-the cladding The results presented in Table 4 are plotted on Fig. lit 
ee ee 16 to indicate the a of 
. : greement between theory and the 
— om a. of 1 _ ees in., oF og data were experimental data of Table 1. be 
reduced by using Eq. (6a). For w = 1 in.: Because of the small number of specimens used im st 
€xy a1 = (€ult../0.925)[1 — 0.5(0.5 — 0.02/1)] the testing program, the results presented herein must Ta 
p = 0.8026 6b be considered qualitative in nature. Definite trends 
ta las iiace — (6b) were indicated by the investigation, but several parame- at 
for w = 100 in.: ters that may be of some importance were neglected. ar 
; s ntietie st 
Sic cos = 0.812 On, (6c) It is felt, however, that the results are sufficiently i “ 
dicative to obviate further time-consuming investiga 
These calculations are also presented in Table 3. tion, although the subject is of sufficient academic he 
These ultimate strain-gage length data are plotted interest to warrant considerable study.* wi 
on Fig. 15. The minimum bend radii for the 0.064-in. ~{— PPO tN eT ETE Toe si ial 
. . . n intensive investigation oO eentire bending probiem 16 
gage materials hatecie determined by the ROO aphic current project of the Office of Scientific Research and Develop- 
solution detailed in reference 4, the results of which are ment being conducted at the Case School of Applied Science co 
presented in Table 4. under Dr. George Sachs. Al 
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Fic. 15. Ultimate strain data in various gage lengths corrected 
for bi-axiality of stress existing in bends of several widths. 


One specimen of each type of material was tested at 
the various bend radii, and thus no attempt was made 
to account for variations in test results which may 
arise from different sheet lots. For a given width of 
bend, this anticipated variation may result in a dif- 
ferent minimum bend radius, but it is expected that 
the indicated trend of increase in bend radius with width 
is unaffected. 

The effect of sheet thickness upon the results was not 
investigated, but this parameter is thought to have 
little influence on the width effect. The thickness is 
of fundamental importance in establishing minimum 
bend radii, but since the width effect arises from the 
stress condition imposed in attaining the minimum 
tadius, the thickness should have little effect per se. 

Because of the small number of specimens used, no 
attempt is made to find any difference in width effect 
among the three alloys tested. However, from stress- 
strain considerations in the plastic range, it appears 
that the material with the smallest degree of strain 
hardening may be the most seriously affected by the 
width effect. 

The agreement between theory and test data (Fig. 
16) is good for Alclad 24S-T, but the theory predicts 
considerably conservative bend radii for R301-T and 
Alclad 75S-T. This may be partially due to errors in 


test data. Such errors would reflect in the theoretical 
results, since it is necessary to establish certain em- 
pirical relationships as a basis for the analysis. 

The theory gives a sharp rise in minimum bend radius 
in small widths and levels off thereafter, as would be 
expected from the trend of lateral contraction data 
with width. Possibly the bend width of 0.1 in., for 
which correlation with tensile coupon data has been 
obtained, is in error. Even a small change in this 
region would cause a large decrease in minimum bend 
radii predictions in the large widths and yield better 
agreement. 


CONCLUSIONS 


The three aluminum alloys tested, R301-T, Alclad 
24S-T, and Alclad 75S-T, indicate a definite increase 
in minimum bend radius of between '/32 and 1/, in. 
for a range of bend widths between 1 and 100 in. This 
effect of width upon minimum bend radius is not too 
serious, but it may be of some consequence when bend 
radii based on small width specimens are used on long 
sections, such as stringers or longerons. 

The practice often used is to establish minimum 
bend radii from results of tests on small width speci- 
mens because of the ease with which such specimens 
can be handled. As a result of work herein, it is sug- 
gested that, when minimum bend radii are established, 
specimens of large widths be included in the test pro- 
gram. If the minimum bend radius varies consider- 
ably over the width range, it may then be desirable to 
establish a set of bend radius minimums for smaller 
and larger widths. 


(Continued on page 170) 
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Fic. 16. Agreement between theory and test data on variation 
of minimum R/t with width. 
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APPENDIX A 

Reduction of Tensile Coupon Elongation Data 

In a tensile coupon, a restriction exists in the direc- 
tion normal to the applied load similar to that in a wide 
bend. The state of strain in the tensile coupon is con- 
sidered identical to that shown in Fig. 14. Hence, the 
strain measured in the x-direction is decreased from 
the ultimate value 


€e = €utt. — Vley — (— veut.)] 
To find the strain under simple tension, 
éult. = (€2 + vey)/(1 — v’) 
or to determine the ratio ¢,/€,1. 


€y/€utt. = e(1 “ar v*)/(€ + Vey) 
Since vy = 0.5, 


€y/Eutt. oad 0.75¢,/ [ee + (€,/2) ] (A-1) 


APPENDIX B 
Correction of Ultimate Strain Data for Cladding Thickness 


Cladding material is a low yield strength alloy 
capable of sustaining much larger ultimate strains 
than the base metal. Consequently, it is assumed 
that when the base material reaches the ultimate strain 
it will fracture underneath the cladding, as shown in 
Fig. 17. With continued application of load complete 
fracture will result. 

From reference 4, the ultimate strain in the gage 
length of the bend is 


éut, = K/(R/tets.) 


The effective thickness subjected to the ultimate 
strain is less than the actual sheet gage, 








Fic. 17. Bend parameters for material of clad aluminum alloy, 


fea, = (1 — 2c)t (A-2) 


where c = ratio of cladding thickness per face to sheet 
thickness. 
Therefore, for clad sheet, 


fur. = [K/(R/t)](1 — 2c) 
letting €ult., = €,/(1 — 2c), 


eult.. = K/(R/t) (A-3) 
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Atomic Energy 


HSUE-SHEN TSIEN* 
California Institute of Technology 


INTRODUCTION 


. pw SPECTACULAR RESULTS achieved by the use of 
atomic bombs toward the end of World War II 
have greatly stimulated the interest in the possibilities 
of atomic energy in other fields of engineering applica- 
tions. The very fact that the energy release by nuclear 
reactions is approximately a million times that of the 
more conventional chemical reactions seems to stagger 
one’s imagination. Although there are wide differences 
in opinion as to the time required to develop this branch 
of new-found knowledge to practical power-plant engi- 
neering, there is no diversity in the belief that a new 
era of technological evolution has begun. To keep 
abreast with the expected rapid advancement in this 
field, the engineers should have a clear understanding 
of the fundamental concepts involved as advocated by 
von K4rman in a recent article.! 

The aeronautical engineers have a further interest in 
the matter due to the belief that the aeronautical power 
plant will probably be the first prime mover to utilize 
the atomic energy. This belief is based upon the 
following fact: For stationary power plants, the cri- 
terion for economic operation is the cost to produce one 
kilowatt-hour and thus the cost of the fuel, be it chemi- 
cal or atomic, is of prime importance instead of the 
weight of the fuel. For automotive applications and 
especially aeronautical applications, the weight of the 
fuel is extremely important. In case of the projected 
supersonic flight at extreme speeds, it seems that only 
the use of atomic energy will make such flight economi- 
cally feasible by reducing the fuel load and increasing 
the pay load. 

This article is a digest of the basic concepts in this 
branch of knowledge as an introductory study of this 
new field. The enormous energy release of nuclear 
reaction is explained through the binding energy of the 
atomic nuclei, with the energy generation in the stars 
and the atomic bombs as typical examples. Frequent 
comparison between molecular reaction and nuclear 
reaction will be made to guide the exploration of an un- 
familiar land. 


THE EQUIVALENCE OF ENERGY AND Mass 


The astonishing amount of energy released by the 
nuclear reactions can be easily understood through the 
still greater energy potential in the atom itself repre- 
sented by its mass. Therefore, one of the fundamental 


Received September 29, 1945. 
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concepts of the atomic energy is the concept of the 
equivalence of energy and mass. This concept was 
first announced by Einstein as a conclusion of his 
special theory of relativity. This theory deals with the 
comparison of measurements made by observers who 
are assumed to be in unaccelerated relative motion and 
in a region of space and time where the action of gravi- 
tation can be neglected. The basis of the theory can be 
conveniently regarded as the two following fundamental 
postulates :? 

(a) It is impossible to measure or detect the un- 
accelerated translatory motion of a system through free 
space or through any ether-like medium that might be 
assumed to pervade it. 

(b) The velocity of light in free space is the same for 
all observers, independent of the relative velocity of the 
source of light and the observers. 

The first postulate is derived from the failure to de- 
tect the ether drift of the motion of the earth through 
the space, and the second postulate can be derived from 
both astronomical evidences and laboratory experi- 
ments. Therefore, although Einstein’s theory is often 
disseminated by popular science writers as an abstruse 
mathematical theory, it really is solidly founded on 
observational facts. In this sense, the relativity theory 
is just as empirical as a stress-strain curve drawn 
through test points obtained on a steel specimen. 

The first postulate of special relativity is already con- 
tained in the Newtonian dynamics and can be accepted 
as part of everyday experience. The second postulate, 
although familiar in wave theory of light, is somewhat 
beyond the daily experience. In fact, one would gen- 
erally expect the velocity of propagation and the 
velocity of the source to be additive. The combination 
of these two postulates of widely different natures 
would certainly lead to entirely new results. This 
expectation is, of course, fully substantiated. One of 
the results of deduction is that, if the mass of a particle 
is my at rest, then its mass m at velocity V is 


m = m/V1 — (V2/c?) (1) 


where c¢ is the constant velocity of light. Since the 
velocity of light is about a million times the velocity of 
sound at ordinary conditions, for all velocities realized 
in the terrestrial macroscopic world the mass m is 
practically the rest mass m%. The kinetic energy calcu- 
lated by the relativity theory is 


E = c?(m — mo) (2) 
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For small velocities—i.e., V < c—the preceding equa- 
tion reduces to the familiar form E = '/gm)V". 

This close relationship between energy and mass, as 
well as other evidences, prompted Einstein to take a 
revolutionary step. He says that there is no difference 
between energy and mass. One is just a form of the 
other. They are closely associated with each other. 
The energy E associated with any mass m is E = me’. 
Any energy E has a mass m = E/c? associated with it. 
In other words, a complete equivalence of energy and 
mass is achieved here. The old concept of separate 
conservation of energy and mass is thus merged into 
one single concept of conservation of “‘energy-mass.”’ 

Since the velocity of light is generally given as 
approximately 3 X 10" cm. per sec., 1- lb. of mass is 
equivalent to 


: 1 \/_1\(8 x 10) 

= +. (2 V2.4 2h) ; 
(3) , a cor” X (3.280) 
3.86 X 10" B.tu. 


1.1380 X 10” kwh. 


(3) 


The problem is then how to transform the mass into 
energy or how to release the energy intrinsically 
bounded with the mass. 


ATOMIC STRUCTURE 


The problem is similar to that of releasing the chemi- 
cal energy, intrinsically bounded with a mixture of 
hydrogen and oxygen. Until a device for igniting this 
gaseous mixture is discovered, the chemical reaction, 
although known to be possible, cannot be started and no 
chemical energy of this mixture will be released. Thus, 
the central problems of converting mass into energy are 
first to discover the type of reactions which can perform 
this mass-energy transformation and then to find the 
means of “‘igniting’’ or starting such reactions. The 
only known reactions that will transform mass into 
energy are nuclear reaction or the reactions with the 
atomic nuclei as the participants. Before entering into 
the discussion of the nuclear reaction, it is thus profit- 
able to study the structure of atoms first. 

The atom is composed of the central core or the 
nucleus and the outer collection of electrons. The 
mass of the atom is mostly concentrated in the nucleus. 
The electrons are extremely light. The mass of an 
electron at rest is 9.035 X 10-% g., and in the physi- 
cal atomic weight* scale it is 0.000549. The total mass 
of the electrons is only several thousandths that of the 
nucleus. The charge of.the nucleus is positive and that 
of the electrons is negative. For a normal atom the 
positive charge of the nucleus is just equal to the nega- 
tive charge of the electrons and thus gives a neutral 
atom. 

*In the physical atomic weight scale, the atomic weight of 
O*, or the oxygen atom with 8 protons and 8 neutrons, is ex- 
actly 16. Because of a slight amount of heavier atoms in natural 


oxygen, the usual atomic weight scale or the chemical atomic 
weight scale is slightly larger (by about 1 part in 4,000). 





TABLE 1 
Electronic Quantum States 








Orbitual Quantum States- 


—, 


Shell _ 
1s(1) 
2s(1) 2p(3) , 
3s(1) 3p(8) 3d(5) 

4s(1) 4p(3) 4d(5) 4f(7) 

5s(1) 5p(3) 5d(5) 5f(7) 5g(9) 

6s(1) 6p(3) 6d(5) 6f(7) 6g(9) 6h(11) 


7s(1) 7p(3) 7d(5) 7f(7) 7g(9) 7h(11) 7i(13) 
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Fic. 1. Energy levels of electronic quantum states. 3d 
states are first occupied by the iron transition elements. 4d 
states are first occupied by the palladium transition elements. 
5d states are first occupied by the platinum transition elements. 
4f states are first occupied by the rare earth elements. 


The electrons are arranged in ‘‘shells’’ around the 
nucleus. In terms of atomic scales there is an enormous 
distance between the nucleus and the first shell. This 
is indicated by the fact that the diameter of the nucleus 
is less than '/,9,000 of the diameter of the atom. The 
shells are named as K, L, M, N,O, P,Q. The states of 
motion the electron might occupy are designated by 
three-positional or orbitual quantum numbers and a 
“nonclassical” spin quantum number. For each orbit- 
ual quantum number there are two spin quantum num- 
bers. These quantum numbers, together with the 
energy of the corresponding states of motion, can be 
determined either experimentally by measuring the 
spectrum of the element or theoretically by calculation 
with the rules of quantum mechanics assuming Cou- 
lomb forces between the electrons and the nucleus. In 
fact, the excellent agreement between the theory and 
the experiment is one of the major triumphs of modern 
physics. It is also found that no two electrons can 
occupy a state with the same four quantum numbers, 
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according to the “Pauli exclusion principle.’ The 
orbitual states are listed in Table 1. The number in 
the parentheses is the number of orbitual states. Since 
there are two spin numbers with each orbitual state, 
the number of electron positions is twice the number 
shown in the parentheses. 

With increasing atomic weight, the nucleus gets 
heavier, and the charge of the nucleus is increased with 
increase in the number of electrons attached to it. In 
other words, the shells of electrons gradually fill up. 
The states or “‘positions’”’ with the loweSt potential is 
filled up first. The potential ladder is shown in Fig. 1. 
Hydrogen uses one of the 1s electron positions. Helium 
uses both of the 1s electron positions and completes the 
K shell. Hence, helium shows little tendency to draw 
electrons of other atoms to form chemical compounds 
and is an inert gas. Neon completes the L shell and be- 
comes an inert gas. Now, the chemical property of 
the atom is mainly determined by the outmost shell. 
Since the Fe, Pd, Pt, transition elements, and the rare 
earth elements are made by filling up the inner shells 
with 3d, 4d, 5d, and 4f positions, respectively, and with- 
out changing the outer shell, their respective similarity 
in chemical property is expected.* For the same rea- 
son, the similarity between the elements U, Np,* Pu,* 
indicate that they may be the elements made by filling 
up the 5f positions in the O shell with the two electrons 
in the Q shell unchanged. These elements then might 
constitute the beginning of a new series similar to the 
rare earth elements.‘ 

The nucleus is now generally considered as composed 
of neutrons and protons. They have about the equal 
mass, with the neutron being slightly heavier, and are 
approximately 1,840 times heavier than electrons. In 
the physical atomic weight scale, a neutron is 1.00893 
and a proton is 1.00757. Similar to the electron, their 
states are determined by the orbitual numbers and the 
two spin numbers. Since only the proton has a unit 
positive charge, the number of protons in the nucleus 
must be equal to the number of electrons. This number 
is the atomic number Z. Hydrogen has the atomic 
number 1 and uranium 92. The total number of 
neutrons and protons in the nucleus is the mass number 
A. Thus VN = A — Z is the number of neutrons in the 
nucleus. Since the mass of neutrons or protons is 
approximately unity in the scale of atomic weight, A is 
also roughly the atomic weight. Therefore, the nu- 
cleus of the hydrogen is a proton. The nucleus of 
deuterium is deuteron, or a proton anda neutron. The 
nucleus of helium is usually called the a-particle and is 
composed of two neutrons and two protons. Elements 
with equal number Z are called ‘‘isotopes.’’ Thus, 
hydrogen and deuterium are isotopes. Similarly, the 
uranium isotopes U-235 and U-238 are composed of 92 
protons and 133 and 136 neutrons, respectively. 





*Np and Pu stand for neptunium and plutonium, two new . 


elements beyond uranium with 93 and 94 electrons, respectively 
(atomic numbers 93 and 94). 


The neutrons and protons in the nucleus are not a 
loose assemblage but are tightly bound together. In 
fact, their bounds are much stronger than the bounds 
between the nucleus and the electrons in the outer 
shells. For instance, while it takes only 13.53 volts to 
remove the electron from the hydrogen atom—.e., a 
binding energy of 13.53 ev.{—the binding energy be- 
tween a proton and a neutron is approximately 2.15 X . 
10° ev. The average binding energy per particle for 
heavier elements is even higher, being approximately 
8.5 X 10% ev. The ratio for these two types of bounds 
is thus roughly 10°. The ordinary chemical energy re- 
leased by molecular reactions is linked with the rear- 
rangement of electrons in the molecules. For instance, 
the electrons in hydrogen molecules and in oxygen 
molecules are bound to the hydrogen nuclei and the 
oxygen nuclei, respectively. By rearranging the elec- 
trons so that they are bound to both the hydrogen 
nucleus and the oxygen nucleus at the same time, water 
molecules are formed. ‘The closer binding between the 
electrons and the nuclei in the water molecule corre- 
sponds to the chemical energy released during the com- 
bustion of hydrogen with oxygen. Therefore, if 
nuclear reaction involving rearrangements of the 
particles in nuclear structure is realized, then the associ- 
ated energy released must be of the order of the nuclear 
binding energy—i.e., approximately a million times that 
of chemical or molecular reaction. This difference in 
the order of magnitude between the nuclear reactions 
and the molecular reactions is basic and will occur fre- 
quently in the subsequent discussions. 


NUCLEAR REACTIONS 


The nuclear reactions known up until 1939 are reac- 
tions produced by the bombardment of the atoms with 
various kinds of projectiles, such as a-particle, proton, 
neutron, and deuteron. The result of this bombard- 
ment is either absorption of the bombarding particle 
with emission of short wave-length electromagnetic 
waves or y-ray, or disintegration of the atom with 
emission of other particles. The reactions can be 
classified as in Table 2. For instance, in 1932 J. D. 
Cockcroft and E. T. S. Walton of Rutherford’s labora- 
tory, bombarded a target of lithium with protons of 700 
kv. energy and found that a-particles were ejected from 
the target as a result of the bombardment. The 
nuclear reaction that occurred can be written sym- 
bolically as 

3Li’ + ,H' — ,He* + 2Het (4) 


where the subscript represents the atomic number Z 
and the superscript the mass number A. Ina chemical 
equation, the number of atoms of any element entering 
into the reaction must be equal to the number of atoms 


t ev. is the abbreviation for electron-volt—i.e., the energy 
spent in moving a unit electron charge through a potential of 1 
volt. lev. = 4.45 X 10°-*kwh. 
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of that element in the product. In a nuclear reaction 
as shown in Eq. (4), the number of the nuclear particles, 
protons, and neutrons must be the same on both sides. 
Thus the sum of subscripts is four and the sum of super- 
scripts eight on each side. 

Neither mass nor energy has been included in this 
equation. In general, the sum of the masses of the 
incident proton and the lithium atom will not be pre- 
cisely the same as the sum of the masses of the a- 
particles produced. According to Einstein’s law of 
conservation of energy-mass, the sums of mass and 
energy taken together should be the same before and 
after the reaction. The masses were known from mass 
spectra. On the left (Li? + H") they totaled 8.0241; 
on the right (2He*), 8.0056, so that 0.0185 unit of mass 
had disappeared in the reaction. The experimentally 
determined energies of the a-particles were approxi- 
mately 8.5 Mev.* each, a figure compared to which the 
kinetic energy of the incident proton could be neglected. 
Thus, 0.0185 units of mass are 3.07 X 10-* g. or 
17.2 Mev. according to Einstein’s equation. Therefore 
the experimental results have completely proved Ein- 
stein’s mass energy equivalence. 

As stated, the protons for the lithium reaction are 
accelerated by 700 kv. or have 0.7 Mev. energy per 
particle. Such an amount of kinetic energy gives the 
proton an appreciable probability of penetrating into a 
lithium nucleus and causing a reaction. However, out 
of a million accelerated protons there still is only one 
that will make a successful hit; the others will be 
slowed down by collisions with atoms without producing 
the desired reaction, and their kinetic energy will be 
lost by transforming into heat. Thus 0.7 X 10° Mev. 
energy is required to accelerate a million protons to 
produce one reaction giving 17 Mev. energy. Hence, 
40,000 times more energy is spent than recovered. It 
is obviously impracticable to use this nuclear process to 
produce energy. ‘The situation could be compared to 
that of a hypothetical molecular reaction where an in- 
efficient igniting device has to be used for every pair of 
molecules. Then the chemical reaction, even if it were 
an energetic one, cannot be practically utilized for 
energy production. 

There is a special class of nuclear reactions which does 
not involve incident particles. This means that the 
nuclei automatically disintegrate with the emission of 
particles. This phenomenon is called the radioactivity 
and is discovered actually earlier than the usual nuclear 
reactions described above. The emitted particles could 
be the a-particles, the protons, the electrons, and the 
positrons, positive electrons, as shown in Table 2. 
The emission of electrons and positrons are found, how- 
ever, to give continuous spectra of the kinetic energy of 
the emitted particles, while the energy states of the 
nuclei are found to be discrete. This is in contradiction 
to the law of conservation energy-mass, which is 

* 1 Mev. is one million electron-volts or 10* ev. and is equiva- 
lent to 4.45 XK 10-” kwh. 





TABLE 2 
Types of Nuclear Reaction 








Incident Emitted 
Particle Particle Z Becomes A Becomes 
a Proton Z+1 Z+3 
a Neutron Z+2 A+3 
= Proton Z+1 A+1 
= 2 Proton a Z-1 A-3 
Ay E Proton Neutron Z+1 A 
Bt Proton Deuteron Z A-1 
4 s Deuteron Proton Z A+1 
‘3 bs Deuteron Neutron Z+1 A+1 
gm Deuteron a Z-1 Z-2 
x Neutron basis Z A+1 
Neutron Proton Z-1 A 
Neutron a Z-2 A-3 
a> a Z-2 A-4 
3 = Proton Z-1 A-1 
2 : Electron Z+1 A 
- Positron Z-1 A 





observed to hold in all other cases. To preserve this 
law, a neutral particle with spin but of practically zero 
mass called neutrino is assumed by W. Pauli. The 
neutrino is, so far, not directly detected. However, 
because of the assumed existence of neutrinos, part of 
the energy could be transferred to them, and the con- 
tinuous spectra of electron energy in radioactivity is ex- 
plained. The kinetic energy of the neutrinos is “‘lost” 
because it cannot be detected. 


NUCLEAR STRUCTURE—BINDING ENERGY 


An entirely new type of nuclear reactions is the 
breaking up of a heavy nucleus into two approximately 
equal fragments. This is called the nuclear fission 
reaction. It was discovered by L. Meitner, O. Hahn, 
and F. Strassmann toward the end of 1938. To under- 
stand the reaction one must look into the structure of 
the nucleus in more detail.® 

In studying nuclear structure one is, however, faced 
with the difficulty of not knowing the character of the 
forces between the nuclear particles. The desired in- 
formation must be deduced from the scanty experi- 
mental data. First, by plotting the atomic number Z 
against the mass number A for all known elements, 
there is an obvious relation between these two quanti- 
ties. This relation is, of course, not precise because of 
the variations by isotopes. For instance, with the 
same atomic number 1 for hydrogen and deuterium, 
the mass numbers are 1 and 2, respectively. However, 
for light elements, the value of A is approximately 
twice that of Z. This is exactly true for »C’, ;N“, 
30"*, etc. This means that the stable light nuclei are 
composed of equal numbers of protons and neutrons. 
Then, the binding force between the nuclear particles 
must be the greatest between the proton and the neu- 
tron. If the force between protons were stronger, the 
stable nucleus would contain more protons than neu- 
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If the force between the neutrons were stronger, 
How- 


trons. 
the stable nucleus would contain more neutrons. 
ever, the fact that the force between the neutron and 
the proton is the strongest does not exclude the possi- 
bility of forces between neutrons themselves and protons 


themselves. But if such binding force between the like 
particles exists, they cannot be widely different, be- 
cause, if the binding force between a pair of neutrons 
were much higher than that between a pair of protons, 
then the stable nucleus would again have a preponder- 
ance of neutron instead of the observed equality in 
numbers of these two kinds of particles. Of course, 
there will be Coulomb repulsion between two protons. 
But this repulsion can be shown to be small in compari- 
In fact, it is 
generally assumed that aside from the Coulomb force, 
the force between a pair of neutrons is equal to the force 
between a pair of protons. 

The a-particle, or the nucleus of the helium atom, 
consists of two protons and two neutrons. The binding 
energy per particle of this nucleus is found to be much 
higher than the deuteron. Furthermore, a nucleus of 
five particles, by adding either a proton or a neutron to 
the a-particle, is never observed. This means that a 
nucleus of five particles is highly unstable and its life- 
time is so short that its detection by experiment is 
highly unlikely. This great stability of the helium 
nucleus and the difficulty of adding another particle to 
it show the saturation of the binding forces in helium 
nucleus. The situation is similar to that of the helium 
atom with the two electrons in the same orbitual quan- 
tum state but different spin states. The resultant 
closed K-shell structure refuses to combine with addi- 
tional electrons from other atoms to form molecules. 
Thus, similarly, the two neutrons must be in the same 
orbitual quantum state with different spin states, and 
the two protons must be in the same orbitual quantum 
state with different spin states. But the large binding 
energy of the helium nucleus requires that there is bind- 
ing force between every pair of the particles in this 
nucleus. Hence, the binding force between a proton 
and a neutron and the binding force between the like 
particles cannot depend to any considerable extent 
upon the relative spin directions of the interacting 
particles. 

If there is binding energy between every pair of inter- 
acting particles, then, since the number of pairs in a 
nucleus is ('/2)A(A — 1), the total binding energy of a 
nucleus must be proportional to the square of A. 
However, experimental data show a practically linear 
relation between the binding energy and the mass num- 
ber A. This is the same as the relation between 
molecular binding energy of a liquid drop containing 
many molecules. It is known that the force between 
the molecules in liquid or solid shows saturation. In 
other words, a molecule is attracted only to the other 
molecules in its immediate neighborhood. Thus, the 
intermolecular forces can be said to be short ranged, 


speaking in terms of the dimension of the drop. This 
character of the binding force greatly limits the number 
of effectively interacting pairs of particles and makes 
the total binding energy proportional to the number of 
particles. The nuclear forces then must show a similar 
saturation character. Here, a neutron or a proton only 
interacts strongest toward another particle of nearest 
orbitual quantum states. As the number of the parti-— 
cles increases with increasing mass number A, the num- 
ber of particles having nearly a given orbitual quantum 
state does not increase, since only two particles of the 
same kind can occupy the same orbitual quantum 
states, according to Pauli’s exclusion principle. Then 
the nuclear forces can also be said to be short-ranged in 
terms ,of the dimension of the nucleus. Thus, the 
number of effectively interacting pairs is limited in a 
similar manner as in the liquid drop, and the total 
binding energy is only proportional to the total number 
of particles instead of proportional to the square of the 
number of particles. 

This liquid drop analogy is also useful in two other 
connections. First, since the volume of the liquid drop 
is proportional to the number of molecules, the volume 
of a nucleus must also be proportional to the number of 
nuclear particles, since they are similarly held together. 
This means that the radius of a nucleus is proportional 
to A’*, This is found to be true experimentally. 
Secondly, a liquid drop exhibits the well-known surface 
tension phenomenon. This is due to the existence of a 
free surface where the molecules are only attracted by 
the molecules within the drop. Therefore, for the 
molecules on the surface, only half of their binding 
force is satisfied, and their contribution to the total 
binding energy of the drop is not so large as the mole- 
cules within the drop. A similar situation exists for the 
nucleus, where the protons or neutrons on the surface of 
the nucleus contribute much less binding energy than 
the particles within the nucleus. For a heavy nucleus, 
the ratio of the particle at the surface to the total num- 
ber of particles is reduced and the binding energy per 
particle should be increased. In other words, the sur- 
face tension effect indicates an increase in binding 
energy per particle with increasing mass number A. 

So far, the effect of Coulomb repulsion between pro- 
tons is neglected. A simple calculation shows that the 
average energy due to this effect is '/, Mev. per pair of 
protons. The average binding energy due to nuclear 
forces is 8'/2 Mev. per particle. Thus the Coulomb 
repulsion is negligible if compared on the basis of a 
single pair of protons. However, the nuclear forces 
show saturation, and the binding energy due to nuclear 
forces is approximately proportional to A. The num- 
ber of protons in the nucleus and, thus, the total proton 
charge are also approximately proportional to A. The 
radius of the nucleus is proportional to A’. Hence, 
the total repulsion energy is proportional to A XK A/A ‘v/s 
or A’*. This difference of the Coulomb repulsion 
energy from the energy due to nuclear forces is, of 
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course, the result of nonsaturation character of the 
Coulomb forces. The ratio of these two energies is thus 
proportional to A”’*. Therefore, for a large nucleus, 
the Coulomb repulsion energy is important in spite of 
its smallness for one pair of protons. Then the total 
binding energy of the nucleus will be increased and the 
nucleus made more stable by substituting some of the 
protons by neutrons. This is actually found to be the 
case: The uranium nucleus has a neutron-proton ratio 
N/Z equal to 1.6, while for light nuclei this ratio is 1 as 
stated before. Furthermore, the effect of the Coulomb 
repulsion will also make the total binding energy in- 
crease less rapidly as the increase in number of particles 
A. This effect is thus opposite to the surface tension 
effect. However, for small nuclei, the Coulomb repul- 
sion is unimportant. Hence, the binding energy per 
particle will increase with A in this range. For large 
nuclei, or large A, the binding energy per particle will 
again decrease as the surface tension effect is overcome 
by the Coulomb repulsion. 

To actually calculate the nuclear binding energy 
from the experimental data, one can proceed as follows: 
An atom of atomic number Z and mass number A con- 
tains Z protons and Z electrons together with N = 
A — Z neutrons. Because of the negligible binding 
energy between electrons and the nucleus, the mass of Z 
protons and Z electrons is equal to Z times the mass of 
hydrogen atom. The atomic mass of hydrogen and the 
mass of neutrons are 1.00813 and 1.00893, respectively, 
in the physical scale. Therefore, the mass of compo- 
nents that go into the structure of the atom is 1.008132 
+ 1.00893N. If the actually measured atomic weight 
of the element is /, then the mass defect is 1.00813Z + 
1.00893N — M. This mass defect must represent the 
binding energy of the atom according to Einstein’s law 
of equivalence of energy and mass. Since the binding 
energy of the electrons to the nucleus is negligible, the 
mass defect is practically the binding energy of the 
nucleus. Generally, a different parameter is intro- 
duced. This parameter is called the packing fraction f 
defined as 





f = (M — A)/A (6) 
The binding energy per particle is then 
1.008132 + 1.00893N — M _ ~f +4 0.008132 r 
A A 
0.00893 =-f+ 0.01706 — 0.00893 (7) 


By neglecting the small variation of the ratio Z/A, 
the binding energy per particle is given by the negative 
of the packing fraction f. Therefore, for increasing 
binding energy per particle, the packing fraction. in- 
creases. According to the discussion given in the 
previous paragraph, the binding energy per particle 
must show a maximum for medium heavy nuclei. 
Thus, the packing fraction must show a minimum for 
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the medium values of the mass number A. This is seen 
in Fig. 2.5 

By examining the graph for the packing fractions 
(Fig. 2), it is seen that very large energy can be re- 
leased by either of the two methods: manufacture of 
the medium heavy elements out of light elements, or 
breakdown of an extremely heavy element into the 
medium heavy elements. In other words, by going 
from the two ends of the periodic table toward the 
middle, closer packing of the nucleus can be achieved 
with the release of the enormous binding energy. 
Nature provides an exainple for each case. The manu- 
facture of helium from hydrogen is of the first type and 
is carried out in the center of the main sequence stars 
such asthesun. The fission of uranium is of the second 
type. This reaction is used in the atomic bomb. 
These two typical examples will be examined in more 
detail presently. 


ENERGY PRODUCTION IN THE STARS 


In the laboratory nuclear transmutations are pro- 
duced by accelerating a few particles to high kinetic 
energy. The resultant poor efficiency is explained 
under ‘“‘Nuclear Reactions.’’ In stars, different condi- 
tions prevail: Because of the high temperature in the 
interior of stars, all protons in the stars have high 
kinetic energies. Moreover, they are not slowed down 
by collisions with other atoms because all atoms have 
equally high energies. Therefore, unlike the laboratory 
nuclear reaction, no energy need be spent in accelerating 
the bombarding particles such as protons, and any 
energy released by the nuclear reaction is a net gait. 
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On the other hand, the kinetic energies of the nuclei in 
stars are small compared with those used in the labora- 
tory. It is true that the temperature is about 20,000,- 
900°C. at the center of sun, and the atomic nuclei are 
therefore traveling with much greater speed than the 
molecules in the terrestrial atmosphere. However, 
their average kinetic energy is still only about 0.003 
Mev.—200 times smaller than the energies commonly 
used in the laboratory. Therefore, their probability of 
penetrating into the nucleus of the bombarded atom is 
much smaller. This is, however, remedied by the pres- 
ence of enormous numbers of bombarding protons, and 
the net reaction rate is not small. 

Through a careful examination of the probable reac- 
tion, Bethe’ picked out a series of the nuclear reactions 
as follows: 


CY! + HI N® + (8) 
™N¥¥ — 6C + 3° (9) 
(C+ sH* > N+ 7 (10) 
Nt + iH! > 408 + (11) 
gO — N® + 1 (12) 
iN + iH! > .C? + :He! (13) 


Here y represents the emission of y-rays and ,e° repre- 
sents the positron that is emitted together with neu- 
trinos, as explained previously. This series of reactions 
is remarkable in that it reproduces the carbon atom 
after six reactions. This is important, since the content 
of carbon in the stars is not high while hydrogen is 
plentiful. The net result is the manufacture of an a- 
particle and two positive electrons from four protons. 
The two positive electrons will react with two ordinary 
electrons to disappear into y-ray, the so-called annihila- 
tion of electron pairs. Then the net energy released per 
cycle is that given by the difference in mass of four 
protons plus two electrons, and an a-particle. Since 
the binding energy between the electron and nuclei 
particles is negligible, this difference is the same as the 
difference between four hydrogen atoms (four protons 
plus four electrons) and the helium atom (a-particle 
plus two electrons). This difference is equivalent.to 27 
Mev. A small fraction 2 Mev. of this energy goes into 
neutrinos and is lost. The remainder goes into radia- 
tion, which is finally emitted from the sun. 

Thus about 25 Mev. are set free for each helium 
nucleus produced in the stars, or 6 Mev. for each proton 
destroyed. One gram of sun material has been calcu- 
lated to contain about 2 X 10?* protons. Therefore, if 
all the protons can be converted into helium, the avail- 
able energy supply is 1.2 X 1074 Mev. per g. At 
present, the sun radiates about 1.2 X 10° Mev. per g. 
persec. At this rate, the energy supply will last about 
30,000,000,000 years. The rate of energy release can 
also be approximated by calculations based upon the 
teactions (8) to (13) and laboratory measured con- 


stants. Thus Bethe’s “carbon cycle’ is indeed the 
energy process in the sun. In fact, the ‘‘carbon cycle’’ 
is found to be the energy process for all the “main- 
sequence stars’’ from ‘‘red dwarfs’ to ‘‘blue giants.” 

On the physical atomic weight scale, four protons and 
two electrons total 4.0215 units. This means that 
4.0215 g. of the reacting material in stars will carry 
out 6.064 X 10?* processes, each of which generates 25 
Mev. energy. The energy generated by 1 Ib. of the 
reacting substance is then 


[1/(0.0040215 X 2.205)] X 25 X 6.064 X 107° X 


4.45 X 10-* & 3,413 = 2.59 X 10" B.t.u. (14) 


By comparing this value with that given by Eq. (3), it 
is seen that the helium manufacturing process is able to 
transform 0.67 per cent of the original reacting mass 
into energy. While this percentage may seem to be 
small, because of an enormous conversion factor from 
mass to energy, the ‘“‘heat value” per pound of reacting 
material is still tremendous. One of the most energetic 
molecular reactions is the combustion of hydrogen with 
oxygen. The heat value of the stoichiometric mixture 
is 6,850 B.t.u. per Ib. The heat value of the nuclear 
reaction given in Eq. (14) is thus 3.78 X 10’ times 
larger than that of the molecular reaction. This con- 
firms the previous surmise that the ratio of energy re- 
leases for the nuclear reaction and the molecular reac- 
tion must be of the order of a million. 


CHAIN REACTION 





NUCLEAR FISSION 


The second method of energy release by breaking 
down the extremely heavy elements depends upon the 
instability of the heavy nuclei under external excita- 
tion. The liquid drop model of the nucleus is also use- 
ful in the stability considerations for this nuclear fission 
reaction. In fact, by assuming that the drop is com- 
posed of an incompressible fluid of radius R and of 
volume (4r/3)R* = (42/3) ro°A, uniformly electrified 
to a charge Ze and possessing a surface tension O, the 
critical excitation energy, E, for separating into two 
equal nuclei is shown by Bohr and Wheeler® to be ex- 
pressed by 


E, = 4nr°0A™{(Z*/A)/(Z*/A)iimiting) (15) 


ro is the radius of a single nuclear particle and 
(Z?/A)timiting = 47.8, limiting value for Z*/A giving in- 
stability for infinitesimal deformation of the “drop.” 
E, is actually the height of the barrier which leads to 
instability. In other words, if a certain excitation has 
an energy of less than E,, the nucleus is stable against 
splitting into two fragments. If a certain excitation 
has an energy in excess of E,, then the nucleus can be 
deformed to such an extent that it can be further de- 
formed into two fragment nuclei without additional 
energy expenditure. The situation is similar to that of 
rolling a ball against a hump of ground. In order to 
pass the hump, one must apply sufficient energy for the 
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TABLE 3 

Excitation of Heavy Nuclei by Neutron Capture 
Initial Produce 
Nucleus Nucleus Z?/A x Ey E, Ey — E, 
920234 2 254 36.00 0.753 5.0My 5.4My —0.4My 
92235 9226 35.86 0.750 5.2 6.4 —1.2 
9228 92 289 35.40 0.741 5.9 5.2 +0.7 - 
aPu2® Pu §=36.80 0.770 4.0 6.2 —2.2 
gTh?3? goTh**3 34.76 0.727 6.9 5.2 +1.7 
91 Pa? Pa? 35.70 0.747 5.5 5.4 +0.1 





ball to reach the hump, then further motion of the ball 
is automatic. 

Consider now what would happen if a neutron were 
added to such a near unstable nucleus. The value of E, 
is calculated by Bohr and Wheeler and the result in- 
cluded in Table 3. After the capture of a “zero” 
kinetic energy neutron (thermal neutron), the equi- 
librium of the collection of neutrons and protons is dis- 
turbed. In general, a readjustment is possible with the 
subsequent release of energy to excite the product 
nucleus. The energy of excitation E, of the compound 
nucleus is then the binding energy between the initial 
nucleus and the neutron. This is also calculated by 
Bohr and Wheeler and the result listed in Table 3. 
Therefore, the threshold energy of neutrons required for 
producing the splitting of fission of a given compound 
nucleus is obtained by subtracting the excitation energy 
of zero kinetic energy neutron from the critical energy 
for fission. The result for this energy of the bombard- 
ing neutrons is listed in Table 3 as E, — E, for different 
bombarded nuclei. The negative values for U?*, U?*, 
and Pu?® mean that these nuclei will undergo fission 
after capture of thermal neutrons. These theoretical 
results are in full agreement with experiments. 

However, there are other competitive reactions that 
could happen to an excited nucleus. These reactions 
are: (a) emission of y-ray, (b) emission of a-particle, 
and (c) re-emission of the neutron. To estimate the 
relative probability of occurrence, one must turn to the 
statistical mechanics of the excited nucleus. It is found 
that, for energies of excited nucleus below, or not far 
above, the critical value, radiation of y-ray will pre- 
dominate. With increasing energy the probability of 
fission will increase rapidly. The probability of fission 
will not increase indefinitely, however, because at still 
higher energies the re-emission of neutrons becomes the 
predominant process. 

Now there still remains the problem of making the 
capture possible. Of course, the nuclei will be bom- 
barded with neutrons, but the probability of actual 
absorption depends upon the kinetic energy of the bom- 
barding neutron and the character of the bombarded 
nucleus. This probability is generally expressed as a 
“cross section”’ o defined by 


N= (16) 


nvo 


N is the number of occurrences of the kind in question 
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(fission) per sq.cm. of bombarded surface per sec.; 1 ig 
the number of bombarding particles per sq.cm. per see, 
(normal incidence around); and yv is the number of 
atoms of the kind involved per sq.cm. of surface. 
situation is as if each atom presented a target area g, 
If the target is struck, the effect in question occurs; if 
not, nothing happens. ~ It is apparent that o gives a 
sort of averaged probability of occurrence of an 
effect for a random spatial distribution of impinging 
particles. 

It is thus evident for the successful fission process 
that two conditions must be satisfied: (a) a large cross 
section of capture of neutrons of proper energy; (b) the 
energy of the neutron must be such that the excitation 
energy of the compound nuclei, after capture of neu- 
trons, is above the critical barrier energy. It is found 
that »U** and yPu?® can be split by either thermal 
neutrons or fast neutrons. Then, as explained before, 
a very large energy release is expected. In fact, the 
measured energy release per fission for 9,.U** is approxi- 
mately 177 Mev.’ Since 1 gram mole or 235 g. of 
2U?* contains 6.064 X 107’ atoms, the energy generated 
by 1 Ib. of 9.U** with complete fission is 


[1/(0.235 X 2.205)] X 177 X 6.064 X 10% x 


4.45 X 10-* X 3,413 = 3.14 X 10" B.t.u. (17) 


By comparing with the value given by Eq. (3), it is 
seen that the fission of U-235 is able to convert 0.081 per 
cent of the original mass into energy. Thus, the reac- 
tion is not so efficient as that of helium manufacture 
from hydrogen, treated in the previous section. This 
is, of course, expected from the diagram for the packing 
fraction (Fig. 2), which shows a much greater change in 
binding energy per nuclear particle from hydrogen to 
helium than from uranium to medium heavy elements. 
Nevertheless, the ‘‘heat value’ of the U-235 is still 
approximately 1.5 X 10® times the heat value of gaso- 
line, which is the most powerful fuel used today. 
However, to make the fission process a really going 
process, a continuous supply of neutrons of proper 
energy is necessary. This is actually accomplished by 
slowing down the fast neutrons released by the fission 
fragments themselves. For »U**, the compound 
nucleus .U?*6 first breaks down into sBa and Kr of 
unequal mass, since this splitting involves smaller 
potential barriers. The product nuclei would be highly 
excited because of the energy released and would con- 
tain more neutrons than stable equilibrium requires, 
because the ratio of neutron to proton is much higher in 
uranium than for atoms of medium atomic weight. 
These extra neutrons will be ‘‘vaporated’’ or emitted 
with high energy and a small time delay of about 1 see. 
For one fission, more than one neutron is thus genet- 
ated. If these-neutrons can be slowed down to proper 
energy for fission of new ».U?* atoms without absorbing 
so many of them as to reduce the desired neutron to less 
than one per fission, the reaction becomes a chain reat- 


tion. The chain carrier is the neutron of proper energy. 
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Then the reaction will tend to build up to explosive 
violence. 

The agents for slowing down the fast netitrons are 
called moderators. They must not absorb the neutrons 
above the desired energy; otherwise the useful yield 
will be extremely small. They must also be light ele- 
ments so that they can effectively share the kinetic 
energy of fast neutrons by collision and thus slow down 
the neutrons with but few collisions. The thickness and 
the bulk of the moderator can thus be reduced. Hydro- 
gen is good because it is the lightest element, but it is un- 
desirable for slow neutron fission because of its large 
capture cross section of slow neutrons. Deuterium is 
best from this aspect, and heavy water is a good 
moderator. 

Besides the absorption of neutrons by the moderator, 
there is another loss due to escape of the neutrons out 
of the reacting mass. This is proportional to the sur- 
face of the mass. Therefore, by increasing the size of 
the reacting mass, this loss can be cut down relative to 
the production, which is proportional to volume. 
Hence, if a small quantity of fission material and 
moderator cannot carry on the chain reaction, a large 
quantity of the same material might be able to do so. 
For any given geometrical arrangement and effective 
composition, there is thus a critical size. Generally, 
single fission can be produced automatically or by cos- 
mic ray; therefore, once the critical size is reached, the 
mass will react by itself or explode. The atomic bomb 
is based essentially on the principle of bringing parts of 
subcritical size together quickly to build up the reac- 
tion. The rate of reaction and particularly the rate of 
building up the reaction are controlled by the diffusion 
of the neutrons through the moderator to the fission 
nuclei. For fast build-up in explosion, fast neutrons 
are preferred. 


ENGINEERING APPROACH TO THE NUCLEAR REACTION 


The previous discussion clearly shows that the funda- 
mental concepts involved in a nuclear reaction—such as 
collision, capture, excitation, and energy barrier—are 
all familiar to students of molecular reactions. The 
study of the molecular reactions or chemical reaction 
from kinetic point of view is, however, only a recent 
development. Its consideration by engineers is even 
more so. For instance, the power-plant engineers 
learned to design large boiler installations burning tons 
of coal per hour without any concept of molecular 
physics. The design of internal combustion engines 
was no exception, until the difficulty caused by detona- 
tion became serious. Then the investigation of the 
slow oxidation of the hydrocarbons by chain reactions is 
helpful. However, the study of nuclear reactions took a 
different road. Here, the effort was concentrated 
immediately on the detailed dynamic processes and on 
the development of a fundamental theory for the nu- 
clear dynamics to explain these processes. This, of 


course, is the natural result of the emphasis by physi- 
cists on the understanding of the fundamentals. 


For an engineer, his assignment is not so much to 
understand a particular phenomenon as to learn to 
utilize it. Electrical engineers knew how to design 
long-distance transmission lines before the development 
of satisfactory physical theory for the conduction of 
electricity. A similar situation will certainly appear in 
the engineering utilization of atomic energy. A satis- 
factory fundamental theory for nuclear structure will 
no doubt be a great help in “atomic engineering,” but 
to wait for such a theory is certainly unwise and un- 
necessary. Then shall the ad hoc experiments without 
coordination be encouraged? Such an approach would 
not only be uneconomical but extremely dangerous be- 
cause of possible uncontrollable explosion with greater 
violence. A better approach would be a semi-empirical 
one. In other words, first the performance of each 
elementary process should be determined experimentally. 
Then the kinetic theory of reaction devcloped for the 
chemical kinetics can be used to predict the perform- 
ance of the overall reaction involving these elementary 
processes. 


For instance, the basic knowledge required for pre- 
dicting the fission reaction is the fission cross sec- 
tion; the elastic scattering cross section; the in- 
elastic scattering cross section of the atoms U-234, U- 
235, U-238, and Pu-239 as a function of the energy of 
the neutrons; the number and energy of neutrons 
emitted by fission and its dependence on the energy of 
incident neutron; etc. With these empirical data, the 
fission reaction in an atomic bomb can be calculated. 
To those familiar with the theory of chemical kinetics, 
the extreme complexity of the treatment due to the 
great variety of possible molecules is well known. The 
application of kinetic theory to nuclear reactions, how- 
ever, may be much simpler, since the total number of 
possible nuclei is only several hundred, while the num- 
ber of known molecular compounds may be hundreds 
of thousands. The reacting components that need to 
be considered in a nuclear reaction will certainly be 
fewer than for the case of molecular reactions. 


According to H. D. Smyth’s report‘ on the develop- 
ment of the atomic bomb, this semi-empirical approach 
is actually followed. Similar methods must be used for 
the development of the general utilization of atomic 
energy to determine the optimum process or processes 
and then finally test the validity of calculation by ex- 
periments. This approach would certainly eliminate, 
as far as humanly possible, the danger and the waste 
that through improper handling might be associated 
with the release of energy a million times as great as the 
conventional combustion process of molecular reac- 
tion. It is seen then that the concept of nuclear 
power engineering is necessarily somewhat different 
from its counterpart in pure science, the nuclear 
physics. 
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k Letter to 


Dear Sir: 

The following comments on the Wright Brothers Lecture 
(‘British Aircraft Gas Turbines,’’ by Dr. H. Roxbee Cox, 
JOURNAL OF THE AERONAUTICAL SCIENCES, Vol. 13, No. 2, p. 
53, February, 1946) are submitted to you as my personal points 
of view and do not necessarily correspond to those of the Bureau 
of Aeronautics. 

As a history of gas-turbine developments in England the paper 
is an extraordinarily interesting document and contains some ma- 
terial with which I was not familiar heretofore. The early history 
of the struggle for favorable consideration of such devices is yot 
unlike similar history in this country. For example, the con- 
clusions of a special panel formed in 1939 to investigate the ap- 
plication of gas turbines to aircraft is worthy of note: ‘‘In its 
present state, and even considering the improvements possible 
when adopting the higher temperatures proposed for the im- 
mediate future, the gas turbine could hardly be considered as 
feasible of application to airplanes mainly because of the diffi- 
culty in complying with the stringent weight requirements of 
aeronautics. 

*.,. The minimum weight for gas turbines even when taking 
advantage of higher temperatures appears to be approximately 
13 to 15 pounds per horsepower.””’ The Navy was fortunate in 
having officers like Captain Kauffman of sufficient vision to pro- 
ceed, despite this lugubrious outlook, to the planning of, and 
final contracting for, a gas turbine, the first such project in this 
country for aviation purposes. These studies led to a joint 
Navy-A.A.F. project under Navy cognizance started early in 
1941. Later in the same year the first United States turbojet 
project was undertaken by the Navy at Westinghouse. I do not 
doubt that this information, or part of it, was known to Dr. 
Cox, but it does illustrate the difficulties involved in getting new 
ideas into development. 

It may have occurred to some persons to wonder why the 
United States did not follow the British example in setting up a 
collaboration committee. Surely the excellent results obtained 
by this method in England is a strong argument in favor of such 
a procedure. However, it was, and still is, my opinion that the 
idea would not have worked so well in this country. The dif- 
ference in the situation here and in England as to urgency, avail- 
able facilities, and methods of handling design and development 
problems are obvious. In this country we depended on liaison 
officers like Lieutenant Commander Bollay of the Navy and 
Colonel Keirn of the Army Air Forces to effect proper distrikution 
of essential design information and development experiences. 
If this system delayed our progress somewhat, it nevertheless 
resulted in several distinct designs, each having its own merits. 
It has long been my contention that competition is an essential 


the Editor 


to proper development in any field where competition can be 
achieved. There is little doubt that collaboration between 
England and the United States in this field was complete and of 
great value to both countries. 

In considering Dr. Cox’s detail comments on troubles experi- 
enced in gas-turbine development, our difficulties seem to corre. 
spond closely with the British with minor exceptions. I do not 
feel that we have completely licked the turbine disc problems 
in this country, although we have had little recent trouble. 
However, I shall not consider the problem completely solved 
until such time as we have suitable turbine discs that can be 
manufactured on a high production basis. The annular com- 
bustion chamber appears to be successful in Westinghouse engines, 
somewhat contrary to British experiences. I do not feel that the 
convenience of the multiple chamber from a development and 
handling standpoint should be a primary design consideration if 
the engine suffers in any respect thereby as an engine. Naturally 
the axial flow engine can less afford the increased size arising 
from the use of multiple chambers, since small diameter is one of 
its major advantages. 

The facilities described as being available to Power Jets are 
comparable to those available in the component test laboratory 
at Westinghouse and at other places in this country. However, 
it is my opinion that more extensive facilities for the test and 
development of full-scale engines is essential to eventual success 
in this field. Although the Navy has used flying test beds for 
some time to determine at altitude such things as starting char- 
acteristics, stable combustion ranges, and minimum speeds, this 
type of test equipment is too severely limited to be considered as 
other than an adjunct to laboratories capable of full-scale alti- 
tude test work. 

The question of thoroughly integrated airplane and engine 
design is one that arises from time to time in everyone’s mind. 
Steps taken in this direction in past years have resulted in major 
improvements, such as new types of cowling, quickly removable 
power plants, etc., as well as revisions of engine shape and ar- 
rangement. However, a completely integrated engine and ait- 
plane has never been built in this country. Nor does it appeaf 
probable that the advent of the gas turbine will appreciably 
change the situation. There are gas-turbine designs that will 
permit some rearrangement of components to fit specific airplane 
designs, but, in general, gas turbines are no more suitable for such 
treatment than are reciprocating engines. In my opinion, the 
desire for standardized engines applicable to a wide variety of aif- 
planes will always prohibit separate engine designs for specific 
airplanes. It is possible that ram jets or rockets may readily lend 
themselves to such treatment, particularly in pilotless aircraft, 
but even here standardization and ease of production must be 
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carefully considered. The views of the airplane and engine 
people are normally at almost complete divergence on this ques- 
tion, and naturally so. The engine manufacturer must build in 
competition the best engine as an engine; the airplane manu- 
facturer realizes that all too frequently the success or failure of 
his airplane depends on the efficiency of his power-plant installa- 
tion. 

It is perfectly true, as the paper points out, that there is in- 
adequate evidence to indicate the superiority of either centrif- 
ugal or axial flow blowers at the moment. However, centrif- 
ugal blowers are limited as to pressure ratios if multiple staging 
js not used, and it would appear that the greater capacity of the 
axial flow blower may limit the use of centrifugal blowers to 
moderate powers in order to remain within reasonable engine 
diameters. On the other hand, the axial blower suffers from two 
plaguing difficulties for which complete cures are not yet avail- 
able. The first is its narrow operating range, which leads to 
instability and surge at altitude, and the second is its starting 
difficulties due, in part, to its narrow operating range and, in part, 
to the little assistance obtained from ram airflow as well as to its 
normally higher design pressure ratio. 

The discussion of the relative merits of turbojets, gas turbines, 
and ducted fans is interesting and, in general, corresponds to 
U.S. Navy ideas. Comparing the three systems alone, the 
turbojet appears to be most generally applicable to the highest 
speeds; the other two, to the greatest ranges. Such a generaliza- 
tion is not particularly pertinent, however. All things considered, 
the gas turbine appears to me to be most suitable for U.S. Navy 
use, offering the most advantages and the fewest disadvantages. 
Present indications are that propeller efficiencies may be in- 
creased markedly at high Mach numbers. If this proves to be 
true the turbojet airplane will offer little advantage in high speed 
and will be inferior in take-off, climb, range, endurance, and 
general flexibility of operation. 

As the paper points out, the consideration of how to select the 


best engine for a particular design is based on a wide variety of* 


factors, only one of which is the engine power curve with its 
corresponding fuel consumption curve. All things being equal the 
commercial ‘‘useful load ton-miles per gallon”’ is a better criterion 
in the selection of an engine for a given design than is pounds of 
fuel per pound of thrust. That number includes consideration of 
take-off power, part throttle performance, differences in engine 
weight and fuel weight, variations in fuel consumption and in 
drag, etc. It cannot be expected to cover such things, however, 
as desirable operating altitudes; penalties acceptable for higher 
operating speeds; or desired airplane characteristics such as 
endurance (as contrasted to range), maneuverability, flexibility of 
loading and operation, stability, etc. Nor can it cover extraneous 
considerations such as whether or not a flying fuel tank can make 
an acceptable military airplane. 

For commercial aircraft, economy would appear to be a prime 
consideration, but economy must include consideration of first 
cost, reliability of operation, maintenance costs, and premiums 
to be expected from higher operating speeds, in addition to the 
usual ‘‘ton-miles per gallon.’’ The military operator must de- 
fine economy primarily in terms of losses of aircraft and crews 
in combat, although he must also consider, from the U.S. Navy 
standpoint, day-by-day operational readiness. It is obvious that 
an airplane that must remain on deck from a lack of engine re- 
liability or durability or from lack of fuel is worse than useless and 
that an airplane is just as surely lost if it has inadequate endur- 
ance or if it cannot be safely operated as it is if it has inadequate 
performance. 

Of interest in this connection might be a summary I prepared 
several years ago near the beginning of the Navy’s gas-turbine 
Program. 

“From the standpoint of thermal efficiency alone, true jet 
Propulsion is at one end of a scale which has the propellered engine 
at the other, most efficient, end. The difference is very marked, 


being on the border of three toone. Assuming equal weights and 
drag, a conventional airplane which has a range of 1,000 miles will 
have a range of only 300 miles as a true jet propulsion airplane. 
The difference in endurance is even more marked, being on the 
border of the ratio of 1 to 4 hours at equivalent percentage powers. 
Basically the true jet propulsion unit is a full throttle engine and, 
hence, suffers in endurance. A possible way of reducing this 
discrepancy would be to have multiple units so that the air- 
plane could be operated at full throttle on whatever percentage 
of units would be required to obtain the desired power. At best, 
however, there must always be consideraple loss in efficiency, 
assuming always, of course, that the propeller does not lose too 
much efficiency at higher speeds (which does not seem unrea- 
sonable at the speeds we are now discussing). 

“In weight again the true jet propulsion is at one end ef the 
scale, possibly as low as half pound per equivalent brake horse- 
power, the propellered engine being 1.1 to 1.6 Ibs. per b.hp. 

“From the standpoint of airplane design the true jet propulsion 
machine offers several advantages. Its small size leads to smaller 
fuselages and possibly better streamlining than the conventional 
propellered airplane. Furthermore, since propeller clearance is 
not involved, the landing gear need only be big enough to keep 
the fuselage off the deck, and therefore the landing gear weight is 
considerably reduced. In this connection it should be pointed 
out that a propeller need not be of the conventional design now 
familiar but could readily be a shrouded (pro) impeller, an axial 
fan, or an exhaust inducer, and therefore part of this advantage 
can be eliminated. 

“From a performance standpoint the jet propulsion unit, for the 
same weight airplane and assuming the airplane is still below the 
critical compressibility range, will give substantially constant 
speed up to its critical altitude. However, for the same weight, a 
turbosupercharged engine will produce more (equivalent) brake 
horsepower at altitudes about 35,000 ft. and, assuming reason- 
able propeller efficiency, more speed. Hence, in speed ranges 
up to perhaps 450 m.p.h., the jet propulsion unit shows tremend- 
ous advantages of speed and climb at low and intermediate 
altitudes but approaches equality at high altitudes. For speeds 
above 500 m.p.h. it is probable that propeller efficiency will be so 
low as to nullify any advantage at high or low altitudes, except in 
The take-off with a true jet propulsion unit will be 
The 


cruising. 
poor in any case and must be augmented in some form. 
range question has already been considered in paragraph 2. 

“From the standpoint of a usable airplane the true jet propulsion 
aircraft must necessarily await the solution of the following prob- 
lems: controllability; balance; handling characteristics at low 
and high speeds; take-off; carrier approach; and gas disposal, 
particularly with regard to damage to other airplanes in carrier 
take-offs and to flame visibility. In this connection a compromise 
solution in which part of the power is absorbed in a propeller will 
reduce the energy in the jet and may reduce the above problems 
considerably. 

‘After due consideration of all the factors involved, this Section 
recommends as an eventual desirable solution a compromise 
power plant in which the division of power between the propeller 
and the jet approximates the theoretical best division—i.e., 75 
per cent to the propeller and 25 per cent to the jet. This setup 
might include a means of adjusting this percentage in the air so 
that the pilot may elect to get high speed low range at sea level, or 
long range, or any compromise between the two. The propeller, 
as previously pointed out, may take the form of a fan or impeller. 
Pending this eventual solution, it is recommended that true jet 
units be installed as accessories to conventional internal com- 
bustion engines in order to obtain better low-altitude performance 
without sacrificing too much in range or endurance. It is the 
present hope of this Section to have one or more such units of 
different types ready for the flight test in early 1943. A true jet 
propulsion installation is not recommended at this time.” Quali- 
tatively, this analysis still appears valid. 
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The lecture presents some ideas relative to possible airplane 
changes because of the advent of the gas turbine. Unquestion- 
ably, some changes will be called for by the advent of long, small 
diameter power plants and by higher airplane speeds. That those 
changes will result in tail first or no tail airplanes is open to con- 
siderable question. We are all more or less familiar with the 
defects of this type of aircraft, particularly from a military stand- 
point. Furthermore, available test data indicate the slipstream 
effect at high speeds to be small either as a pusher or as a tractor. 
The maximum slipstream drag is probably no more than 2 per 
cent of the total drag in the extremely high-speed condition. The 
tractor arrangement is preferable at low speeds, as in the ap- 
proach condition, because of the much lower stalling speeds avail- 
able with tractor propellers with power on. 

One effect not covered as such in the paper is the radical change 
in center of gravity location, which may be expected in airplanes 
whose major useful load is in fuel, between the full load and no 
load condition. Stability and controllability at low speeds is 
likely to be an extremely serious problem in turbojet airplanes 
in which the designer tries to obtain ranges and endurances com- 
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parable to those obtainable in propellered aircraft. The Me 163 
is an exaggerated example of this effect, since it is a rocket air. 
plane, but that airplane is almost uncontrollable at low speeds 
with no fuel. 

It is assumed that the paper omits a discussion of the so-calle 
compound engine, since that arrangement is not considered g 
true gas turbine. However, combinations of reciprocating 
engines and gas turbines have much to recommend them. The 
use of an Otto or Diesel cycle engine as compressor and com. 
bustion chamber offers higher overall efficiencies than appear 
obtainable in other ways. The free piston compressor, if applied 
to a gas turbine, would be, in effect, a combination of recipro. 
cating engine and turbine. Admittedly, the result of such g 
martiage is a complicated and bulky engine, but, if past experi. 
ence is a suitable criterion, simplicity will always give way to 
performance. 

S. B. SPANGLER 
Captain, U.S.N, 
Bureau of Aeronautics 
Navy Department 








Sciences offers the facilities of: 


years of age. 


ties will be sent on request. 


Library Facilities of the 
Institute of the Aeronautical Sciences 


To serve Institute members and others in the aeronautical industry, the Institute of the Aeronautical 


The W. A. M. Burden Library 
2 East 64th Street 
New York 21, N.Y. 


The facilities of this library are available for reference study at the Institute. 


The Paul Kollsman Lending Library 
2 East 64th Street * 
New York 21, N.Y. 


This library loans books without charge to members and others in the United States over eighteen 
Full information will be sent on request. 


The Pacific Aeronautical Library 
6715 Hollywood Boulevard 
Hollywood 28, Calif. 


This is a service and reference library for West Coast organizations. 


Fuli information as to its facili- 

















> Me 163 
cket air. 
W Speeds 


So-called 
idered a 
rocating 
m. The 
nd com- 
| appear 
f applied 
recipro- 
' such a 
t experi- 
- way to 


NGLER 

, U.S.N. 
onautics 
artment 








Flapping Characteristics of Rigid Rotor 
Blades 


RALPH W. ALLEN 
Consulting Engineer 


SUMMARY 


The influence of blade stiffness on rotor blade flapping is in- 
vestigated and equations are developed which describe the blade 
motion. The blade lag angle is also obtained and an example 
calculation is given. 

It is found that blade stiffness has a marked influence on blade 
lag and flapping amplitude and that there exists a possible reso- 
nant condition for the semirigid type of rotor blades. Increasing 
the blade stiffness decreases the magnitude of flapping and re- 
duces bladelag. The stiffness effects are especially noticed in the 
transition from a freely hinged blade to a semirigid blade. 


INTRODUCTION 


ep THE AERODYNAMICS or structural quantities 
of the present helicopter rotor blade systems can 
be properly determined, it is necessary to know the 
characteristic motions of the blades as they rotate. 
An important motion of the blade system is flapping, 
which occurs mostly when the helicopter is in motion. 
Blade flapping is generally brought on by the differences 
of lift distribution between the advancing blade and the 
retreating blade. Each blade tends to seek its equilib- 
rium position in the vertical direction and in doing so 
causes the blade to flap. The flapping changes the 
angle of attack of the blade thereby changing the lift 
distribution. These changes occur in a cyclical man- 
ner around the azimuth. 

The majority of the present helicopters have their 
blades hinged somewhere near the center of rotation, 
thus permitting the blade to flap freely and seek its 
own equilibrium position in the vertical direction. 
The investigation of the flapping motion of the hinged 
rotor blade has been developed not only in theory but 
has been experienced in many flight tests.'. Recently, 
there has been a definite swing toward rigid-type rotor 
blade construction, and the flapping motion of this 
type of blade has not been greatly investigated. In the 
tigid-type blades, such as is used in many coaxial-type 
rotor blade systems, the flapping motion of the blade is 
restrained by the elasticity of the blade material and 
will, therefore, have a different flapping pattern about 
the azimuth. The positions of maximum and minimum 
flapping will be different in the rigid blade than in the 
hinged blade. This paper analyzes the flapping mo- 
tions of a rigid-type rotor blade and determines the 
angle of lag, which is defined as being the azimuth 
angle between the maximum thrust position and maxi- 
mum flapping position of the rotor blade. 
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DEVELOPMENT OF THE FUNDAMENTAL EQUATIONS 


In order to simplify the equations, it is assumed that 
the blade in this analysis is perfectly rigid and develops 
all of its bending at the root section of the blade. The 
elasticity of the material at the root section is repre- 
sented by a torsional spring of stiffness factor K. A 
rotor blade system that mechanically duplicates an 
actual rotor blade motion is schematically shown in 
Fig. 1. 
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Fic. Schematic setup of blade flapping mechanics and forces 

acting on blade. 
NOMENCLATURE 

8 = flapping angle 

y = azimuth angle measured from downwind direction 

ao, %, 5b; = flapping constants 

I = moment of inertia of blade in flapping 

t = time 

K = blade stiffness constant 

w) = natural frequency of blade in bending 

w = rotational velocity of rotor blade 

M, = thrust bending moment 

M, = blade weight bending moment 

R= blade radius 

c = blade chord 

a = lift slope factor 

X = down-wash factor 

6 = blade pitch 

6, = blade twist factor 

uw = tip speed ratio 

B- = blade tip loss factor 

¥, = blade lag angle 


The moment due to the centrifugal force is 
S Bw'rtdm = Bw*l 
Assuming cos 8 = l, 
sin B = B 
The moment produced by the inertia force is 
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JS (D°8/DE)r'dm = (D*B/D#2)I 


The restoring moment of the spring is K8. 

These moments sum up to the left side of the dif- 
ferential Eq. (4). 

The flapping angle 8 is expressed in the form of a trig- 
onometric series 


B = a — a, cosy — sin y — azcos 2y — b, sin 2y (1) 
Neglecting second harmonic terms 6 becomes: 
= a — aqcosy — h siny (2) 
Eq. (2) is sometimes written in the form 
B = a — a, cos (¥ — Ys) (3) 


The differential equation expressing the relationship 
of the aerodynamic moments to the mechanical mo- 
ments is 


I(d*B6/dt?) + KB + Iw?B = M, — M, (4) 


The aerodynamic moment , is as follows: 


M: = 5B + 70 (BY + 2B2) + 
V/apcaw*R* 3 ° M 


5 
= (a + 6 u*B*) + +( HOB * + = 5 HOB! — 
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for its magnitude is far less than the aerodynamic 
moment and therefore negligible error should be in- 
troduced. 


SOLUTION OF THE DIFFERENTIAL EQUATION 


The solution of Eq. (4) is found by substituting the 
value of B into Eq. (4) and equating coefficients of like 
trigonometric functions; thus 


d?B/dt? = w*(a, cos y + b, sin y) (6) 
Substituting Eqs. (6) and (3) into Eq. (4) yields 
w?(a, cos Y + b, sin Y) + 27(a — a, cosy — db) sin y) = 
M,/1 (7) 
2? and wo? = K/I. 


where w? + w? = 


Equating coefficients of like trigonometric functions, 
the following values for the constants are obtained: 


e a 


me ee 


E AB? + - vs (B* + p2B?) + 


6, (a +2 u*5) (8) 


orl 


y = ‘/2pcaw*R* 


‘ ] 1 1 
; 1) a (wT = — 5 uBio + GHB + 5 OB (0) 
4 uBs +- 9 MAB? + - ae, sin y + : : 
By ime 1 1 
bi(w? — 2?) — = 5 whB* + = uO, Bt — — a, B4 
(-3 = MaoB* + = 7 DBs +5 u’b,B ) cos io y y¥ 3 a a tp + 
1 1 
2 — 42 2 
The gravity moment M, is neglected in this analysis, 2 © 3" eS 
From Eqs. (8), (9), and (10) a; and 2; are obtained. 
Ne | eee a ge Ce ee ee, ns aan 
a, = — 3 HB ‘Ao > — 9?) + & ¥%B + 3 #48 + 5 wAB 7 (w? — 07)? + 16 P (11) 
For a blade with no twist and neglecting u* terms 
- Sips re (wo? — 92) me 1 ps4 typst /[F (co? — 02)? + 1 Bs (12) 
\ 3 pat Be | gs f/ Ly 16 . 
Neglecting u* terms 
re eee 1 1 a | w mines 
bh = = 2")\ 3 HOB + 5 uXB® + 5 uiB*) + 75 uB 0 | -- (w? — 0%)? + 75 B® (13) 


with yu? terms included, b; becomes 


I 2 1 ig 
b, = E (w? — 0y(; pO,B* + = ; * XB? +; 5 HOB ‘) + BD usr |/[ (w? — 07)? +7 ee 6 a pBs + — 39 5 mB] 


y 


For a blade with no twist 


of as cign (rep a2 bp P tg: 
bh =|B oe — 2?) 3 % +54 + yo Beto aT 


(14) 


9°)? + — 6 | (15) 


Substituting the equations for }; and a; and into Eq. (2) completes the solution of the flapping motion of the 


blade. 
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The minimum flapping position is found by differentiating Eq. (16) and setting to zero; thus 


_ (—I/*/2 pcaR*) (w/w?) (?/s3 u00.B? + '/2 udB*) + */12 uB"do 





tan y;, = 


(1/3 uB*Ido/1/2 pcaR*)(w?/w*) + 1/6 u0o.B" + 1/3 urB® 


This equation is shown in curve form in Fig. 3, for untwisted blade. 


EXAMPLE CALCULATION 


An example calculation is made to show the effects 
in curve form of the parameter w/w. The necessary 
data for the blade are as follows: 


I = 334 slugs ft.? 

R = 22.5 ft. 

C = 1.25 ft. 

6 = 0.06 rad. 

zp = 0.5 

p = 0.0021 slug per ft.* 
a = 5.58 


The equations that develop from substituting the 
data into Eq. (16) are 


“ — 0, 8 = 0.12190 — 0.10000 cos y — 
0.08193 sin y 
“ — 1, 8B = 0.06090 — 0.15020 cos y — 
: 0.0079640 sin y 
~ = 2, B = 0.02430 — 0.17550 cos y — 0.02357 sin y 


” — 10, 8 = 0.001219 — 0.001675 cos y — 
w 
0.0002 4996 sin y 


These equations are shown in curve form in Fig. 2. 

The lag angle is found by substituting the necessary 
data into Eq. (17). The results are shown in curve 
form in Fig. 3. 


PROCEDURE IN DESIGN ANALYSIS 


In the preliminary design of rotary aircraft it is neces- 
sary to know the flapping characteristics of the rotor 
blades. Aside from structural reasons, interference be- 
tween blades, as in the coaxial-type rotor system, and 
aerodynamic efficiency are some of the important infor- 
mation needed in the proper design of present-day rotary 
aircraft. 

The procedure necessary for the analysis of blade 
flapping consists of, first, setting up the design criteria— 
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measured from the down-wind direction. 


that is, such fixed quantities as weight of aircraft, speed 
of rotation, blade radius, number of blades, airfoil 
section, maximum ceiling, and maximum speed. From 
this information the factors %, A, and uw can be found. 
Knowing these factors, the flapping coefficients do, a, 
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and }; can be obtained. By substituting the coefficients 
into the flapping equation for an assumed J and w/Q, 


B =a — aqjcosy — hsiny 


the flapping amplitude can be found and the maximum 
and minimum position along the azimuth obtained. If 
the results are not desirable in accordance with the par- 
ticular design, the process is repeated using different 
values of J and w/Q until the necessary and proper 
blade flapping characteristic is determined. The mass 
moment of inertia and blade frequency will then be 
found and, from this, the blade weight. The blade is 
then designed to meet these values. 


CONCLUSIONS 


Blade flapping is found to be a function of the fre- 
quency ratio w/w, which is a measure of the blade 
stiffness; thus a higher ratio indicates a stiffer blade. 
It has been found that the greater this ratio, the less 
becomes the lag between maximum thrust and flapping 
position. 


Apparently, the transition from a freely hinged blade 
to a semirigid blade produces marked changes in the 
lag angle, as shown in Fig. 3; however, as the blade 
stiffness is greatly increased the lag angle change is 
small. 

The amplitude of flapping is reduced by increasing 
the blade stiffness, other factors remaining the same, 
From the curves in Fig. 2, the flapping amplitude is 
negligible under a frequency ratio of w/w = 10. It is 
also noted that the flapping angle reaches fairly high 
negative values, indicating that there exists a critical 
flapping condition that occurs near a frequency ratio 
of 2 for this particular example. It is, therefore, im- 
portant to investigate, in the preliminary blade design 
stages, the flapping characteristics of the rotor blades 
so that a critical resonant flapping condition will not 
persist. This is especially needed for semirigid blades, 
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Buckling of a Flanged Plate 


HENRY L. LANGHAAR* 
Consolidated Vultee Aircraft Corporation 


ABSTRACT 


This investigation treats the buckling of an axially compressed 
isotropic elastic rectangular plate that is hinged on the ends, sup- 
ported but subjected to elastic rotational restraint on one longi- 
tudinal edge, and flanged but otherwise free on the remaining 
edge. A numerical procedure with a rigorous mathematical 
basis is proposed for solving elastic stability problems of this type. 


INTRODUCTION 


3 then GENERAL PROBLEM under consideration presents 
no inherent mathematical difficulties, but it leads 
to a finite equation that is too complicated to be directly 
useful for engineering stress analysis. The following 
work relieves the labor that is encountered in individual 
solutions of the problem by the use of numerical tables 
of the transcendental functions that constitute the 
snarl in the algebraic result. The tabulations may be 
regarded as an extension of data charted by Timo- 
shenko! for solving special problems of the present type 
—namely, problems wherein the elastic rotational re- 
straint of the supported edge is infinite (clamped edge). 
It should be observed that the condition of simple sup- 
port at the ends of the plate is not an essential restric- 
tion if the plate is long, since the buckling stress of a 
long plate is nearly independent of the end conditions. 

The results that are presented are a logical basis for 
treating certain aircraft structural problems, such as 
the buckling of the outstanding leg of a lipped or a 
bulbed angle. With a view to these applications, it may 
be mentioned that a practice that has been adopted by 
many aircraft stress analysts is to present plate-buck- 
ling data throughout the elastic and the semiplastic 
ranges by a single curve of o,, versus b/(t\/K), the 
abscissa b/(t~/K) being derived from the universal 
plate-buckling formula 


Ger = KE(t/b)? (1) 


just as the column-curve abscissa L/(p+/C) or L’/p is 
derived from the Euler column formula. It is not 
within the scope of this study to pursue theoretical or 
experimental investigations of the accuracy of this pro- 
cedure nor to inquire whether variation of the coefficient 
K is required when the proportional elastic limit is ex- 
ceeded. Nevertheless, it is pertinent to note that there 
is, at present, no more reliable method for predicting 
buckling stresses of highly loaded plates than to ascribe 
to the coefficient K a significance beyond the confines 
of elasticity theory. 


Received July 17, 1945. 
* Structural Research Engineer. 


GENERAL EQUATION FOR THE BUCKLING COEFFICIENT 


A plate (Fig. 1)—hinged on the edges x = 0 and x = 
a, supported and subjected to elastic rotational re- 
straint on the edge y = 0, and flanged on the edge y = 
b —is buckled by an axial compression stress o,,. The 
thickness of the plate is ¢, and its flexural rigidity is D = 
Et®/12(1 — v?). The cross-section area of the flange 
is w, and the moment of inertia of the flange cross section 
with respect to its centroidal axis parallel to the plane 
of the plate is J. The depth of the flange is assumed to 
be small compared to the width b of the plate. 

















Fic. 1. Diagram of flanged plate showing nomenclature. 


The deflection of the plate due to buckling is w. The 
spring constant for the rotational edge restraint is C; 
i.e., the moment M, per unit length on the edge y = 0 
is M, = —C(0w/Ody). The number of half waves in the 
buckled plate is m. 

The following dimensionless nomenclature is em- 
ployed: 


g 


mx/a, = y/b 

(mrb/a)? 

12(1 — v?) [I/(6t*)] 

w/ (bt) 

= Cb/D (2) 

O-7tb?/D 

functions of \ and y defined by Eq. (15) 
and tabulated in Tables 1, 2, and 3 

F(A) = a function defined by Eq. (4) 


Y 
A 
B 
K 
rN 
f, gh 


By virtue of its definition, the variable ) is related to 
the coefficient K of the buckling formula (Eq. (1)) by 
the equation 


K = 4/12(1 — ?) (3) 
The theory of elasticity (see Appendix) furnishes the 


following equation for the coefficient \: 


Ay - r+ rt8t/7_ may Je? &® 


The quantities f, g, and / are functions of the two vari- 
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TABLE 1 ie cil 
The Function F(A, y). (v = 0.30) 
my -. ai — 
4 3 5 10 15 20 25 30 35 40 45 50 55 60 ; 
0.0 -—3.00 -—3.00 -—3.00 -—3.00 -—3.00 -—3.00 -—3.00 —3.00 —3.00 —3.00 —3.00 —3.00 —3.00 0 
0.2 —3.15 -3.07 —2.83 -2.59 -—2.34 -2.10 -—1.86 -—1.62 -—1.388 —1.13 —0.85 —0.62 —0.57 0 
0.4 —3.33 -3.15 -—-2.68 -—2.19 -—1.69 -—1.21 -—0.70 -0.21 +0.28 +0.80 +41.38 +1.86 +2.03 0 
0.6 -—3.52 -—3.26 -2.54 -—-1.80 -1.04 -—0.30 +0.46 +1.23 +2.01 2.81 +3.69 +4.47 +4.81 0 
0.8 -—3.74 -—3.38 -2.42 -—1.42 -0.40 +0.60 +1.63 2.71 +3.79 +4.90 +6.09 +7.21 +7. 0 
1.0 —3.97 —3.52 -—2.31 -—1.06 +0.23 +1.51 2.84 +4.21 +5.62 +7.06 +8.55 +10.1 +11.0 ¥y 
1.2 —4.23 -3.68 -—2.22 -—0.71 +0.86 2.44 +4.08 +5.77 +7.49 +9.28 411.2 +13.2 +14.6 1. 
1.4 -—4.51 -3.86 -—2.14 -—0.37 +1.48 +3.36 +5.32 +7.36 +9.43 +11.6 +14.0 +16.5 +18.4 1. 
1.6 -—4.81 -—4.05 -—2.09 —0.04 +2.09 +4.28 +6.58 +8.98 +11.4 +14.0 +16.9 +19.9 +22.4 1. 
1.8 -—5.12 -—4.26 -—2.05 +0.27 +2.69 +5.20 +7.84 +10.6 +13.5 +16.6 +20.0 +23.5 +26.7 1, 
2.0 —5.45 -—4.50 -—-2.03 +0.56 +3.28 +6.11 +9.11 +12.3 +15.6 +19.2 +23 .2 27.3 +31.4 2. 
2.2 -—5.80 -—4.75 -—2.02 +0.85 +3.88 +7.01 +10.4 +14.0 +17.8 +21.9 +26.6 +31.4 +36.6 2. 
2.4 -6.17 -—5.02 -—2.04 +1.12 +4.46 +7.92 +11.7 +15.8 +20.0 +24.8 +30.1 +35.8 +42.1 2. 
2.6 -6.55 -5.31 —2.07 +1.37 +5.02 +8.82 +13.0 +17.5 422.3 27.7 +33.8 +40.4 +47.9 2. 
2.8 -—6.95 -—5.62 -—2.12 +1.60 +5.56 9.71 +14.3 +19.3 . +24.7 +380.8 +837.6 +45.2 +54.1 2. 
3.0 —7.37 -—5.94 -—2.19 +1.81 +6.09 +10.6 +15.6 421.1 +27.1 +33.9 441.5 +50.4 +60.8 3. 
4.0 -9.72 —7.82 —2.82 +2.58 +8.47 +15.0 22.2 +30.4 +40.0 +651.2 +64.4 +82.4 +106 4. 
5.0 —12.5 —10.1 —3.95 +2.85 +10.4 +18.9 +28.6 440.1 +54.0 +71.5 +94.9 +128 +181 5. 
6.0 —15.8 -—12.9 -—5.57 +2.57 +11.7 22.2 +34.6 +49.7 +69.0 +94.8 +4133 +197 +332 6. 
7.0 —19.3 —16.1 —7.68 +1.75 +12.5 +25.0 +40.0 +59.0 +84.3 +121 +181 +305 +750 7 
8.0 —23.3 -—19.7 -—-10.3 +0.37 +12.6 +27.0 +44.7 +67.6 499.4 +149 +241 +493 +6,030 8. 
9.0 —27.6 —23.7 -—13.3 —1.53 +12.1 28.2 +48.4 +75.2 +114 +178 +315 +869 9. 
10.0 —32.2 —28.1 -—16.8 -—3.97 +10.9 +28.7 +51.2 +81.6 +127 +207 +402 +1,880 10. 
11.0 —37.2 -—32.8 -—20.6 —6.88 +9.11 +28.4 +53.0 +86.8 +138 +234 +499 +9,948 11.! 
12.0 —42.5 . —37.7 —24.9 -10.2 +6.77 +27.5 +53.8 +90.6 +148 +258 +600 12.1 
ease ‘ TABLE 2 f 
The Function g(d, y). (v = 0.30) 
——— : a 
4 3 5 10 15 20 25 30 35 40 45 50 55 60 
0.0 —4.20 -—4.20 -4.20 -—-4.20 -—-4.20 -4.20 -4.20 -—4.20 -—4.20 —4.20 —4.20 —4.20 —4.20 
0.2 —4.42 -—4.38 -—4.26 -—4.14 -—4.02 -3.92 -3.77 -—3.65 -—3.52 -—3.41 —3.27 —3.14 —3.01 
0.4 —4.65 -—4.58 -—4.34 -4.09 -3.84 -3.64 -3.34 -—3.09 -2.84 -—2.59 —2.31 —2.03 —1.76 
0.6 —4.89 -—4.78 -—-4.42 -4.05 -3.68 -—3.35° —2.92 -2.53 -2.13 -—1.74 —1.30 —0.88 —0.44 
0.8 —5.14 -4.98 -—4.52 -4.02 -3.52 -3.05 -—2.50 -—-1.96 -—1.41 —0.86 —0.26 +0.34 +0.96 
10 —5.40 —5.20 —4.62 -4.00 —3.38 -2.74 -2.08 -—1.39 —0.68 +0.05 +0.83 41.62 42.46 
12 -5.67 -5.42 —4.74 -—3.99 -3.25 -2.46 -1.66 —0.81 +0.07 +0.98 +1.96 2.99 +4.16 
1.4 -5.95 -5.65 -—4.86 -—4.00 -—3.12 -—2.20 -—1.24 —0.23 40.83 41.95 43.15 +4.42 46.91 
1.6 6.24 -5.89 —5.00 —4.02 -—3.01 -1.95 —0.82 +0.36 41.61 42.95 44.388 +5.93 +7.73 
1.8 —6.52 -6.14 -—5.14 -—-4.04 -2.91 -—1.70 -—0.42 +0.94 +2.40 +3.97 +5.67 +7.51 +9.65 
2.0 —6.82 -6.39 -—5.29 -4.08 -2.82 -—1.47 -0.02 41.53 +3.20 +5.03 +7.02 +9.20 +11.7 
2.2 —7.12 -—6.66 -—5.45 -—4.13 -—2.74 -—1.24 +0.38 +2.13 +4.02 +6.12 +8.43 +411.0 +14.0 
2.4 —7.43 -6.93 -—5.62 -—4.19 -—2.66 -—1.02 +0.77 +2.72 +4.85 +7.23 +9.90 +13.0 +16.5 
2.6 —7.74 -—7.21 -—5.80 -—4.26 -—2.60 -—0.82 +1.15 +3.31 +5.70 +8.38 +11.4 +15.0 +19.1 
2.8 -—8.06 -—-7.49 -—5.98 -—-4.34 -2.56 -—-0.62 +1.52 +3.90 +6.55 +4+9.57 +13.0 +17.1 +21.8 
3.0 —8.38 -—7.78 -—-6.18 —4.48 -2.52 -—0.44 +1.88 +4.48 +7.41 +10.8 +14.7 +19.3 +24.8 
4.0 -—10.1 —9.32 -—7.29 -—5.04 -—2.53 +0.381 +3.57 +7.36 +11.9 +17.3 +24.2 +33.0 +45.1 : 
5.0 -—11.9 -—11.0 -8.59 -—85.90 2.83 +0.74 +4.96 +10.1 +16.4 +24.7 +36.0 +52.5 +79.2 BF 
6.0 -—13.8 -—12.8 -—10.1 —6.98 -3.40 +0.84 +6.01 +12.5 +21.0 +32.8 +850.6 +81.1 +147 
7.0 —15.7 -—14.6 -—11.7 -—8.25 -4.24 +0.62 +6.66 +14.5 425.3 +41.4 +68.6 +126 +332 . 
8.0 —17.8 —16.6 —13.4 -9.70 —5.32 +0.06 +6.90 +16.1 +29.2 +50.2 +490.3 +202  +2,680 ae 
9.0 —19.8 -—18.6 —15.2 -—11.3 -6.61 -—0.79 +6.74 +17.1 +382.4 +58.8 +116 +352 curve 
10.0 —22.0 -—20.7 -—17.2 —13.1 —8.09 -—1.91 +6.18 +17.5 +34.9 +66.6 +145 +752 ti 
11.0 —24.1 —22.8 -19.2 -14.9 -9.74 —3.28 45.25 +17.4 +36.6 +73.2 +177 +3,920 gin 
12.0 -—26.38 -—24.9 —21.2 -—16.8 -—11.5 —4.87 +3.99 +16.7 +37.3 +78.1 +207 interf 
on th 
: aes i ; cedur 
ables \ and y. These functions have been tabulated m. This circumstance requires that several determina- 
for the case in which Poisson’s ratio is equal to 0.30 tions of \ be performed with different integral values of 
(Tables 1, 2, and 3). m until the value is found that furnishes a minimum to ene 
In any specific problem the numerical values of the A. The number of trials may be reduced by a good 
quantities A, B, and « are determined by Eq. (2). It intuitive estimate of the wave length of the buckled Thi 
may here be noted that the areaw is introduced in the form. In a specific trial calculation, the numerical restra 
theory only for the determination of the net load value of y is known, and therefore there is essentially Le., K 
wo, that exists in the flange. Hence, Eq. (4) also only one unknown in Eq. (4)—namely, ». Thus, Eq. condit 
furnishes the solution in cases where the plate and the (4) may be symbolically written F(A) = 0. long, | 
flange are unequally stressed tf the parameter B is modified The nature of,the function F(A) is shown by Fig. 2. theref 
by a factor that represents the ratio of the flange stress to In general, F(A) is a decreasing function that is positive * Mi 
the plate stress. for small positive values of \ but which approaches a which 1 
In general, the quantity y is not precisely determined vertical asymptote as \ increases. Values of \ that lie Milne-” 
by Eq. (2) because it depends upon an unknown integer beyond the asymptote have no practical significance, Ed., p. 
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TABLE 3 


The Function h(A, y). (» = 0.30) 
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35 40 45 50 55 60 
.00 —3.00 -—3.00 —3.00 —3.00 —3.00 —3.00 
97 —2.95 -—2.92 —-2.91 —2.89 —2.87 —2.84 
.93 -—2.90 -—-2.85 —2.81 —2.77 —2.72 —2.68 
90 -2.84 -2.77 -—2.72 —2.65 —2.58 —2.51 
87 -—2.79 -—2.70 —2.61 —2.52 —2.42 —2.33 
84 -—2.74 -2.62 —2.51 —2.39 —2.26 —2.13 
81 -—2.69 -—2.54 —2.40 —2.25 —2.09 —1.91 
78 -—2.64 -—2.46 —2.29 —2.11 —1.90 —1.69 
75 —2.58 -—2.38 -—2.18 —1.96 —1.71 —1.45 
72 —2.53 -—2.30 —2.06 —1.80 —1.51 —1.20 
70 —2.48 -—2.22 —1.94 —1.64 —1.30 —0.93 
68 -—2.48 -2.13 —1.82 —1.47 —1.07 —0.62 
66 —2.38 -—2.05- —1.69 —1.30 —0.83 —0.29 
64 -—2.32 -—-1.96 -—1.56 —1.11 —0.58 +0.05 
62 -—2.27 -—1.88 —1.43 —0.92 —0.31 +0.41 
60 —2.22 -—-1.79 —1.30 —0.72 —0.03 +0.80 
53 —2.00 -—1.36 —0.58 +0.41 +1.69 +3.46 
52 —-1.82 -—0.95 +0.20 +1.78 +4.12 +7.93 
55 —1.69 —0.56 +1.03 +3.46 +7.65 +16.7 
63 —1.63 -—0.23 +1.88 +5.48 +13.1 +41.0 
76 —1.62 +0.03 2.71 +7.88 +22.4 +845 
93 -—-1.68 +0.20 +3.47 +10.7 +40.4 
14 —1.81 +0.27 +4.10 +13.7 +88.2 
38 -—2.00 +0.23 +4.54 +16.8 +464 
65 —2.24 +0.09 +4.77 +19.8 





























3 5 10 15 20 25 
0.0 -—3.00 -—3.00 -—3.00 -—3.00 -3.00 -—3.00 — 
0.2 -3.07 -—3.06 -3.05 —3.03 -—3.01 -2.99 — 
0.4 -—-3.14 -3.12 -3.09 -—3.06 -—3.02 -2.97 — 
0.6 -—3.20 -3.19 -—3.14 -3.08 -3.03 -—2.96 — 
0.8 -—38.27 -—3.25 -3.18 -—3.11 -—3.04 -2.95 — 
1.0 —3.34 -3.31 -—3.23 -3.14 -3.05 -2.94 — 
1.2 -—3.41 -—3.37 -3.28 -3.17 -—3.06 -—2.93 — 
1.4 -3.48 -—3.44 -—-3.33 -3.20 -3.08. —2.93 — 
1.6 -—-3.55 -3.50 -3.38 -—3.24 -3.09 -—2.92 — 
1.8 -—3.62 -3.57 -—3.48 -—3.27 -3.11 -2.92 — 
2.0 -3.69 -—3.63 -—3.48 -—3.30 -—3.12 -2.92 — 
2.2 -3.76 -3.69 -3.53 --3.34 -—3.14 -2.92 — 
2.4 -—3.83 -3.76 -3.58 -3.38 -—3.16 -—-2.92 — 
2.6 -—3.90 -3.82 -3.64 -3.41 -3.18 -2.92 — 
2.8 -—-3.97 -3.89 -3.69 -—3.45 -—3.20 -—2.92 — 
3.0 —4.04 -3.95 -—-3.74 -—3.49 -—3.22 -2.93 — 
4.0 —4.38 -—4.28 -4.01 -—3.71 -—3.37 -2.98 — 
5.0 —4.71 -—4.60 -4.30 -3.95 -—3.55 -—3.08 — 
6.0 —5.04 -4.92 -4.59 -4.21 -—3.76 -—3.22 — 
7.0 —5.36 -—5.23 -—4.88 -—4.48 -—3.99 -3.39 — 
8.0 -—5.67 -—-5.54 -—5.18 -4.75 —4.24 -3.59 — 
9.0 -—5.97 -—5.84 —5.47 -5.04 -4.50 -3.82 — 
10.0 —6.26 -6.13 -—5.76 -—5.32 -—4.77 -—4.07 — 
11.0 -—6.55 -6.41 -6.04 -5.60 -—5.05 -—4.34 — 
12.0 —6.82 -6.69 -6.32 -—5.88 -—-5.33 -4.62 — 
Fa) 
! 
! 
' 
1 
1 
' 
! 
! 
! 
H 
TO A 
! 
! 
' 
! 
! 
! 
! 
' 
' 
' 
! 
' 
! 
! 
Fic. 2. Curve showing character of the function F(A). 


since the object is to find the first intersection of the 
curve with the \ axis. In nearly all cases this intersec- 
tion may be located with sufficient accuracy by a linear 
interpolation* between two computed points that lie 
on the curve on opposite sides of the \ axis. The pro- 
cedure will be illustrated by some elementary examples. 


EXAMPLE 1 
All Edges Simply Supported 


This condition is realized if the elastic rotational 
restraint is zero and the stiffness of the flange is infinite; 
Le.,k = OandA = ~, Kq. (4) is consistent with these 
conditions only if h = 0. If the plate becomes infinitely 
long, the sequence of values of - becomes dense, and 
therefore \ is the smallest value for which h can vanish 





* Milne-Thompson has proposed a refinement of this method 
which rests upon Newton’s general interpolation formula. See 
Milne-Thompson, L. M., The Calculus of Finite Differences, 1st 
Ed.,p.4; The Macmillan & Co. Ltd., London, 1933. 


without restriction of y. Inspection of Table 3 shows 
that this value of \ lies between 35 and 40, and, by linear 
interpolation, it is found that X = 39.3; whence, by 
Eq. (3), K = 3.60. 

If the plate is short, then Eq. (2) restricts y to a set 
of discrete values. For example, if a/b = 1.5 andm = 
1, Eq. (2) yields y = 4.38. With this value of y, the 
following values are interpolated from Table 3: 


(A = 45, h = —0.28), (A = 50, h = 0.93) 


Hence, by interpolation, it is found that h vanishes 
when \ = 46.1, and it follows from Eq. (3) that K = 
4.22. Logically it should be shown that a larger value 
is obtained with m = 2, but it is intuitively evident 
that a plate with aspect ratio 1.5 buckles into a single 
lobe. 


EXAMPLE 2 
One Edge Simply Supported and the Other Edge Free 


This condition is realized if the elastic rotational re- 
straint is zero and the flange is absent; i.e., x = 0 and 
A = B = 0. Then Eq. (4) becomes F(A) = (A + 
g)/h = 0, whence \ + g = 0. Inspection of Table 2 
shows that the smaller the value of y, the smaller is 
the value of that satisfies this equation. Hence, a 
plate with these edge conditions always buckles into a 
single lobe (m = 1). If y is unrestricted (i.e., if the 
plate is infinitely long) the smallest solution of the 
equation A + g = 0 is evidently \ = 4.20; whence, 
K = 0.38. 

On the other hand, if, for example, a/b = 1, Eq. (2) 
yields y = 9.87, and the solution of the equation \ + 
g = O evidently lies between 10 and 15. Interpolation 
of data from Table 2 yields: 


(A = 10,A +g = —6.9), (A = 15, +g = 2.1) 
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Hence, by interpolation, the expression \ + g is found 


to vanish when A = 13.8. It follows that K = 1.27. 


EXAMPLE 3 
One Edge Clamped and the Other Edge Free 


This condition is realized if the elastic rotational re- 
straint is infinite and the flange is absent; i.e.,x = © 
andA = B=0. Then Eq. (4) becomes F(A) = —f/y = 
0; whence, f = 0. It is seen from Table 1 that 
when y is unrestricted, the smallest value of \ which 
satisfies this equation lies between 10 and 15. Inter- 
polation of data from the table yields \ = 12.6. Thus, 
for an infinitely long plate, K = 1.15. 

If, instead of an infinitely long plate, one considers a 
plate with aspect ratio a/b equal to 7, the admissible 
values of y are the squares of the integers; i.e., y = 
m*, For the values y = 1, 4, and 9, the following solu- 
tions of the equation f = 0 are obtained by interpola- 
tion of data from Table 1: 


Mi = 14.1, de = 12.6, A3 = 15.6 


Thus, the smallest buckling stress is obtained with a 
two-loop mode for which the coefficient \ is 12.6. Since 
this is the same coefficient that was obtained for an in- 
finitely long plate, it may be inferred that an extremely 
long plate breaks up into lobes whose complete wave 
lengths are approximately times the width of the plate. 


APPROXIMATE ALGEBRAIC SOLUTION 


If the data of Tables 1, 2, and 3 are plotted, it is 
seen that the functions f, g, and 4 are approximately 
linear in the range y < 1. Coefficients determined by 
least-square adjustment of the data corresponding to 
y = 1 furnish the following approximate expressions 
for the functions f and g: 


f = —3.00 — 2.107 + 0.270Ay7 (5) 


g = —4.20 — 1.907 + 0.139A7 
The function h may be roughly approximated by the 
first term in its MacLaurin expansion; i.e.,h = —3.00. 


These linear formulas for the functions f, g, and h de- 
termine, with Eq. (4), a linear equation in \ which has 
the following solution: 


x = Ale +3) + 1.90] y + (2.10k + 4.20) + 3x/y 6) 
1+ Bix + 3) + 0.1397 + 0.270« 


If the plate is short, the parameter y is limited to a 
few discrete values, and the value that furnishes a 
minimum to A is most quickly determined by trial. 
On the other hand, if the plate is long, the range of y 
may be considered to be a continuum, and the minimum 
value of \ is obtained when y has approximately the 
following value: 


y = V3x/[A(x + 3) + 1.90] . (7) 


This result is obtained by ignoring the term 0.139 in 
the denominator of Eq. (6) and then equating the 
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derivative to zero. Observe that y is determined by 
Eq. (2) rather than by Eq. (7) if the plate is short. 
Eqs. (6) and (7) are not always sufficiently accurate 
for engineering calculations, but, when y is small (say 
y < 2.00), these equations form a useful adjunct to the 
numerical methods described in the preceding section, 


EXAMPLE 
Zee Stringer 


As an example of the application of the preceding 
equations, it is proposed to ascertain if the outstanding 
flange of a 24S-T84 Zee-type stringer (Fig. 3) is stable 
when the stringer is subjected to 50,000 Ibs. per sq.in, 
stress. An experimentally determined buckling curve 
for this material shows that the required stability exists 
only if the parameter }/ (tA/K) is less than 11.0. It 
seems reasonable and consistent with methods that 
have been successfully utilized for hat-stringer design’ 
to measure the width b between the centroid of the e- 
bow and the centroid of the lip. This dimension is indi- 
cated on Fig. 3 to be 1.00 in. The thickness of the 
flange is 0.080 in., and, accordingly, the criterion 
b/(t~/K) < 11.0 signifies that K > 1.29. 

The flange of the stringer is itself a plate having the 
following flange (lip) area and moment of inertia; 
w = 0.0286 and J = 0.000151. Hence, by Eq. (2),4 = 
3.23 and B = 0.358. 

If coupling within the web of the stringer is disre- 
garded, the elbow may be considered to furnish elastic 
rotational restraint that is equal to the couple per unit 
length which must be applied at the elbow to produce 
one radian rotation. By the theory of cylindrical bend- 
ing of plates, this is equal to D/1.5, since the web is 
11/.in. wide. Hence, by Eq. (2), x = 0.667. 

With the foregoing constants, Eqs. (7) and (6) yield 
y = 0.38 and \ = 6.32. Since vy is small, it may be 
expected that this result is fairly accurate. 

More refined calculations may be based upon Eg. 
(4) and Tables 1, 2, and 3. The algebraic results indi- 
cate that y lies in the neighborhood of the value 0.40 
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Cross section of Zee stringer. 
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and that A lies between 5and 10. Hence, for any speci- 
fied value of y, the zero of F(A) may be determined by a 
linear interpolation between F(5) and F(10). Thus, the 
following values are obtained: 








+, O23 04 0.6 
» 7.12 6.138 6.52 


0.42 
6.16 


0.38 
6.15 





These results indicate that the minimum value of X is 
approximately 6.13, and it follows that K = 0.56. 
Since it has been shown that the value K = 1.29 is re- 
quired for stability, the flange is evidently improperly 
proportioned. 

It is interesting to observe that if the lip is eliminated 
from the flange (i.e, A = B = 0), then a numerical 
analysis similar to the preceding shows that K is in- 
creased to 0.61. This demonstration that a lip is not 
invariably a stabilizing agent does not appear paradoxi- 
cal if one considers that a lip in itself is a column with 
little inherent stability and that, consequently, an ac- 
tion may result in which the tendency for the lip to 
buckle actually reduces the stability of the combination. 

No general inferences concerning the structural use- 
fulness of a lip should be drawn from the present ex- 
ample, but the numerical methods that have been illus- 
trated should serve to determine, in any specific case, 
whether or not a lip increases the stringer crippling 
stress. It should be observed that, even though a lip 
slightly reduces the crippling stress, it may effectively 
increase section properties that control the flexural or 
torsional buckling loads. 


Appendix 


MATHEMATICAL THEORY 


Since Timoshenko! has expounded the differential- 
equation theory of buckling of axially loaded rectangu- 





lar plates, it suffices here to outline the course that has 
been followed. 

The mathematical formulation of the problem re- 
quires the establishment of conditions under which the 
differential equation 


My 
m2 


7? 


2 
os* Weeee + on? Wetnn > baum = — We (8) 


has a nontrivial solution that conforms to the boundary 
conditions 


(y/m*) we: + Voy = 0 (9) 
on the edges § = O and é = -; 
Wan ok (vy/m*) wee = KW, (10) 


on the edge 7 = 0; 
Y Ay? Bhy 
Way + (2 — ») yt ete —o Weeee - =a ee (11) 


on the edge n = 1. 

In the derivation of Eq. (11), the torsional stiffness 
of the flange has been neglected. 

It may be shown that the deflection function w must 
have the form 

w = f(n) sin még (12) 

With Eq. (8) this relation yields an ordinary differential 
equation in 7 that has the general solution 


f(n) = Ge™" + Cre" + C3 cos Bn + Cysin Bn (13) 
in which 


a=VVivty 8=VVa—7 (1) 


The boundary conditions—Egs. (9), (10), and (11)—furnish, with Eqs. (12) and (13), four linear, homogeneous, 
algebraic equations that admit nonzero values for the constants C;, C2, C3, Cy if, and only if, the determinant of 


the equations vanishes. 


The expansion of the determinant then yields a finite equation that can be expressed in 


the form of Eq. (4) when the functions f, g, and # have the following definitions: 


fy [A — (1 — v)*y] + [A + (1 — »)?4] cosh a cos 8 + y[A(1 — 2%) — y(1 — »)4J 


\ 


sinh a sin 8 


6B 








ae 


(sinh a/a) cos 8 — cosh a (sin 6/8) 
(a? — vy)*(sinh a/a) cos B — (8? + vy)? cosh a (sin B/8) 





eee is 


y[(sinh a/a) cos 8 — cosh a (sin 6/8)] 
2Vry (sinh a/a)(sin 8/8) 


(15) 





= (sinh a/a) cos 8B — cosh a(sin 8/8) 





Tabular values of these functions for the case v = 0.30 have been given in Tables 1, 2, and 3. 


REFERENCES 


' Timoshenko, S., Theory of Elastic Stability, 1st Ed., pp. 337- 
350 and Fig. 185; McGraw-Hill Book Company, Inc., New York, 
1936, 


2 Langhaar, H. L., Design of Hat-T ype Plate-Stringer Combina- 


tions, Automotive and Aviation Industries, Vol. 91, No. 11, pp. 
28-32, December 1, 1944. 











Aileron Reversal Speed by Influence 
Coefficients and Matrix Iteration 


WILLIAM T. THOMSON* 
Ryan Aeronautical Company 


ABSTRACT 


_ Methods most frequently used in determining the aileron 
reversal speed are briefly reviewed. An alternative method of 
determining the reversal speed, making use of influence coeffi- 
cients and matrix iteration, is presented. 


INTRODUCTION 


| es OF AILERON EFFECTIVENESS in developing roll 
results from the torsional flexibility of the wing. 
When an aileron is deflected downward, the increase in 
lift of the wing is partially compensated by the twisting 
of the wing resulting from the nose-down pitching 
torque. Since the torsional stiffness of the wing is 
independent of the speed, while the twisting torque 
increases as the square of the speed, it is evident that 
there is a critical speed at which the effect of the aileron 
will be completely canceled by the twisting of the wing, 
thereby resulting in a zero rolling moment. Beyond 
this critical speed the effect of the aileron will be the 
reverse of that for normal operation and, hence, the 
name ‘‘aileron reversal speed.” 

The aileron reversal speed is used as one of the design 
criteria of an airplane. In order to ensure sufficient 
aileron control at high speeds, present specifications re- 
quire that the reversal speed exceed the maximum dive 
speed by at least 15 per cent. The reversal speed is 
also used in calculating the rolling performance of the 


airplane. 


METHODS OF CALCULATION 


Several methods are available for calculating the 
aileron reversal speed. The British use a method 
based on the semirigid theory developed by Ptgsley' 
and Cox.? The equation for the reversal speed de- 
veloped by Pugsley may be expressed as follows: 


Ve =K V mp/*/epsc? 


where mz, atid ¢ are the torsional stiffness and chord at 
the mid-aileron section; s is the semispan; p is the air 
density; and K is a constant depending on the aero- 
dynamic coefficients determined by the wing-aileron 
configuration. Both the two- and three-dimensional 
cases are covered by the proper evaluation of K. For 
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the three-dimensional case, Hirst’s* simplified method 
can be used for evaluating K. 

In the United States, a method presented by Shor- 
nick* of the Army Air Forces has been adopted. In 
this method the wing is divided into a number of rigid 
spanwise segments joined elastically. By assuming a 
mode of wing twist for the condition of zero rolling 
moment, the integral set up for the twisting of the wing 
is solved by tabular method and successive approxima- 
tion. The aerodynamit coefficients used in Shornick’s 
tabular form are for open-gap ailerons hinged at the 
leading edge, and modifications are necessary for aero- 
dynamically balanced ailerons with seal. 

In the following section a third method making use of 
influence coefficients and matrix iteration is suggested 
by the author. The procedure is analogous to Biot’s® 
method of determining the wing divergence speed, and 
the equations reduce to that of divergence when 8 = 0. 
Both the open-gap aileron hinged at the leading edge 
and the aerodynamically balanced ailerons with seal 
are taken care of by the proper evaluation of the aero- 
dynamic coefficients. These aerodynamic coefficients 
should be obtained from wind-tunnel data whenever 
possible; however, theoretical expressions reduced 
from Theodorsen® and Kussner’ are presented in the 
Appendix. 


DEVELOPMENT OF EQUATIONS 


Consider a wing divided into m spanwise segments. 
Letting 8 be the aileron deflection and a the change in 
the torsional angle due to 8, the increment of lift and 
pitching torque of any section produced by 8 can be 
written as 

lL, = '/opV°c(aza; + b,B) Ax, ~ (1) 

T; = ‘/2pV°c,?(ara; + b78) Ax, (2) 
where a and 6 are aerodynamic coefficients depending 
on the wing-aileron configuration, and c; the average 
chord of the 7th section. 

At the reversal speed the rolling moment obtained by 


summing up the lift times lever arm over both wings 
must vanish 


Mr = 2 >, Lex; Ax, = 0 (3) 


#=1 


which results in the following relation 
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n s” 
Br = — (x2 Cxaryx, Axs)/ (bp Do Cx, Ax;) (4) 
be s=s 
Eq. (4) assumes 8 to be constant over the aileron span, 
and, since 8 = 0 over the rest of the wing, the summa- 
tion indicated in the denominator extends only over the 
aileron span. 
Substituting Eq. (4) into Eq. (2), the twisting torque 
of the 7th section at reversal may be expressed as 


n 
T; = '/op VP lara, = G > cx, Ax,} Ax; (5) 
t=1 
where 


s” 


G = azbr/by D> cx; Ax; (6) 
i=s’ 
a constant depending only on the wing-aileron con- 
figuration. 

The change in the angle of twist at any section k can 
be expressed in terms of the torsional influence coeffi- 
cients ¢,; defined as the angle of twist at k due to a unit 
twisting torque at 7. 


ay, = Tidm + Tobe + .--Trdin = 2d Titi: (7) 
Substituting Eq. (5) in Eq. (7) 


n n 
ao = 1/sp V2 Zz c2{ara, a G De Creeses Ars} des Arg (8) 


‘= 


By expanding and rearranging, the above equation may 
be presented in the form 


n 
ar = */2pV? Y Ques (9) 
i= 
where 
n 
Qei = Arli2by, Ax, — Ge, Ax, i Crd, Ax; (10) 
t=1 


Eqs. (9) and (10) may be interpreted as follows. 
When 8 = 0, the second term of Eq. (10) vanishes and 
Eq. (9) reduces to that of divergence given by 


n 
a, = '/opV? p> Orly? i AX; (11) 
je 
The second term of Eq. (10) is therefore the contribu- 
tion due to the aileron deflection which alters the 
divergence speed determined from the set of equations 
expressed by Eq. (11). . 

The equations represented by Eq. (9) will be satis- 
fied by n values of V: however, only the lowest root 
corresponding to the lowest speed is of practical inter- 
est. 


SOLUTION BY MATRIX ITERATION 


The solution of the equations indicated by Eq. (9) 
can be conveniently carried out by matrix iteration. 
Assuming a unit angle at the tip section, the m equations 
may be arranged in matrix form as follows 


| ax, QiQr . a Oral ay lay) 
| a2 (Qe1Qe2 . 2 - Qen| a | |) 


ws i} | .epen 


| =]. | (12) 


|| 


-Qna||1.0} jan! 


< oe 
On| QniQne- 


where the operations to be performed are indicated by’ 
Eq. (9). Assuming any reasonable twist distribution 
with unit tip value, the column matrix obtained at the 
right side is normalized and compared to that on the 
left side. Using the new deflection mode so obtained, 
the process is repeated until the modes stabi- 
lize. Generally, three or four repetitions are suffi- 
cient because of the rapid convergence of the proc- 
ess. 

The determination of the critical reversal speed is 
obvious by examination of Eq. (12). Taking the re- 
sults of the final iteration 


_ 
1/spV? | 


1 1/spV? = Un 
the reversal speed is determined from 
Ve = V2/pdn (13) 


The matrix process will always result in the lowest 
mode, since the convergence to higher modes is possible 
only by additional constraint of the orthogonality 
principle. 


COMMENTS ON NUMERICAL COMPUTATION 


The numerical computations may be carried out by 
the following procedure. Number the wing sections, 
starting with 1.0 at the root and ending with m at the 
tip. Because of the reciprocity theorem, the influence 
coefficients are symmetrical about the diagonal of the 
matrix—i.e., d:; = oi. Also, o;; will be the same for all 
values of 7 greater than k, since the same torque is trans- 
mitted to section k when the torque is applied outboard 
of k. Itis evident then that only » different values of ¢ 
appear in the matrix, all of which may be determined 
from a single test by applying the torque at the wing 
tip. Numerical values for Q,; may then be determined 
by arranging the numbers in tabular form corresponding 
to the position in the matrix equation. Having set up 
the matrix equation, the iteration process may be 
carried out in a relatively short time with a computing 
machine. 

The reversal speed is always lower than the diver- 
gence speed, since by is negative for the usual wing- 
aileron configuration. For aerodynamically balanced 
ailerons with seal, the coefficient by is considerably 
smaller than that of the open-gap aileron hinged at the 
leading edge, and- the reversal speed approaches the 
divergence speed as by approaches zero. The diver- 
gence speed can be determined from Eq. (11) in matrix 
form. 
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APPENDIX 


Theoretical expressions for the aerodynamic coeffi- 
cients of Eqs. (1) and (2) are given for two types of 
ailerons. 


Type I Aileron 
Open gap hinged at leading edge 


a, = 2x 

by = 2Ti0 

Ar = @(1/2 + a) 

br = Ty(*/2 + a) — [(Ts + Tw) /2] 


Type II Aileron 


Aerodynamically balanced with seal 


a, = 2r 
by = 2[Tw — (c — e)¢gis] 
ar = 1(1/2 + a) 
1 
br = [Ti — (¢ — e)¢s] € + s) = 
[e+ T+ on] 
where 


e = distance from mid-chord of wing to leading edge 
of aileron as a fraction of the semichord of 
wing 

distance from mid-chord of wing to elastic axis 
of wing as a fraction of the semichord of wing 
(positive if elastic axis is aft of mid-chord) 
¢ = distance from mid-chord of wing to aileron hinge 

as fraction of semichord of wing 
These coefficients may be calculated by referring to 
Figs. 1, 2, and 3. 


a 
lI 
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A Method for the Stress Analysis of 
Helicopter Blades 


DANIEL O. DOMMASCH* 


Civil Aeronautics Administration 


ABSTRACT 


The method of analysis presented herein is applicable to both 
fixed and hinged rotors. Although the equations developed ap- 
ply directly to rectangular blades, the method may be extended 
to cover blades having other shapes.‘ 

Up to this time, no ‘‘exact” general solution has been found 
for even the simplest form of the differential equation governing 
rotor bending moments, although solutions can be determined 
for certain combinations of constants (these solutions are not 
general). Because of the lack of a general exact solution, vari- 
ous approximate methods such as series, perturbation methods, 
etc., were tried, but convergence was poor in most cases, neces- 
sitating the use of an unwieldy number of terms. 

The solution presented in this paper is based on the method of 
finite differences and depends on the evaluation of a number of 
simultaneous linear equations. The method leads to a rapid 
determination of the bending moment at any point on the blade. 


SYMBOLS 


= Young’s modulus (Ibs. per sq.in.) 
= moment of Inertia (in.*) 

= cross-section Area of Blade (sq.in.) 
= unit weight (Ibs. per in.*) 

= 386 in. per sec.? 

= rotational velocity (rad. per sec.) 
centrifugal force 

differential mass 

= cone angle 

= axis of rotation 

= axis normal to Z 

= undeflected blade axis 

= axis normal to é 

= length of blade 

= axial load 

shear 

moment 

running coordinates 

(Aw/g)w* cos 3, (Ibs. per sq.in.) 


1 /w, B sin 3 P 
BS bh. J , (1/in.4 
a(e 6 ) Bi? 


Ma = hinge damping moment 
= Bcos #/EI, (1/in.*) 

tip loading (Ibs. per in.) 
(é/S) 

dg/dé 

boundary conditions 
w,/12EI S?, (1/in.5) 

B sin #/6EI, (1/in.*) 
3CS*4 

DS*/2 

4FS* 

38GS4 

value of function at a 
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fz’ = first derivative of f relative to x 
n = S' nAw/EI 
n = load factor 


INTRODUCTION 


— EQUATIONS DEVELOPED BELOW hold for blades 
of rectangular plan form; the development for 
other plan forms follow the same general outline. At 
this point it is noted that the assumption normally 
made in beam theory that EJd’¢/dt? = M, is also made 
here. This assumption implies small deflections and 
negligible lengthening of the blade under load. Be- 
cause of the small deflection notion, the angle for freely 
hinged blades must be chosen carefully (see Eq. 24). 

In applying this method, certain preliminary work is 
required, namely: 

(1) Determine the blade principal axes. 

(2) Resolve the air load relative to these axes (the 
total thrust may be assumed to cause bending about 
the principal axis of minimum inertia). 

In line with the assumption generally made in pro- 
peller analysis, only the air load normal to the prin- 
cipal axis of minimum inertia need be assumed as act- 
ing.? 


Derivation of the Differential Equation 


Referring to Fig. 1 and terminology: 






lected Blade Axis, 


Undeflected (§$) - 
i 9% 


at 
Pla & & 
Xe) Abr ‘O, ! 


/?.. 


~— 








Fic. 1. Relation of axes. 
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dm = (Aw/g)dé 
dC.F. = rw*dm 
dC.F. = (Aw/g)rwdé 


Let B = (Aw/g)w? cos 3, then, since r = £ cos J, we 
have: 
dC.F. = Bédé 
At any point &, the total C.F. is 
F.gé = BSS tdé = B/2(e — S*) (1) 
where the negative value of C.F. is due to the point S 
being the origin of integration. 
Moment Due To C.F. ABout ANY POINT 


To determine the moment about the point (£, ¢) due 
to C.F., the C.F. may be broken into components 
normal to and parallel to the £ axis. (See Fig. 2.) 


s 


f 





Deflected Blade Ax/s~, 




















yOrigin for U___ 
(6) 








(a) 
' 
i 


§ 


u and v coordinates. 








Fic. 3. 


Considering the normal (relative to £) component of 
cr.. 


(S; o)o = B sin 8 ff tdé 
(M..o), = Bsin d Jy fi tdédé 


Sz oo = (B sin 3/2)(# — S?) (2) 
(M; 0)» = (B sin 8/6)(2S* + & — 35S%E) 3) 


Considering now the parallel component of C.F. 


P = axial load = C.F. cos 3? = (B cos 8/2)(£? — S?) 
(4) 
Moment about (&, ¢) of dC.F. cos # is (see Fig. 3): 
dM, = dC.F. cos #(v) 
but in terms of u and v, dC.F. = Budu, hence: 
dM, = B cos 3 vudu 
M, = Bcos 3 f¢ vudu (5) 


MoMENT DUE To AIR LOAD 


The moment due to a known air-load distribution may 
readily be found. Depending on the down-wash dis- 
tribution, blade twist, etc., the air loading may vary 
from an approximately triangular distribution to some 
nth degree parabolic distribution or may be defined by 
a trigonometric or power series. Two simple cases are 
considered here: (1) triangular and (2) second degree 
parabolic. 


(1) Triangular Distribution (See Fig. 4) 


If tip load = w,, then w;, = (w,/S)é, and 


2 (&2 — S$?) (6) 


yA 


Ma = af It idédt = T+ (25° + # — 35%) (7) 


(2) Parabolic Distribution 


UW 
Ss = > 


2 

Str 

A 

Q, 
wt 

a 

ra) 
—~ 


Again let w, = tip load, then w;, = (w,/S*)&, and 
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Fic. 4. Air loadings considered. 
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€ 
Ws 97% 
&= 5 MS 


M==a oh fi Edédt =e (3S4 + £4 — 48%) (7a) 


ToTaAL MOMENT 


ie — oo (6a) 


(1) Triangular Air-Load Distribution 


Ws B sin 3 





) (2S? + & — 3S°§) + 


B cos 0 foudu (8) 
s 


Since M/EI = (d*¢/dé*), we have 





oe ..4. w, Bsind O¢8 3 202 
Toa mite ar) ee 
é 
——e vudu (9) 
Let 
_1l/w Bsnd\ - 
c= (Xs as) (20) 
Bcos # 
ine lhies (11) 
Then: 
Ht _ cross + #8 — 352 ‘ 
age ~ COS + #- 35%) +D f oudu (12) 


(2) Parabolic Air-Load Distribution 
M = —(3) + (5) + (7a) 











B sin 8 
M = ‘ ‘— 453 
ogi (35! + & — 45%) — FAX x 
€ 
(2.S% + & — 3S?) + Bcos off vudu—- (8a) 
s 
a* Ws 
B sind Beosd [fF 
ae. ae Fos ae t . ‘O- 
GET (2S* + 3S?) + ET < vudu (Qa) 
Let 
B sin 3 
F = w,/(12EIS?); = 
w,/( ); G SET (10a) 
Then: 
ad? 
wa = POS! + & — 45%) — 


€ 
G(2.S% + § — 387) + D vA vudu (12a) 
Ss 


TRANSFORMATION 


Prior to reduction of Eqs. (12) and (12a), note that 
the B.C. are as follows: 
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(A) Fixed Blades 
§(0) = ¢,'(0) = §"(0) = 


(B) Hinged Blades 
¢(0) = ¢,"(0) = &"(S) = 0 
no hinge damping. 
¢(0) = ¢,"(S) = 0; &"(0) = M,/EI 
with hinge damping. 
Eqs. (12) and (12a) may be simplified by considering 


the following :* 


— [fou == 4-3 (vu)du + [2 (D815 [u(S)é IS 


Now: v = ¢(u) — ¢(£). Hence, o(&) = ¢(£) — ¢()) = 
0; also, dS/dé = 0, since S is constant. 
Since u does not depend on &, 
(d/dé)(vu)du = udu(dv/dé) 
dv/d& = d[g(u) — §(&)]/d& = —df/dé 


Hence, 


€ 
@ fondu = - FE f'udu = - 
dg Ss 


Making use of Eq. (13) and differentiating Eqs. (12) 
and (12a), one obtains: 


-(4) ace - $9 
(13)* 


OT awvss . sit 
_> = — S$?) — D/2 SF 14 
of = sce — 5) — D/AE- YZ (14) 
at . ~— . 
—~ = 4F($ — S*) — 3G(é — S*) — 
de (é ) (g ) 
p/2e — sy (14a) 
dé 
Let 
g = dt/dé (15) 
Then: 
oe + Die — Se = aCe — $4 (16) 
Te + D/2(e — So = AF(E — S*) —3G(F — *) 
ag* 
(16a) 
Let 
t{//S =x (17) 
Note: 
(dp/dé)(dt/dx) = de/dx = S(dg/dé) 
(d/dx)S(de/dt) = (d/dx)(dg/dx) = d*g/dx?* 
But 


S(d/dx)(dg/dé) = S*(d*e/dé?) 


* This equation, when multiplied by B cos 3, gives the equiva- 
lent shear due to the centrifugal load component parallel to the 


E axis. 
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Hence, 
(d*p/dx*)(1/S*) = d*e/dé? (18) 
Let 
e” = d*o/dx* (19) 
Then Eqs. (16) and (16a) become 
ge” + (DS*/2)(x? — 1)e = 38CS*(x? — 1) (20) 


e” + (DS*/2)(x* — 1)e = 4FS°(x* — 1) — 
3GS*(x? — 1) (20a) 


Finally, let 
B = DS*4/2, y= 4FS (21) 
a = 3CS‘4, X\ = 3GS* 
And Eggs. (20) and (20a) become 
ge” + B(x? — 1)e = a(x? — 1) (22) 


e” + B(x? — 1)p = y(x* — 1) — Aw? — 1) (22a) 
The B.C. on ¢(x) are: 


(A) Fixed Blades 
9(0) = ¢’(1) = 0 


(B) Hinged Blades 
g'(0) = ¢’(1) = 0 
no hinge moment. 
¢’(0) = SM,/EI; ¢’(1) = 0 
with hinge damping. 
Note that 


a 
iid’ hed tm F* 
Also 

go = d/di; § = Sie dt = Sf\*o dx 


At this point it is noted that the angle # is given for 
fixed blades as the initial dihedral angle (if any); for 
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freely hinged blades, # is determined by the condition 
that the moment at the root must be zero. The angle, 
8, which satisfies this condition is given by Eqs. (24) 
and (24a). 


(1) Triangular Loading 
sin 28 = 2gw,/SwAw? (24) 


(2) Parabolic Loading 


sin 28 = */,(w,g/SwAw?) (24a) 


It is interesting to note that Eq. (24) automatically 
gives zero bending over the entire blade. This follows 
from the fact that the assumed air-load distribution in 
this case is geometrically similar to the distribution of 
centrifugal load and it is possible for the blade to align 
itself with the total load resultant without flexure. 


DEVELOPMENT OF FINITE DIFFERENCE EguarTions? 


A number of methods were tried in an effort to deter- 
mine an “exact’’ solution to Eqs. (22) and (22a), 
However, except for the fictitious case of 8 = —1 
(implies negative dihedral, very stiff rotor), none was 
found. Series solutions were found but the convergence 
was generally considered too poor in the working range 
of 8 to make the series solution of practical interest. 

As a consequence, the method of finite differences was 
tried, and it was found that this method gave accurate 
results with minimum labor. 

As an aid to those not acquainted with difference solu- 
tions of differential equations, the development of the 
difference equations is outlined below. (See Fig. 5.) 

The slope at the point ¢, for any continuous function 
¢(x) may be presented by 


Pn. nd (Ag/ Ax), = (Pn41 


The rate of change of slope, A*g/ Ax?, is then ex- 
pressed: 


— Y—1)/2Ax (25) 








TABLE 1 
1 2 3 4 5 6 7 8 
Point x x? x?—1 @ X Ax? -a X©® X 108 “BX® @ —-—2 
0 0 0 — 1.0000 —0.00250000 5.00000 0. 1750000 —2.1750000 
1 0.05 0.0025 —0.9975 —0.00249375 4.98750 0.1745625 —2.1745625 
2 0.10 0.0100 —0.9900 —0.00247500 4.95000 0.1732500 —2.1732500 
3 0.15 0.0225 —0.9775 —0.00244375 4.88750 0.1710625 —2.1710625 
4 0.20 0.0400 —0.9600 —0.00240000 4.80000 0. 1680000 —2. 1680000 
5 0.25 0.0625 —0.9375 —0.00234375 4.68750 0. 1640625 —2.1640625 
6 0.30 0.0900 —0.9100 —0.00227500 4.55000 0.1592500 —2.1592500 
7 0.35 0.1225 —0.8775 —0.00219375 4.38750 0. 1535625 —2. 1535625 
8 0.40 0.1600 —0.8400 —0.00210000 4.20000 0. 1470000 —2.1470000 
9 0.45 0.2025 —0.7975 —0.00199375 3.98750 0. 1395625 —2.1395625 
10 0.50 0.2500 —0.7500 —0.00187500 3.75000 0. 1312500 —2.1312500 
11 0.55 0.3025 —0.6975. —0.00174375 3.48750 0. 1220625 —2. 1220625 
12 0.60 0.3600 —0.6400 —0.00160000 3.20000 0. 1120000 —2.1120000 
13 0.65 0.4225 —0.5775 —0.00144875 2.88750 0.1010625 —2.1010625 
14 0.70 0.4900 —0.5100 —0.00127500 2.55000 0.0892500 —2.0892500 
15 0.75 0.5625 —0.4875 —0.00109375 2.18750 0.0765625 —2.0765625 
16 0.80 0.6400 —0.3600 > —0.00090000 1.80000 0.0630000 —2.0630000 
17 0.85 0.7225 —0.2775 —0.00069375 1.38750 0.0485625 —2.0485625 
18 0.90 0.8100 —0.1900 —0.00047500 0.95000 0.0332500 —2.0332500 
19 0.95 0.9025 —0.0975 —0.00024375 0.4875 0.0170625 —2.0170625 
20 1.00 1.0000 0.0000 0.00000000 0 0 0 
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Xe, Xn Xow 
Fie. 5. 


Ao Pn+1 — Yn Gn —~ Pn-1 
‘= = — xA 
¥e (ss ), ( Ax Ax y/ 


(A’y/ Ax*)n = (Gari — 2¢n + ¢n—1)/ Ax? (26) 


As Ax — 0, Eqs. (25) and (26) approach the exact 
derivatives; however, for a continuous function with- 
out singularity, even large values of Ax give excellent 
approximations to the exact solution. 

Noting that y” = A*y/ Ax*, Eqs. (22) and (22a) may 
be written: 


¢ati + Pn-1 + ¢n[B Ax?(x? — 1) pon 2] mid a Ax?(x? ried 1) 
(27) 








Outi + On—-1 + On[B Ax?(x? — 1) — 2] = 


y Ax?(x? — 1) — XAx%(x® — 1) (27a) 


The method of solving these equations is illustrated 
below by means of three related problems. 


ILLUSTRATIVE EXAMPLES 


Consider the following rotor blade. 


S = 20ft. = 240 in. 

chord = 2 ft. = 24 in. 

thickness (maximum) = 15 per cent C = 3.6 in. = d 

A! = 0.74 bh = 0.74(3.6)(24) = 63.94 sq.in. 

I’ = 0.0472 bh? = 0.0472(24)(3.6)3 = 52.85 in.4 

Solid wood construction, E = 1,500,000 Ibs. per sq.in. 

w, = 700 lbs. per ft. = 58.2 Ibs. per in. (assumed) 

w = 0.60(62.4) = 37.44 lbs. per ft. = 0.0216 lbs. per 
in.? 

g = 32.2(12) = 386 in. per sec.? 

w,/S = 58.2/240 = 0.242 Ibs. per sq.in. (may vary 
from 0 to 0.25 approximately) 

# = 300 r.p.m./60 (27) = 30 rad. per sec. 

w* = 900 

sin 28(0 to 0.5); #8, 0° to 15° 


(63.94) (0.0216) 


7 = 9 
Sab (000)(1) = 3.22 





B = 4* 3 cos 8 2 
g 


Ibs. per sq.in. 
C = (1/EI)[(w,/6S) — (B sin 8/6)] 


‘<. hityeins : 
EI ~ (1.5)(5.285) * 107° = 0.1261 X 10 








Since w,/S may vary from 0 to 0.25 and B sin # from. 
0 to 0.483 (approximately 0.5), 


C ranges from +0.53 X 10~* to —1.1 x 10-, 1/in.* 


pam ani = 3.22(1)(0.1261)(10-) = 40.6 x 107 


a = 3CS*; S* = (2.40)* X 108 = 33.20 X 108 
a = +1.76 to —3.65 


a 
DS* _ 40.6 x 10” X 3.32 & 10° = 67.4 


dag, hae 2 
rs a 





For a hinged blade: 
= arc sin (0.242/3.22) = arc sin 0.0752 = 4° 19’ 


Example 1 


It is desired to determine the bending moment in a 
fixed blade having 6 = 70, a = 2, and loaded by a tri- 
angular air-load distribution. 

The quantities [a( Ax?) (x? — 1)] and [8(Ax*) X 
(x? — 1)— 2] are calculated below for Ax = 0.05. 

From Table 1 and Eq. (27), the following set of linear 
equations is obtained. 


(1) ge +0 — 2.174563¢, = —0.00498750 
(2) vs +1 — 2.173250¢. = —0.00495000 
(3) ge +2 — 2.171063¢; = —0.00488750 
t 
Ny 
sil 
SX 
ef 
4 £0 
SX Z for AX=.0S. 
—4*—s— forAX=.10 
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Fic. 6. Curves showing relation between ¢ and x and ¢’ and x. 
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(4) os + ¢3 — 2.168000¢, = —0.00480000 (12) gis = —95.613207¢, + 0.8960256 + 
(5) os +4 — 2.164063¢; = —0.00468750 2.112000(136.815544¢, — 1.2941609) — 0.0032000 
(6) oz +5 — 2.159250¢. = —0.00455000 gis = 193.3412229, — 1.8404422 
(7) ¢s + — 2.153563¢; = —0.00438750 (13) gu = ~ 196. 8155449, + 1.2941609 + 
(8) v9 +7 — 2.147000¢; = —0.00442000 2.101063(193.341222y, — 1.8404422) — 0.0028875 
(9) gi + ¢s — 2.139563¢9 = —0.00398750 gu = 269.406554¢, — 2.5756116 
(10) gu + ¢9 — 2.131250¢10 = —0.00375000 (14) gs = —193.341222¢9, + 1.8404422 + 
(11) ¢2 + gio — 2.122063¢, = —0.00348750 2.089250(269.406554¢, — 2.5756116) — 0.0025500 
(12) gis + gu — 2.112000¢12 = —0.00320000 ¢is = 369.516421¢, — 3.5432043 
(13) gu + ¢2 — 2.101063¢:3 = —0.00288750 (15) os = ~ 269.406554¢1 +2 5756116 + 
(14) gis + ¢13 — 2.089250¢4 = —0.00255000 2.076563 369.5164219, — 3.5432043) — 0.00218750 
(15) ge + gu — 2.076563¢15 = —0.00218750 gis = 497.917574¢, — 4. 7842 2629 
(16) giz + ¢15 — 2.063000¢16 = —0.00180000 m, giz = —369.516421g, + 3.5432043 + 
(17) gis + gis — 2.048563¢17 = —0.00138750 063000(497.917574¢, — 4. 7849629) — 0.0018000 
(18) ¢i9 + giz — 2.033250¢:13 = —0.00095000 ¢17 = 657.6875340, — 6.3285301 
(19) goo + gis — 2.017063¢19 = —0.00048750 (17) gis = —497.917574¢, + 4.7842629 + 


(20) ie Yi9 — 2420 = 0 

The B.C. are 9(0) = 0; o’ = 0 
(0) =0—> go = 0 
g’(1) = 0 (¢a 


2 = $19 


Ax = 0 


= ¢19) 2 


Eq. (20) gives 2¢i9 — 2¢20 = 0; or gio = gro. 
The above system of equations is easily solved as 


follows: 


(1) ¢. = 2.174563¢, — 0.0049875 
(2) os = —g + 2.173250(2.174563¢, — 
0.0049875) — 0.0049500 
gs = 3.725869¢, — 0.0157891 
(3) gs = —2.174563¢, + 0.0049875 + 
2.171063(3.725869¢, — 0.0157891) — 0.00488750 
vs = 5.914514¢, — 0.0341791 
(4) gs = —3.725869¢, + 0.0157891 + 
2.168000(5.914514¢, — 0.0341791) — 0.0048000 
vs = 9.096797¢, — 0.0631112 
(5) gs = —5.914514y, + 0.0341791 + 
2.164063(9.096797¢, — 0.0631112) — 0.0046875 
gs = 13.771528¢, — 0.107085 
(6) v7: = —9.096797¢: + 0.0631112 + 
2.159250(13.771528¢, — 0.107085) — 0.0045500 
¢7 = 20.639375¢1 — 0.1726621 
(7) vs = —13.771528¢; + 0.107085 + 
2.153563(20.639375¢, — 0.1726621) — 0.0043875 
¢s = 30.676666¢, — 0.2691412 


6 


(8) v9 = —20.639375¢; + 0.1726621 + 
2.147000(30.676666¢, — 0.2691412) — 0.0042000 
yy = 45.223427¢, — 0.4093841 
(9) gio = —30.676666¢, + 0.2691412 + 


2.139563 (45.223427¢, — 0.4093841) — 0.0039875 
gio = 66.081705¢, — 0.6107494 
(10) vu = —45.223427¢, + 0.4093841 + 
2.131250(66.081705¢, — 0.6107494) — oe 
¢y = 95.613207¢, — 0.8960256 
(11) ¢2 = —66.081705¢; + 0.6107494 + 
2.122063(95.613207¢, — 0.8960256) — 0.00348750 
¢2 = 136.815544¢, — 1.2941609 


2.048563 (657.687534¢, — 6.3285301) — 0.0013875 
gig = 849.3967749, — 8.1815172 

(18) ¢i9 = —657.687534¢, + 6.3285301 + 

2.033250(849.396774¢, — 8.1815172) — 0.00095000 
gig = 1069.348457¢, — 10.3074897 

(19) 849.396774¢, — 8.1815172 — 

1.017063(1069.348457¢9, — 10.3074897) + 
0.0004875 = 0 


2.30233¢ 
= 3023367 _ 0.0096656 


- 23 8.197976 
Note: It can be shown that for this type of loading 
¢1 is a direct function of a, and the following can be 
written. 


g:/a = 0.0048328 (8 = 70) 


Further, g at any point is linearly dependent on a, 
with the consequence that ¢’ and MM are also linearly 
dependent on a. It follows that, knowing J for one 
a, M can be determined directly for any other a. 

Recapituiating: 


1 = 0.0096656 


g 
¢: = 2.174563(0.0096656) — 0.0049875 = 0.016031 
¢s = 3.725869(0.0096656) — 0.0157891 = 0.020224 
gs = 5.914514(0.0096656) — 0.0341791 = 0.022988 
¢s = 9.096797(0.0096656) — 0.0631112 = 0.024815 
¢s = 13.771528(0.0096656) — 0.107085 = 0.026025 


20.639375(0.0096656) — 0.1726621 = 0.026830 
vs = 30.676666(0.0096656) — 0.2691412 = 0.027367 
vy = 45.223427(0.0096656) — 0.4093841 = 0.027728 
¢10 = 66.081705(0.0096656) — 0.6107494 = 0.027970 
gu = 95.613207(0.0096656) — 0.8960256 = 0.028133 
v2 = 136.815544(0.0096656) — 1.2941609 = 


6 
S 
Il 


0.028243 
vis = 193.341222(0.0096656) — 1.8404422 = 
0.028317 
vu = 269.406554(0.0096656) — 2.5756116 = 
0.028364 
vs = 369.516421(0.0096656) — 3.5432043 = 
0.028394 


¢16 = 497.917574(0.0096656) — 4.7842629 = 0.028409 
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ie oc 
l 2 3 4 

Point x x? x? — 1 G 

0 0 0 —1 - 

1 0.1 0.01 —0.99 - 

2 0.2 0.04 —0.96 _ 

3 0.3 0.09 —0.91 _— 

4 0.4 0.16 —(0).84 _ 

5 0.5 0.25 —0.75 — 

6 0.6 0.36 —0.64 — 

7 0.7 0.49 —0.51 _ 

8 0.8 0.64 —0.36 _ 

9 0.9 0.81 —0.19 —_ 
10 1.0 1.00 0 





gi7 = 657.687534(0.0096656) — 6.3285301 = 

0.028415 

gis = 849.396774(0.0096656) — 8.1815172 = 

0.028412 
gig = ¢20 = 1069.348457(0.0096656) — 

10.3074897 = 0.028405 


The moments at the stations considered will be de- 
termined later; however, in order to determine the ef- 
fect of larger Ax, ¢ will be determined first for Av = 0.1. 


Example 2 


Same as Example 1 except for Ax = 0.1. 
From Table 2, the following system of equations results: 


(1) g + ¢0 — 2.693¢, = —0.0198 
(2) vs +41 — 2.672¢. = —0.0192 
(3) os + 2 — 2.637¢3 = —0.0182 
(4) os + ¢3 — 2.588¢, = —0.0168 
(5) gs + os — 2.52595 = —0.0150 
(6) or + ¢5 — 2.448¢6 = —0.0128 
(7) a%8 + ge — 2.35707 = —0.0102 
(8) oo + ¢7 — 2.25293 = —0.0072 
(9) gio + vs — 2.133¢9 = —0.0038 
(10) gnu + ¢9 — 2.000¢i0 = —0.0000 


The B.C. give go = O and gu = ¢9 = ¢i0 = 0. 
Solving: 


(1) g. = 2.693¢, — 0.0198 
(2) v3 = —¢g: + 2.672(2.693¢, — 0.0198) — 0.0192 


3 = 6.195696¢, — 0.0721056 
(3) gs = —2.693¢: + 0.0198 + 

2.637(6.195696g, — 0.0721056) — 0.0182 
1 = 13.645050¢, — 0.1885425 


6 


? 
(4) gs = —6.195696¢, + 0.0721056 + 
2.588(13.645050¢, — 0.1885425) — 0.0168 
¢s = 29.117693¢, — 0.4326424 


(5) gs = 13.64505¢: + 0.1885425 + 
2.525(29.117693¢1 — 0.4326424) — 0.0150 
vs = 59.877125¢, — 0.9188796 
(6) o, = —29.117693¢, + 0.4326424 + 


2.448(59.877125¢1 — 0.9188796) — 0.0128 
¢7 = 117.461509¢, — 1.8295749 
(7) gs = —59.877125¢; + 0.9188796 + 
2.357(117.461509¢, — 1.8295749) — 0.0102 
s = 216.979652¢, — 3.4036284 


6 


TABLE 2 
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5 6 4 Ss 

) *& Ax? aXG BXG 7) —2 
0.100 —0.0200 —0.7000 —2 7000 
0.0099 —0.0198 —0.6930 —2. 6930 
0.0096 —0.0192 —0.6720 —2.6720 
0.0091 —0.0182 —0.6370 — 2.6370 
0.0084 —0.0168 —0. 5880 —2. 5880 
0.0075 —0.0150 —0.5250 —2. 5250 
0.0064 —0.0128 —().4480 —2.4480 
0.0051 —0.0102 —0.3570 —2.3570 
0.0036 —0.0072 —(0).2520 — 2.2520 
0.0019 —0.0038 —0.1330 —2.1330 
0.0000 0.0000 +0.0000 —2 0000 

(8) go = —117.461509¢, + 1.8295749 + 
2.252(216.9796529, — 3.4036284) — 0.0072 
g9 = 371.176667¢—, — 5.8425963 


(9) 216.979652y, — 3.4036284 — 
1.133(371.176667¢, — 5.8425963) + 0.0038 = 0 
_ 3.2198332 


203.563512 


1 = 0.0158173; 


*! _ 900790865 
Qa 


Recapitulating: 
= 0.0158173 


> 
¢g. = 2.693(0.0158173) — 0.0198 = 0.02280 
gs = 6.195696(0.0158173) — 0.0721056 = 0.02589 


4 = 13.645050(0.0158173) — 0.1885425 = 0.02729 
29.117693(0.0158173) — 0.4326424 = 0.02792 


6 


eS = 
¢g6 = 59.877125(0.0158178) — 0.9188796 = 0.02821 
¢7 = 117.461509(0.0158173) — 1.8295749 = 0.02835 


¢s = 216.979652(0.0158173) — 3.4036284 = 0.02840 
¢9 = 371.176667(0.0158173) — 5.8425963 = 0.02842 
Yio = Yo = 0.02842 


Comparing the above with the results obtained for 
Ax = 0.05, one notes that there is a maximum error of 
11/3 per cent in ¢ caused by an increase in Ax from 0.05 
to 0.1 (considerably less work is involved in using 
Ax = 0.1). 

Note: g, for Ax = 0.1 = gy for Ax = 0.05, etc. 

A plot of ¢ versus x for Av = 0.05 and Ax = 0.1 is 
given below. 

Since the two curves for ¢ fall almost on top of one 
another, g’ = SM/EI is determined for Ax = 0.05 
only and plotted in Fig. 6. 

Eq. 25 gives 

On = (Gnt1 — ¢n-1)/2 Ax 
go’ is taken directly from the plot of ¢. 

.004 p ke 
= — = 0.222(check go’= “* = 0.0096656+ 
0.018 Ax 

0.05 = 0.193312) 


/ 


yy. = 








2 — oo - 0.01603 . 
gf «22s 0.016031 _ 9 16031 
2Ax 0.1 
ae 020224 — 0.009665 
sal _ gi _ 0.020 0.00 656 _ 9 10558 
2 Ax 0.1 
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a A Gs — Ga _ 0.022988 — 0.016031 = 0.06957 
2 Ax 0.1 
sb a 9092 
aisle 5 Gs 0.024815 — 0.020224 = 0.04591 
2 Ax 0.1 
oe 96095 — 99 
jiptg 6 G4 0.026025 0.022988 = 0.03037 
2 Ax 0.1 
a 26830 — 0.024 
fad ¢7 Ps _ 0.026830 0.024815 = 0.02015 
2 Ax 0.1 
_ 27: _ 26025 
oe $38 oo 0.027367 0.026025 = 0.01342 
2 Ax 0.1 
aie 97793 — 26 
a Qs ag’ 0.027728 0.026830 = 0.00898 
2 Ax 0.1 
a 97970) — 927 
page gio — 8 _ 0.027970 0.027367 = 0.00603 
2 Ax : 0.1 
— 28133 — 2772 
siden gu — 99 _ 0.028133 0.027728 = 0.00405 
2 Ax 0.1 
ae 92942 — 
eas giz — $10 _ 0.028243 — 0.027970 = 0.00273 
2 Ax 0.1 
_ 25 = 25 
eat « S22, 0.028317 — 0.028133 _ 6 goi¢4 
2 Ax 0.1 
.028364 — 0.0282: 
ous! Pe 14 P12 _ 0.0 8 6 0.0: 8 43 nes 0.0012 
2 Ax 0.1 
- 0.028394 — 0.028: 
ou) = SM, 0.028317 _ 9.00077 
2 Ax 0.1 
, vi — ¢14 0.028409 — 0.028364 " 
= = => 0.0004. 
sis 2 Ax 0.1 P 
. 417 — ¢15 0.028415 — 0.028394 
= = = 0. 21 
= 2 Ax 0.1 ” 
gis — ¢i¢ 0.028412 — 0.028409 
gr = a= = 0 
2 Ax 0.1 
eailcne Yi9 — gir 0.028405 — 0.028415 0 
s = = = 
: 2 Ax 0.1 
° Bia g20 — Fis 0.028405 a 0.028412 Li 0 
3 2 Ax 0.1 
$20" = 0 
Example 3* 


Same as Example 2 except for freely hinged rotor. 
The system of equations is the same as that given in 
Example 2; however, the B.C. are: 9’(0) = ¢’(1) = 0; 
these give go = ¢1; ¢9 = ¢10. 


Solving: 
(1) ¢2 = +1.693¢, — 0.0198 
(2) gs = —¢1 + 2.672(1.693¢, — 0.0198) — 0.0192 
¢3 = 3.523696¢, — 0.0721056 : 
(3) gs = — 1.692¢; + 0.0198 + 0.0198 + 
2.637( 3.523696¢, — 0.0721056) — 0.0182 
gs = 7.59898649, — 0.1885425 


* Eq. (24) for this case gives a = 0; however, to determine the 
effect of an incorrect choice of 3, a is taken as 2. Note that 
a = 0, with the B.C. ¢’(0) = ¢g’(1) automatically gives ¢g =0 
at all points. 
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(4) os = —3.523696¢, + 0. 0721056 + 
2.588(7.598986¢, — 0.1885425) — 0.0163 
vs = 16.142480¢, —0.4326424 


(5) gs = —7.598986¢, + 0.1885425 + 
2.525(16.142480¢, — 0.4326424) — 0.0159 
vs = 33.160776¢, — 0.9188796 
(6) gr = —16.142480¢ + 0.4326424 + 


2.448(33.160776¢, — 0.9188796) — 0.0128 
v7 = 65.035010¢, — 1.8295749 
(7) vs = —33.160776¢, + 0.9188796 + 
2.357 (65.035010¢, — 1.8295749) — 0.0102 
vs = 120.126743¢, — 3.4036284 
(8) ¢ 9 = —65.035010¢; + 1.8295749 + 
2.252(120.126743¢, — 3.4036284) — 0.0072 
v9 = 205.490415¢, — 5.8425963 
(9) 120.126743¢, — 3.4036284 — 
1.133(205.490415¢, — 5.8425963) — 0.0038 


1 = 3.2198332/112.693897 = 0.0285715 
Y1/q = 0.0142857 


Recapitulating: 


go = g = 0.0285715 
1.693000(0.0285715) — 0.0198 = 0.0285715 


g = 
vs = 3.523696(0.0285715) — 0.0721056 = 
0.0285717 
gs = 7.598986(0.0285715) — 0.1885425 = 
0.0285719 
¢s = 16.142480(0.0285715) — 0.4326424 = 
0.028725 
¢s = 33.160776(0.0285715) — 0.9188796 = 
0.0285735 
¢7 = 65.035010(0.0285715) — 1.8295749 = 
Q.0285729 
¢s = 120.126743(0.0285715) — 3.4036234 = 
0.0285728 
gio = oo = 205.490415(0.0285715) — 5.8425963 = 
; 0.0285731 


The above result is interesting, since it indicates that 
an incorrect choice of # did not give an incorrect answer. 
In the above problem flapping inertia, coriolis and gyro- 
scopic loads were not considered. If the net distri- 
bution of these loads combined with air load were geo- 
metrically similar to the centrifugal load distribution, 
the condition of zero moment all along the blade would 
still exist even if these loads were considered. If the 
load distributions were not similar, small bending 
moments would exist, since, in this case, the blade would 
have to flex in order to satisfy the boundary conditions. 
The effect of a hinge damping moment is considered i 
the Appendix. 


CONCLUSIONS 


(1) The method of finite differences leads to a rapid 
determination of rotor bending moments. 

(2) For constant chord (rectangular plan form) 
blades of uniform section and weight, loaded with 4 
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triangular load, the bending moments are directly pro- 
portional to .a and, hence, to the tip loading. 


APPENDIX 


], Summary of Formulas for Untapered Blades 


Axial Load: 
_ Boosd,,, _ » _ Bsind —_ 2 af 
ee ree S?) aoe S) 


Shear with triangular air load: 





- w, Bsind .. ow _ Beosd,,, _ 24 
5,= ($5 — 2822) ce — sy — 282 — soft 
Shear with parabolic air load: 
_ (33 _ ox, _ Bsind 2_ $2) — 
S= 25 (@ — 9) — Ses) 
Boos? .., _ dg 
ed 


Bending Moment: 
M = (EI/S)(de/dx); 9’ = de/dx 
Deflection : 
f= Sf," odx 
Il. Effect of Blade Dead Weight 


For a rectangular blade of uniform section and weight, 
the dead-weight inertia loading is uniform and has a 
running value of: mwdW, where » = load factor and 
dW = differential blade weight; but, dW = Awdé; 
hence, ndW = nAwdt. 


S, = inertia shear = nAw fe'dt = nAw(t — S) 
tt 
nAw "4 y age = "A™ (g2 + S* — 25%) 
S y4 


The inertia moment acts in opposition to the air-load 
moment; hence Eqs. (22) and (22a) become: 


g” + B(x? — 1)p = a(x? — 1) — a(x — 1) 


M, 


(II—22)" 


e” + B(x? — 1)e = y(x* — 1) — A(x? — 1) — n(x — 1) 
(I1I—22a) 
where n = S*2Aw/EI. 


Ill. Effects of Hinge Damping 


Since the magnitude of a hinge damping moment is 
dependent on the angular flapping velocity and not 
alone on the blade position, the angle, # (and, hence, a), 
is indeterminate when hinge moments due to damping 
are present. The moment distribution can be deter- 
mined, however, if an assumed value of # is 
used. 

The method of determining the bending moment at 
any point with hinge damping follows exactly the pro- 
cedure of the illustrative examples with the exception 
of consideration of the B.C.—in this case, 

EI 


‘Ss go = Ma; go’ = SM,/EI 


Assuming a straight line between ¢; and gp: 


(¢1 — ¢o)/ Ax = SM,/EI 
ee SM, Ax/EI + 0 


Knowing the relation between ¢ and ¢), the remaining 
¢n are readily determined. 
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Influence of Different Factors on Buckling 
Loads of Curved Thin Aluminum-Alloy 


Sheets for Monocoque Constructions 
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Ecole Polytechnique de Montreal 


ABSTRACT 


Aerodynamic requirements of smooth surfaces and fair con- 
tours for both wings and fuselage have made aeronautical engi- 
neers conscious of the buckling problems associated with the de- 
sign of semimonocoque (sheet-stringer) construction. This paper 
discusses the fundamental principles involved in increasing the 
buckling strength of the sheet. Test results show a 50 to 70 per 
cent increase in ultimate buckling load obtained by changes in 
the configuration of the sheet, which resulted in more load being 
carried in the supported edges and less in the unsupported central 
region. Although several suggestions are made, the practical 
method of achieving the desired results is left to the ingenuity of 
the designer. As a sidelight, these tests emphasize the necessity 
for well-controlled end and side conditions if consistent results 
are to be obtained in panel testing. 


INTRODUCTION 


I ORDER TO THROW some additional light on the 
problem of pure monocoque construction, which is 
intimately related to the buckling of sheets in compres- 
sion, a few problems concerning the behavior of thin 
curved aluminum-alloy panels of different shape and 
width, as well as the influence of different loading condi- 
tions, were investigated. In connection with an in- 
vestigation titled ‘‘Curved Aluminum-Alloy Sheets in 
Compression for Monocoque Constructions,’’: this ad- 
ditional work was carried out on duralumin sheets 
(17ST) of 0.032-in. thickness with a radius of curvature 
of 24 in. and a free length between grips of 9in. The 
following points were investigated. 

(1) Influence of different shapes of panels on buck- 
ling loads. 

(2) Clearance between the ends of the panels and 
the loading tables. 

(3) Width of panels and width of supported edges 
in relation to buckling and maximum loads. 

(4) Influence of clearance between the supported 
edges of the sheet and the width between the teeth of 
the guides. 

The guides used had, if not otherwise stated, teeth of 
1/,-in. depth, spaced 1/2 in. apart (see Fig. 5 of refer- 
ence 1). The tests were carried out on a 60,000-lb. 
Baldwin-Southwark-Tate-Emery testing machine using 
a loading range of 12,000 Ibs. 
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INFLUENCE OF DIFFERENT SHAPES OF PANELS ON 
BUCKLING LOADS 


Thin sheets used as structural elements in stressed 
skin constructions show, in general, low buckling 
stresses compared to the ultimate loads they can carry, 
It can be supposed that the compression stress on the 
loaded edge is more or less uniformly distributed in the 
pure elastic range up to the critical buckling load. But 
after buckling takes place, the greatest part of the sheet 
(the central part), representing more than two-thirds 
of the width of the sheet, is, as already shown! for 
panels without guides, not able to carry an additional 
load. After buckling, this part of the sheet has prac- 
tically no additional load-carrying capacity. Depend- 
ing upon the stiffness of the testing machine, the maxi- 
mum buckling load suddenly falls down a certain per- 
centage (22-50 per cent) of this load and, despite this 
automatic unloading of the panels, their supported 
edges are not in the same proportion unloaded; on the 
contrary, the edge stress increases during this auto- 
matic unloading. This can be easily observed by the 
readings of the extensometers recording the edge strain. 
Furthermore, the edges of the panels cannot take part 
in the buckling process itself since they are supported 
and are forced by the guides to remain straight. 

The total deflection of the sheet takes place normally 
in two distinct parts: 

*(1) The completely elastic deformation up to the 
maximum buckling load; this is the straight stress- 
strain part AB of the usual load-deformation diagram, 
as shown in Fig. 1, diagram a. 

(2) From the buckling load to the ultimate load 
the stress-strain diagram (part CD) is curved and shows 
permanent deformations which gradually increase with 
higher loads. ' 

The airplane designer is concerned with (a) initial 
buckling and whether it is elastic or permanent and (b) 
ultimate strength. Although elastic buckling of wings 
in normal flight condition has been permitted in the 
past, there is considerable evidence that the buckles 
have produced adverse pressure gradients on the wings 
and control surfaces, both at high speed and at high 
load factor, and designers are realizing the necessity 
for eliminating buckles. 


204 





temp 
the s 





“essed 
*kling 
carry, 
n the 
n the 

But 
sheet 
hirds 
1 for 
ional 
prac- 
end- 


per- 
this 
ted 
1 the 
uto- 

the 
rain, 
part 
rted 


ally 








BUCKLING LOADS OF CURVED THIN SHEETS 











Fic. 1. Influence of shape of panels. 


If in a loaded thin sheet, a distribution of the com- 
pressive stresses could be obtained in such a way that 
the ultimate load is reached without buckling taking 
place, then a favorable solution to this basic problem 
could be found. If it were possible to transform the dis- 
continuity of the diagram during the loading process 
(Fig. 1, diagram a) into an uninterrupted line with an 
elastic range up to maximum load, this embarrassing 
problem would no longer exist for the constructor, 
whether or not buckling can take place. This would 
mean a stress-strain diagram as shown in Fig. 1, dia- 
gram 6, having possibly an elastic range A’—B’ as high 
as in point B, diagram a, for the usual panels and about 
the same maximum load in D. It is quite natural that 
the slope of the elastic part A’-B’ would in this case be 
somewhat smaller (a’ < a). It seems that this would 
represent a more elastic and stable skin construction as 
opposed to the fairly stiff and instable body materials 
used in present panels. 

As shown in Fig. 1, diagrams ¢ to h, preliminary at- 
tempts have been made to attain thisaim. By changing 
the shape, especially that of the central part of the 


to 
i=) 
or 


panels, it may be expected that some improvement in 
this respect can be obtained. About 30 panels, each 
one with a different distribution of the central longi- 
tudinal internal forces, were tested for this purpose. 
This was achieved by making a few horizontal jigsaw 
incisions of different lengths located at different spots 
in the central part of the sheets, as well as by incisions 
located at an angle of 45° to the loading direction. In . 
addition, some discharging central openings of varying 
sizes, differing in number and orientation, vertical slots 
to limit the spreading of the buckle, transverse slots, 
slots at the base and at the top of the plates, were made, 
as well as longitudinal overlapping slots, as represented 
in Figs. 2-4. 

It can be stated that by this means the shape of the 
load deflection diagram can be changed appreciably in 
the required direction, so that a more or less uninter- 
rupted compression curve can be obtained. Some 
sheets, especially those with horizontal slots, show a 
fairly straight load deflection diagram up to about 
1,500 Ibs. (for sheets of 0.032-in. thickness and 6-in. 
width, as shown in Fig. 1, c, g, kh). Some of these panels 
had only an unpronounced buckle and showed a small 
difference of slope between the pure elastic range and 
the parts of permanent deformation in the diagram 
(Schemes Ic 5 and Ic 6, Fig. 1). 


CLEARANCE BETWEEN THE ENDS OF THE PANELS AND 
THE LOADING TABLES 


Investigations of greater practical interest seem to 
be those conducted in a similar direction on normal pan- 
els without, however, appreciably changing their geo- 
metric shape. A more regular distribution of the unit 
stress over the whole width of the panels can be at- 
tained in the following manner. If the panel is fastened 
so that the unsupported central parts are not loaded 
uniformly over the whole width, that is, less intensely 





Fics. 2, 3, and 4. 
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Fic. 5. Influence of shape of panel. 


in the middle than in the restrained edges, the buckling 
will probably take place at a higher load compared to 
the normally uniformly loaded panels at their ends. 
In fact, by introducing the end loads gradually into the 
panels, starting from the two restrained edges toward 
the central part, the maximum buckling load can be 
shifted appreciably toward much higher values. With 
this in mind, the central part of the panel -was loaded 
gradually and less intensely than the supported edges. 
A clearance of only a few thousandths of an inch on 
both ends of the panels, decreasing from the center part 
toward the edges, which were loaded over a width of 
about 1'/, in., plainly had the expected effect. As 
shown in Fig. 5, with a clearance of about 0.003 in. 
only (diagram IB 13, IB 11), a buckling load of 2,220 
to 2,490 Ibs. with an ultimate load of 2,550 and 2,660 
Ibs. was reached; for a clearance of 0.01 in., this load 
was extremely regular for two panels—in one case it 
was as high as 2,630 Ibs. (panel H-41), in the other 
(panel H-42) it was 2,640 Ibs. It is interesting to note 
that this last buckling load is only 20 Ibs. lower than the 
ultimate load of 2,660 lbs. reached during the buckling 
test. By providing the panels with a still greater clear- 
ance of 0.015 in. at both ends, the buckling load was ex- 
ceptionally high; for panel H-47 it reached 2,810 Ibs., 
while the collapse load was the same (2,810 Ibs.). It 
is remarkable to note that with 0.02-in. clearance 
(panel, H-48), a buckling load could be reached which 
was even higher (about 180 lbs.) than the ultimate load 
of 2,650 lbs. For these sheets, tested under normal 
loading conditions, the average buckling load was about 
1,660 Ibs. and the collapse load about 2,700 to 2,800 Ibs. 


WIDTH OF PANELS AND WIDTH OF SUPPORTED EDGES IN 
RELATION TO BUCKLING AND Maximum LOADS 


The results of this investigation are shown in Fig. 6. 
Sheets having widths of 1 in. up to 8 in., in steps of 1 
and 2 in. were tested. The free edges of these panels 
were simply supported in the usual manner by guides of 
1/,-in. pitch and '/2-in. depth. The thickness of, the 
sheet tested (H) was 0.032 in. and the specimens were 
designated H 39a, H 39b, H 45b, H 39c, H 24 and H 40. 
As will be seen by the copies of the original records in 
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this figure, a panel of 1-in. width (H 39a), that is, with. 
out free space between the guides, and having '/;-in, 
depths of the teeth, supported a maximum load of 949 
Ibs. This diagram is composed of a completely elastic 
reversible straight part up to about 750 Ibs. without 
permanent deformations and a smooth curve of gradual 
curvature indicating that permanent deformation took 
place. The edge stress in this narrow sheet, supported 
over the whole length and width in a sort of sheath of 
1-in. width and provided with teeth separated by 1/2 in, 
has, in view of these ideally supported edges, a maxi- 
mum value. It was found for these conditions to be as 
high as 30,000 Ibs. per sq.in. 

It is curious to note that a sheet of twice this width 
(2 in.-H 39b), having a free space between the guides of 
1 in., showed a maximum load of nearly double the 
value measured for the sheet of 1l-in. width; about 
1,840 Ibs. were recorded, corresponding to an edge 
stress of 28,700 lbs. per sq.in. This means that the 
guides of '/.-in. depth can be supposed to act ona 
depth of more than '/, in. A free space of at least 
another 1/2 in. outside the guides seems to protect the 
panel from buckling in its middle part; this record may 
be considered to be straight up to a load of about 900 
Ibs., where a discontinuity in the chart took place. For 
a greater width of 3 in. (H 45b), this protection of the 
guides was already less effective; the diagram shows, at 
about 1,200 lbs., the beginning of a small buckle and 
the maximum load with 2,370 lbs. is only about 30 per 
cent higher than that for the panel of 2-in. width. A 
check test for another panel of the same width gave a 
maximum load of about 2,380 Ibs. (H 45a). A more dis- 
tinct buckle, shown in the stress-strain record, was ob- 
served with panel H 39c of 4-in. width. At 1,230 Ibs. 
the first buckle took place and a maximum load of 
2,580 Ibs. was reached. The usual shape of the stress- 
strain diagram was obtained in this test series with a 
sheet of 6-in. width (H 24) showing a pronounced buckle 
at 1,790 lbs. and an ultimate load of 2,720 lbs. Further- 
more, a sheet of 8-in. width, according to record H 40, 
had two buckles, one at 1,860 lbs. and another at 1,880 
Ibs. The ultimate load of this sheet was 2,930 Ibs. 

From the results of another test completing the pre- 
ceding investigation, it can be deduced what part of 








Influence of width of panels. 


Fic. 6. 





the t 
buck] 
6-in. 

2-, al 
porte 
Ib 3, 
Fig. “ 
shows 
Ibs. (1 
ried 2 
lbs. 

lower 


shape 
the se 


INFLI 
EDGE 


Th 
in, thi 
ness, 
+().0( 
and a 
the us 
of the 





is, with. 
y '/s-in, 
1 of 949 
r elastic 
without 
gradual 
on took 
ported 
eath of 


ates 
/2i,, 


1 maxi- 
0 be as 


> width 
lides of 
ble the 
about 
n edge 
iat the 
t ona 
t least 
ect the 
‘d may 
ut 900 
>. For 
of the 
ows, at 
le and 
30 per 
th. A 
ave a 
re dis- 
as ob- 
30 Ibs. 
vad of 
stress- 
with a 
yuckle 
rther- 
H 40, 
1,880 
S. 
€ pre- 
art of 





BUCKLING LOADS OF CURVED THIN 























Ib I Ib 2 Ib 3 
ia tt i 
Ic9 IclO Icll 














Influence of width of panels. 


Fic. 7. 


the total load the supported edges may carry without 
buckling in the central part of the sheet. Sheets of 
6-in. width, provided with a central free space of 1-, 
2-, and 4-in. width, were tested. The remaining sup- 
ported edges of these panels designated Ib 1, Ib 2, and 
Ib 3, were of 21/2-, 2-, and 1l-in. width. According to 
Fig. 7, a panel with a central free space of 4-in. width 
showed no buckle and carried a maximum load of 1,740 
lbs. (Ib 1); a panel with a free space of 2 in. (Ib 2) car- 
ried 2,450 Ibs., and with 1-in. width (Ib 3) carried 2,540 
Ibs. All three sheets did not show any buckling at 
lower loads. Check tests made with three identical 
shaped panels (Ic 9, 10, and 11) showed practically 
the same results (Fig. 7). 


INFLUENCE OF CLEARANCE BETWEEN THE SUPPORTED 
EDGES OF THE SHEETS AND THE TEETH OF THE GUIDES 


The clearance of the guiding slots for sheets of 0.020- 
in. thickness as well as for sheets up to 0.050-in. thick- 
ness, has been established by other investigators with 
+0.002 in. That means +10 per cent for thin sheets 
and about +6 per cent for thick sheets. However, if 
the usual fairly appreciable variation of the thicknesses 
of these sheets is taken into consideration, where dif- 
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ferences of about 10 per cent in a single sheet have been 
measured the clearance between the guides and the ef- 
fective thickness of the individual panels becomes 
much greater and equally of much greater importance. 
Furthermore, since commercial tolerances of +0.0025 
in. are admissible for these sheets, which may measure 
from 0.0295-in. to 0.0345-in. thickness, it is obvious 


_ that the clearance between the guides and the supported 


edges of the sheets may also vary in fairly wide limits. 

Since in the preceding test series it could be observed 
that the test results at buckling loads and at ultimate 
loads seemed somewhat influenced by the clearance of 
the guides and their adjustment, as well as by the thick- 
ness of the sheets tested, the following supplementary 
investigation has been carried out. 

The normal guides with '/.-in. pitch were adjusted 
with '/2-in. depth in three different ways. Three tests 
(N 14, 15, and 16) were made with guides clamped 
strongly on the edges of the sheets, so that not only was 
there no clearance between the sheet and the guides 
but the latter were fixed so tightly that the sheet could 
not glide in the guides. 

Fig. 8 shows stress-strain records of these three pan- 
els. The buckling loads, varying between 1,600 
and 1,640 Ibs., are fairly régular and the ultimate loads 
which these panels can carry are from 2,960 to 3,200 
Ibs. 

Four other tests with panels cut out of the same sheet 
N (39, 40, 31, 32) were carried out in the normal way. 
That is, the guides were adjusted on each panel with a 
minimum of clearance, sufficient, however, to permit 
the panel to glide easily in the guides. Since the panels 
tested varied appreciably in thickness (0.031 to 0.034 
in.) the adjustment of the guides was made individually 
for each sheet in such a way that the guides could just 
glide by their own weight on the edges of the panels. 
Apart from one panel (N 32 with 1,580 Ibs.), the buck- 
ling loads for these tests were higher than the preceding 
values. Three panels showed buckling loads of 1,720, 
1,810, and 1,990 Ibs., with an average value of 1,775 Ibs., 
that is, about 150 Ibs. higher than those of the preced- 
ing series. On the other hand, the ultimate loads, 
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varying between 2,660 and 2,850 lbs., were distinctly 
lower than those for the panels with strongly clamped 
guides. 

The third test series was made with a somewhat ex- 
aggerated clearance between the guides and the edges 
of the panel. The guides were preliminarily adjusted 
on a sheet of 0.045-iy. thickness, and afterward used on 
panels with 0.033-in. thickness, so that a clearance of 
about 0.012 in. was used during these tests. 

As may be seen by Fig. 8, the records of panels N 6, 
7, and 8 show a buckling load that is distinctly lower 
than that for panels tested with normally adjusted 
guides, as used in the preceding series (N 39, 40, 31, 
32). The buckle took place between 1,400 and 1,620 
Ibs. The average value for this test series is only 
1,480 lbs. compared with 1,775 lbs. for normally ad- 
justed guides with a maximum value of 1,990 Ibs. The 
buckling loads of sheets with exaggerated clearance are 
distinctly lower than those for strongly clamped guides 
as shown by the first series. 

An even more pronounced difference between these 
three ways of adjusting the guides on the panels could 
be observed for the ultimate loads. From an average 
ultimate load of 3,093 Ibs. for tightly clamped guides, 
this load falls down to an average of 2,770 lbs. for nor- 
mally adjusted guides and for guides adjusted with a 
clearance of 0.012 in. the average ultimate load is as 
low as 2,353 Ibs. 

From these tests and records it becomes obvious that 
the influence of the clearance of the guides supporting 
the free edges of the panels cannot be neglected, either 
in further investigations or in theoretic calculations 
concerning the maximum buckling loads of panels. It 
seems that it may become necessary to adjust the 
guides individually to each sheet and even to each edge 
of the panel, in order to obtain reliable results in com- 
pression tests of this sort. Guides with a nonadjust- 
able slot and a given tolerance do not appear to give 
the highest possible accuracy. Guides with variable 
slots, enabling them to be adjusted to the edge thick- 
ness of each individual panel, should be used in further 
investigations. 


TEST RESULTS 


The results of the preceding tests may be summed up 
as follows: 


(1) The stress-strain record of panels under com- 
pression loads can be greatly influenced by the shape 
of the panel. By changing the shape of the central 
part of the panel, a certain improvement can be ob- 
tained because of a different distribution of the longi- 
tudinal forces in the panels. Using discharging central 
openings of varying sizes and orientation, an uninter- 
rupted load deflection diagram without the typical 
buckling phenomenon was obtained. The elastic range 
of the diagram changed over smoothly to the plastic 
range without any sudden interruption in the graphs. 


(2) By testing the panels with a clearance of only 
3 to 20 thousandths of an inch at both ends betweg 
the loading tables, the buckling load could be distinctly 


- modified. By only a few thousandths of an inch clear. 


ance, the buckling load could be increased materially 
compared to panels supported over the whole width; 
with 0.003-in. clearance, the buckling load was aboy 
35 to 50 per cent higher, with 0.01 in. the buckling loag 
was about 60 per cent higher, and with 0.015 in. th 
buckling load, ameliorated about 70 per cent, reached 
the ultimate load. In a further test with a clearang 
of 0.02 in. the curicus phenomenon took place that the 
buckling load (2,830-3,890 Ibs.), registered by the re. 
corder, was higher than the collapse load (2,650 Ibs), 
It is believed that by more elaborated tests of this 
kind even better and more interesting results can be ex. 
pected. 

(3) The width of the panels is of the greatest in. 
fluence on their load-carrying capacity, as well as on 
the shape of their load deflection diagram. Because of 
the restrainers, the smaller the panels are in width, the 
more resistant they are. The collapse load of the 1- 
in. width panel was as high as 940 Ibs., reaching nearly 
30,000 Ibs. per sq.in.; no buckling effect took place and 
the proportional limit in compression was attained for 
this panel at about 23,000 lbs. per sq.in. For the pane 
of 2-in. width, the collapse load reached was nearly 
twice that of the l-in. sheet, while the proportional 
limit was only about 25 per cent higher. With a further 
increase in the width of the panels (3 and 4 in.) the col- 
lapse load increased more slowly and a distinct discon- 
tinuity in the diagram indicated that the buckling load 
was reached at 1,230 lbs. for the sheet of 4-in. width. 
Panels with a width of two-thirds of their free length 
showed the characteristic upper and lower buckling 
load, while for larger panels (8 in.) a double buckling 
effect at 1,860 and 1,880 lbs. tcok place with only a 
small increase for the collapse load. Compared to the 
4-in. panel, this load was only 15 per cent higher for an 
increase of the width of 100 per cent. 

Quite different loading diagrams were obtained with 
panels of different width by taking the middle part of 
the panels out. In all these cases, no buckling took 
place, the shape of the diagram was quite different, and 
collapse loads were in practically fairly good agreement 
with the results of the preceding series. 

(4) The tests concerning the clearance of the guides 
showed that a fairly high collapse load (average 3,09 
Ibs.) was reached with guides clamped tightly to the 
edges of the panels. On the other hand, the buckling 
load for these panels was not high (1,623 Ibs.). A bet 
ter buckling load was obtained with guides adjusted 0 
the panels without clearance; the average value wa 
about 10 per cent higher than for the preceding series 
while the collapse load was distinctly lower (2,770 lbs). 
With the same guides fastened on the panels with 4 
clearance of 0.012 in., the average buckling load was 

(Continued on page 217) 
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A Flight Investigation of the Effects of 
Surface Finish on Wing Profile Drag 


W. H. DAY* anp J. M. SCHWARZBACH* 
Curtiss-Wright Corporation 


ABSTRACT 


An investigation in flight was conducted on a Curtiss P-40F 
airplane to determine the effect on the section profile drag of 
smoothing the wing surface. A total and static pressure rake, 
located aft of the wing trailing edge in the wake of the test sec- 
tion, measured the loss in momentum from which the sectional 
profile drag coefficient (cg,) was calculated. The airfoil section 
at the plane of measurement was approximately an NACA 2212.8 
profile. For each surface finish condition, data were taken at 
values of airplane lift coefficient varying between C, = 0.15 and 
C, = 0.35. The corresponding Mach Numbers were between 
M = 0.28and M = 0.41. 

The surface finish smoothness has a large effect on the profile 
drag coefficient. At an airplane lift coefficient of C, = 0.16 
(RN = 16.5 X 108, M = 0.405), four degrees of wing surface 
smoothness produced the following values of c,,: 


Fresh camouflage paint (standard produc- 
tion type finish) 
Above paint after several days weathering 


Cap = 0.0088 
Ca, = 0.0083 


Above paint sanded and pumiced Ca, = 0.0078 
Wing surface filled and smoothed to obtain 
as aerodynamically smooth a surface as 
= 0.0064 


possible Cap 


Calculations made on the basis of these data indicate that an 
appreciable improvement in performance could be obtained 
merely by. smoothing the factory-finished wing. 

Two flights, one with a thick mud layer on the wing and the 
other with a somewhat lighter layer, indicated the necessity of 
flying with clean wings. The thin layer of mud caused a drag 
increase above the production wing of about Ac,, = 0.0035, while 
the addition of the heavy mud caused the profile drag to be 
slightly more than double that of the factory-finished wing. 
These surfaces were tested to show what might be expected in 
actual combat service. 

This investigation also indicated that calculations made with 
the simplified Silverstein-Katzoff integrating equation, corrected 
for compressibility, yielded values of profile drag coefficient sub- 
stantially equal to those calculated by use of a more exact equa- 
tion that required lengthy point-by-point integrations across the 
wake to determine c,,. 


INTRODUCTION 


R= TESTS HAVE SHOWN that over one-third of 
the total parasite drag of a modern fighter-type 
airplane is due to the wing profile drag. A series of 
measurements in flight on a representative section of 
the Curtiss P-40F wing (NACA 2212.8 airfoil) at 
moderate Reynolds and Mach numbers was recently 
made to determine whether the wing profile drag 


“coefficient, c,,, could be reduced merely by improving 
“the wing surface finish. Additional tests were made 


Received August 10, 1945. 
* Aerodynamicists, Airplane Division, Buffalo Plant. 


with the wing surface covered with two different thick- 
nesses of mud to approximate service conditions often 
encountered. 

During the investigation, a comparison was made of 
the drag coefficients calculated by Silverstein and 
Katzoff’s simplified’ and exact equations in an attempt 
to determine the accuracy of the simplified equation. 


INSTRUMENTATION 


The flight measurements were made on a representa- 
tive section of the left wing of a Curtiss P-40F airplane. 
The wing profile drag coefficients were calculated from 
measurements of the momentum loss in the wake. The 
loss in momentum was obtained from measurements of 
the static and total pressure in the wake. A rake 
mounted 0.16c aft of the wing about 108 in. outboard 
of the airplane centerline on a cantilever strut fastened 
through the upper surface to the ribs near the trailing 
edge, as shown in Fig. 1, was used to measure the pres- 
sures. The wing chord at this station was approxi- 
mately 78.15 in. The pressures were recorded by a 
photo-observer, located in the baggage compartment, 
which contained a liquid manometer, altimeter, air- 
speed indicator, and chronometer. The pilot recorded 
the free air temperatures, which were measured by a 
calibrated thermocouple located on the underside of 
the wing. 


FLIGHT PROCEDURE 


Twelve flights were made at approximately 5,000-ft. 
density altitude. The first flight was made to locate 
the rake in the center of the wake and to determine, by 
means of long light yarn streamers, the direction of air- 
flow at the rake. Small tufts on the wing surface indi- 
cated an absence of noticeable transverse flow. The 
rake was centered in the wake for subsequent flights, 
and the total and static pressure tubes were adjusted to 
be parallel to the direction of airflow. 

On each of the remaining flights, data were taken at 
stabilized, true air speeds varying from 200 to approxi- 
mately 330 m.p.h. in increments of about 20 m.p.h. 
In addition, the first and last runs during each flight 
were made at 100 m.p.h. to enable the pilot to record 
the calibrated outside air temperature. 

The flights were made in air as smooth as possible to 
minimize the acceleration effects on the manometer 
fluid, but some fluctuation of the liquid persisted. In 
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Fic. 1. Profile drag rake mounted on the Curtiss P-40F wing. 0 Carll 1 
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most cases, the amplitude was small and was ¢,, com- —_—, ; AILERON 
puted from averaged liquid levels. This averaged "ga betas ig 
liquid level was determined from the theoretical de- FLAP HINGE 
flection of the free stream static pressure below the free bait on 
stream total head liquid level which should have existed th 
to correspond to the indicated air-speed reading. The sh 
compressibility and air-speed pitot-static tube interfer- flt 
ence corrections were included in the calculations. of 
This averaged, theoretical value checked the actual 
liquid deflections used for computation, and it had been 4 
determined previously that the manometer level 
fluctuation occurred simultaneously in all of the tubes. th 
In the few cases where the liquid level fluctuation was | 
: ‘ ot 
appreciable, the data were discarded. 
The inclusion of “dampers,” consisting of porous 5 
plugs of copper or glass wool in the lower sections of the 
manometer tubes, helped to alleviate the fluctuation of - 
the liquid levels with no effect on the manometer read- a2 
ings. oe 
The photo-observer instruments were calibrated at = 
frequent intervals. Tests in the Curtiss-Buffalo 8-ft. ju 
wind tunnel showed negligible interference effects on 6 
the rake static-pressure readings from the supporting 
strut. At the start and finish of each flight the actual 
rake tube locations were measured, and these values hi 
were used to plot the corresponding recorded pressures. co 
The take-off weight and fuel consumption were meas- sal 
ured during the flight to allow calculation of airplane wi 
lift coefficient (Cz). we 
pe 
SURFACE FINISHES 
Flights were made with six different types of surface 
finish on the wing test panel. The test section of the Su 
wing is shown in Fig. 2. The first surface tested was 
the regular camouflage paint as applied to the produc- . 
2 tion airplane in the factory. This surface was modified ‘ 
in three ways. Finally, to simulate the worst condi- be 
- that igi encounters’ mm, actual combat Fic. 3. Curtiss P-40F wing test section with a thick irregular ca: 
operation from small advanced air bases, two tests were mud layer (Finish 2). 
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run with mud applied to the surface—first in a thick 
layer and then in a thin layer. The six finishes are 
described as follows: 


() Freshly Painted Production Camouflage Paint 


This was the actual production finish that was ap- 
plied to the airplanes leaving the factory. The tex- 
ture of the surface was rough and sandpapery. No 
attempt was made to fill or smooth the skin joiuts, 
ammunition door and flap piano-type hinges, rivets, or 
other protuberances. A Q-shaped wire (plane of wire 
perpendicular to the flight path), extending about 1/, in. 
above the surface at the ammunition door hinge, was 
not removed. 


(2) Production Paint with Thick Layer of Mud 


A thick irregular layer of mud was applied over Finish 
1 as shown in Fig. 3. 


(3) Weathered Production Paint 


The mud of Finish 2 was removed by water and the 
surface left to weather for several days, leaving a finish 
that was not so rough to the touch as Finish 1. The Q- 
shaped wire at the ammunition door hinge was made 
flush with the surface and remained so for the remainder 
of the program. 


44) Thin Mud Layer on Production Paint 


A light layer of mud was applied to Finish 3. It was 
thinner than Finish 2, but it did have small pebbles and 
other irregularities imbedded in it as shown in Fig. 4. 


45) Sanded ai d Pumiced Paint 


The production paint was sanded smooth with No. 
320 (wet) sandpaper followed by a rubbing with wet 
rags and pumice. There was no filling of the irregulari- 
ties or smoothing of protuberances; the paint was made 
just as smooth as possible. 


(6) Aerodynamically Smooth Surface 


All of the surface irregularities including piano 
hinges, skin joints, rivets, ejection chutes, and the 
covered gun ports were filled with glazing putty and 
sanded. The surface was pumiced with long chord- 
wise strokes until it felt perfectly smooth. Every care 
was taken to keep the surface clean and free from im- 
perfections for all flights with this particular surface. 


RESULTS AND DISCUSSION 


Surface Finish Effect on Wing Profile Drag 


The variation of section profile drag with airplane 
lift coefficient is shown in Fig. 5. The test Mach Num- 
bers and Reynolds Numbers are also presented. 

The wing with the production finish of freshly applied 
camouflaged paint (Finish 1) had a profile drag co- 





Fic. 4. Curtiss P-40F wing test section with a thin mud layer 
(Finish 4). 


efficient of about cg, = 0.0088 in the high-speed range 
(C, = 0.15 to 0.18). Allowing the paint to weather 
and making the 2 hinge wire flush with the wing surface 
(Finish 3) caused a drag decrease of about Acg, = 0.0005 
as shown by the curves in Fig. 5. This slight profile 
drag decrease with weathering is confirmed by results of 
speed-power calibration-type flight tests on several 
fighter airplanes. An examination of the two surfaces 
revealed that the weathered paint looked and felt 
somewhat smoother to the touch. 
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A further decrease of Acg, = 0.0005 was obtained by 
wet-sanding and pumicing the paint surface (Finish 5). 
The resultant paint was extremely smooth, but the un- 
filled rivet depressions, hinges, and lap joints were still 
present. A decrease in wing profile drag of about 
Aca) = 0.0010 seems to be all that would be obtained by 
merely smoothing the production finish paint. 

Filling the irregularities, fairing the surface contour, 
and then carefully painting and polishing the wing 
(Finish 6) reduced the drag coefficient to about ca, = 
0.0064. This value was checked by four separate 
flights. 

The addition of a comparatively thin coating of mud 
(Finish 4) to the wing increased the drag coefficient to 
Cao = 0.0125. A heavy, irregular mud application 
(Finish 2) more than doubled the drag of the production 
finished section. 

These results show the relation of surface finish to 
the drag coefficient for an NACA 2212.8 airfoil. It had 
been long thought that comparatively small drag de- 
creases could be obtained by smoothing a conventional- 
type airfoil. These results show, however, that a 
smooth wing is essential to the attainment of low drag, 
even on a conventional airfoil. Of course, the need for 
a smooth wing is even greater if the new N.A.C.A. low- 


drag airfoils are employed. 


Estimation of Change in High Speed Due to Smoothing 
the Wing 


As an example of the application of the flight cz, data’ 
the change in performance of a Curtiss P-40N airplane 
due to changes in wing finish is estimated. All of the 
P-40 series airplanes have identical wings; conse- 
quently, the incremental drag coefficient values ob- 
tained from these tests on an early model P-40F may be 
applied to the more recent P-40N airplane. Several 
assumptions are made in order to apply the drag 
coefficients, measured only at a single spanwise station, 
to the entire wing. 

The section profile drag coefficient increments could 
probably be applied directly to the entire area outboard 
of the guns, since, as shown in Fig. 2, there are no major 
protuberances that could not be eliminated. If only 
the wing area outboard of the guns is refinished, calcu- 
lations indicate that sanding and pumicing (Finish 5) 
will add 3.4 m.p.h. to the high speed, while application 
of the aerodynamic smooth finish (No. 6) will add an 
additional 4.3 m.p.h. or a total increase of 7.7 m.p.h. 
over the high speed obtained with the production finish. 

The parasite drag coefficient for the airplane with the 
smoothed wings was calculated from the following 


equation: 


Coy = Cdy, — (AcayS2/S) 


where 
Cp,, = the original airplane parasite drag coefficient 
Cp, = the airplane parasite drag coefficient with a 


changed wing finish 
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Sa = the wing area affected by the smoothing 
process 
S = the total wing area 


Regular methods of calculation were then used to 
estimate the high speed. 

It is more difficult to judge the drag reductions that 
may actually be attained by smoothing the entire wing, 
since the exact effect of the slipstream is unknown. If, 
because of turbulence in the slipstream, transition 
occurs at or close to the wing leading edge, then the 
effect of smoothing the surface in the turbulent boun- 
dary layer area can be estimated only roughly. The 
presence of the guns and wheel fairings and other large 
protuberances on this portion of the wing is another 
factor preventing the application of the measured sec- 
tion drag data. 

Flight tests made on the airfoil with Finish 6 and 
transition induced at 0.10c on both the upper and lower 
surfaces yielded a drag coefficient in the high-speed 
range of about cg, = 0.0077. 

Two assumptions were made in order to estimate the 
high speed increase that might be obtained by applying 
the aerodynamically smooth finish (No. 6) to the entire 
P-40 wing. 

(1) The low profile drag coefficient obtained from 
these tests of Finish No. 6 apply to the entire wing area 
outboard of the guns. 

(2) The drag coefficient obtained from the tests with 
transition fixed at 10 per cent chord is realized over the 
inboard wing area except for those areas affected by 
protuberances, wheel pockets, landing gear, etc., where 
there is no drag decrease. 

The calculated increase in high speed, then, is 9 m.p.h. 

The validity of these calculations has been approxi- 
mately verified by flight tests of various airplanes, in- 
cluding several P-40’s with the entire wing finished in 
various ways. 

The wing finish on each production airplane is not 
exactly the same. It is dependent upon the painter, 
atmospheric conditions, care taken in handling the air- 
plane, and other factors. Overspray, a condition that 
results in a rougher surface, may be present to a greater 
or lesser degree. All these variables tend to make the 
exact prediction of the gain in high speed or improve- 
ment of other performance items of a specific airplane 
difficult and prevent an exact comparison of increments 
obtained from different tests. 

The values estimated herein are of use as an indica- 
tion of the gain that might be expected on the P-40 air- 
plane, but they also illustrate the general tendency 
toward improved performance that might be obtained 
with any airplane merely by constructing smooth 
wings. It is obvious that, if the production wing 
finish is poor to start with, more of a gain will be realized 
by smoothing the wing than if a consistently good finish 
is always produced by standard shop practice. 
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Comparison of Methods 


Several equations have been derived to calculate the 
section profile drag coefficient from measurements of 
the total and static pressures in the wake and in the free 
stream. 

Jones? derived an equation that neglects the effect 
of compressibility. 


_2 ” f= Bi (1 - 
Cay c J Ay anal Po : 


His = Br) dy (1) 


Hy — po 

where 

H, = free stream total pressure 

H, = total pressure in wake 

bo = free stream static pressure 

fi = Static pressure in wake 

c = local wing chord 

J” = integrate across the wake 

dy = elemental distance across wake 





Silverstein and Katzoff! derived a simplified equation 
for calculating c,, based on the assumption of a cosine 
squared wake shape but also neglecting the effects of 
compressibility. They substitute for the point-by-point 
calculation of the integrand in the Jones equation and 
subsequent integration across the wake, a “summary”’ 
equation including a factor, F;, which is dependent on 
the maximum and average pressure losses in the 
wake. 

- Ww (Ho fir Hh) average 9 
Cao Fy - (Ho — po) (2) 
where W is the width of the wake. 

Values of F; were obtained’ by simultaneous solution 
of Eqs. (1) and (2) for:wakes having cosine squared 
shapes and various maximum and average pressure 
losses in the wake. 





They also derived an equation requiring an integration across the wake which takes into account the effects of 
compressibility, including the density decrease in the wing wake due to the temperature rise at high Mach Num- 


bers. 














w= 2 f y" Ca — PNG Ems + 02a Ch BG + a (1 + 0.202M,?) 








(Za tahy @ 
(Ho — po)(1 + m) 


where 
y = ratio of the specific heat at constant pressure to that at constant volume 
1+ mm = 1 + '/4Mo? + '/Mot + ... 
l+m=1+ 7M? + Voli + ... 
1+ m= 1+ '/4M2? + /oMe* +... 
Mo = Mach Number of the free stream 
M, = Mach Number in the wake at the measuring plane 
M, = Mach Number in the wake at an infinite distance downstream from the wing 








TABLE 1 


TABLE 2 





—Section Profile Drag Coefficient, ca— 


Calculated Calculated 
with with 
Airplane C, Eq. (1) Eq. (2) 
0.347 0.0094 0.0096 
0.290 0.0090 0.0092 
0.249 0.0092 0.0092 
0.209 0.0087 0.0087 
0.196 0.0090 0.0089 
0.159 0.0087 0.0092 





A correction for the effect of compressibility, (F./F,), 
which may be applied to drag coefficients calculated 
with Eqs. (1) or (2), was derived by Silverstein and 
Katzoff from Eq. (3) on the basis of the cosine squared 
wake shape. The corrected drag coefficient is ex- 
pressed as c 


——Section Profile Drag-——— 
Coefficient, ca, 





Calculated Calculated 
with with 
Mach Eqs. (2) Eq. 
Airplane C, Number and (3) (4) 
0.347 0.280 0.0093 0.0092 
0.290 0.306 0.0088 0.0088 
0.249 0.329 0.0088 0.0089 
0.209 0.359 0.0083 0.0084 
0.196 0.396 0.0085 0.0082 
0.159 0.408 0.0084 0.0083 
Cay = (F./F:) (incompressible cay) (4) 


At low Mach Numbers, Eqs. (1) and (2) give accurate 
Cay Values and differ only by the error in the assumption 
of a cosine squared wake shape. A comparison of the 


(Continued on page 217) 








Static Displacement and Coupled Frequencies 
of a Twisted Bar 


R. J. DUNHOLTER* 


Unwersity of Cincinnati 


INTRODUCTION 


EEASPIC_ BODY CONSIDERED is a Straight bar 
with constant reetangular cross section and an ini- 
tial structural twist of @ (eemstant) radians per unit of 
length. The object of the paper is twofold: (1) to 
calculate the static displacement components u(z) and 
v(z) in the direction of the principal axes u and v at 
each section due to the uniform loads w, and w, also in 
the direction of the principal axes; (2) to approximate 
by the energy method, using the displacement functions 
mentioned in (1), the first two natural frequencies of 
the bar. The method used is applicable to any end 
conditions but only the fixed-free case is considered. 
In the case of a fixed-free bar without twist the two 
fundamental frequencies corresponding to vibrations 
normal and parallel to the minor principal axis v are 
given by the formulas 


3.52V EI,/yL‘ 
3.52V rv El,/yL4 


where r = J,/I, and y is the mass per unit length of 
bar. In the case of the bar with an initial structural 
twist these two fundamentals will be coupled. 

The coefficients replacing those given above for the 
untwisted bar are obtained by the energy method. It 
is known that the functions representing the normal 
modes are solutions of the variational equation 


8(T — U) = 0 (1) 


where 7 and U are the kinetic and potential energies, 
respectively. If J and U are evaluated for an assumed 
displacement function satisfying the boundary condi- 
tions, then Eq. (1) furnishes approximate values for the 
natural frequencies. * A function that can be ex- 
pected to give good values for the two coupled funda- 
mental frequencies of the twisted bar is the function 
representing the static displacement due to uniform 
loads w, and w, in the direction of the principal axes 
u and v, respectively, at each section. 

The fundamental formulas from which the compo- 
nent displacement functions are obtained are 


EI, [(d*u/dz*) — 20(dv/dz) — 6u] = — M, (2) 
EI ,,[(d*v/dz?) + 20(du/dz) — 0v] = +M, (8) 


These are a special case of more general formulas for 


Wy 


Wp 


Received June 29, 1945. 
* Assistant Professor of Mathematics. 


bars with initial variable twist and initial curvature dis- 
cussed and derived in detail by A. E. H. Love’ and H. 
Reissner. For completeness a brief derivation of the 
formulas is given in the next section. 


DERIVATION OF FUNDAMENTAL FORMULAS 


The differential equation for the displacement of a 
bar may be written in the form 


Elx = M (4) 


where « is the curvature of the plane displacement curve 
and M is the bending moment. In the case of an 
untwisted bar with loads along both principal axes it 
is customary to apply Eq. (4) to the two projections of 
the space displacement curve. For a twisted bar with 
small rate of twist @ it is assumed to be a good approxi- 
mation to apply Eq. (4) to the projections of the space 
curve on the wu, z and v, z axes at each section (Fig. 1). 








Fic. 1. End view of bar showing reference axes. 


The curvatures « of the two projections at each section 
can be calculated directly or, more simply, as follows. 
Let r be the position vector of a point on the space 


displacement curve. Then 
r = u(z)i + v(2)j + ck (5) 


where i, j, R is a right-hand set of unit vectors. If s 
is distance measured along the curve and if the dis- 
placement is assumed to be small so that ds/dz = 1, 
then 


dr/dz = (dr/ds)(ds/dz) = T (6) 
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DISPLACEMENT OF A TWISTED BAR 


d’r/ds? = (dT/ds)(ds/dz) = dT/dz = «N_ (7) 


where T and N are the unit tangent and unit principal 
normal vectors, respectively. 

Before applying Eqs (6) and (7) to Eq. (5) it should 
be noted that 


di/dz = 0j and dj/dz = —06i 


Now by differentiating Eq. (5) twice with respect to z 
it is easily shown that 


du dv d*y du )i 
= — 20— — ou)i 20 — — 6% 
an (= dz . u) ind (2 + dz in ad 
(8) 








It is apparent now that the components on the right of 
Eq. (8) are the curvatures of the projections of the 
space curve on the u, 2 and », 2 planes, respectively, at 
section z. Substitution of these values in Eq. (4) gives 
the formulas (2) and (3). 


BENDING MOMENTS AND SHEARS 


M, and M, are the moments of the loads between 
section z and the free end with respect to the u and » 
axes, respectively, at section z. M, is positive if, as a 
vector, it points in the direction of the u axis. The same 
right-hand rule applies to M,. 





U aris ab free erd 


End view of bar showing axes at section z = z and loads 
Ww, and w, at section z = ¢. 


Fic. 2. 


The following formulas are obtained with reference 
to Fig. 2. 


M, = w, f,*(z — t) cos 6 (z — t) dt — 

Wy Jo (z — t) sin 6(z — 2) di 
M, = —w,/J/,"(z — t) cos 0(z — t) dt — 

w,J,*(z — #) sin 0(z — t) dt 
V, = uw, fo* sin 6(z — t) dt + wy J,* cos 0(z — t) dt 
V, = w, f,* cos 0(z — t) dt — w, f,* sin 0(z — #) dt 


bo 
— 
or 


The second derivatives with respect to z can be expressed 
in the forms 

D°M, = —@M, — 20V,+w, 

D°?M, = —#M, — 20V, — wy 

D*V, = —@V, + Ow, 

D*V, —#V, — Ow, 
where D = d/dz. 


If V, and V, are eliminated from the above expres- 
sions, one obtains M, and M, in the operational forms 


M, = —6*w,/(D? + 6)? (9) 
M, = +6°w,/(D? + 62)? (10) 


DISPLACEMENT FUNCTIONS 


The expressions for M, and M, given in Eqs. (9) and 
(10) are now substituted in the Eqs. (2) and (3). The 
differential equations of the component displacement 
functions u(z) and v(z), expressed in operational form, 
are 


(D? — @)u — 26Dv = —6@*%a/(D* + 67)? (11) 

20Du + (D* — 6?)v = —6*b/(D? + 6*)? (12) 
where a = w,/EI,, b = w,/EI,,. 

There are eight boundary conditions. The slope 


and displacement at the fixed end of the bar are zero. 
Thus, 


(1) w=0,2 
(2) v=0,2 


L (3) Du=0,2=L 
L (4) Dv =0,2=L 


Since M, and M, are zero at the free end it follows from 
Eqs. (2) apd (3) that 
(5) (D? — 6*)u — 20Dv 
(6) 206Du + (D? — 6?) 


0,2z=0 
0, z=0 


The shears V, and V, are also zero at the free end. 
Then from the formulas for M, and M, one observes 
that 
DM, = +0M,+ V,=0,2=0 
DM, = —0M, — V, =0,z=0 
The last two boundary conditions are now obtained by 
differentiating Eqs. (2) and (3) once with respect to z. 
Thus, 
(7) D(D* — 6*)u — 20D% = 0,2 =0 
(8) 26D*%u + D(D? — @)\v =0,2 = 


Eqs. (11) and (12) can now be solved. From Eq. (11) 


_ (Dt + 6)*(D? — 6%)u + 6a 
m" 26D(D? + 62)? 





Substituting this value of v in Eq. (12), one obtains 
(Dt + 6)4u = 64a (13) 


The complete solution of Eq. (13) can be expressed in 
the form 
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u = R(A + Be + Co? + Dae + 0/04 (14) 


where R indicates the real part of the expression follow- 
ing and A, B, C, D are complex constants. It is easily 
shown that corresponding to u given by Eq. (14) 


v = Ri(A + Bs + Co? + Dz)e™ + 5/04 


where, as in Eq. (14), 7 is the imaginary unit. 


(15) 





P= 7 1 — cos AL — z) + ; 6>(L — 2)? cos 2 — 0(L — z) sin 0(L — 2) + = OL — 2)?(2L + 2) sin | (18) 
Q= LY sin 6(L — z) — OL — 2) cos AL — z) + = OL — 2)? sin 6z — = ONL — 2)?(2L + 2) cos se | (19): 


If in Eqs. (16) and (17) one substitutes 6 = 0 the 
expressions P and Q have the indeterminate form 0/0. 
However, it is easily shown that 


lim P = (38L4 — 4L%z + 2*)/24 
6—>0 

lm Q=0 

6—>0 


Thus, u and v will reduce when 6 — 0 to the familiar 
formula for the displacement of a cantilever bar due to 
a uniform load. 


COUPLED FUNDAMENTAL FREQUENCIES 


The total potential energy of the oscillation, con- 

sidering bending only, is 
U = fo" (M,2/2EI,) dz + f* (M,2/2€1,) dz 
where M, and M, are determined by Eqs. (2) and (3). 
Using the displacement functions u(z) and v(z) given 
by Eqs. (16) and (17) one obtains 
M, = —EI,(P’a + Q"b — 0?Pa — 0b — 20P’b + 
200Q’a) 


M, = +EI,(P") — Q"a — 6°Pb + 6°0a + 20P’a + 
260'b) 


where the primes denote differentiation with respect to 
z. Or 


M, = —EI,(aC, + bC2) 








M, = +£I,(0C, — aCz) 
where 
C; = P” + 20Q’ — #P (20) 
C, = Q” — 26P’ — OQ (21) 
Thus, 
U= St (aC, + bC,)? dz + aS (0C, — aC)? dz 


The total kinetic energy of the oscillation is 


1946 


The eight real constants are determined by the eight 
boundary conditions given. The final solutions for y 
and v are 


u = Pa + Qb (16) 
v = Pb — Qa (17) 


where 





T = Yayot fy (ut + 02) dz 
= Vay? fr (P? + Q)(a? + b4) ds (22) 


where w is the frequency of the oscillation. The quan- 
tity (J — U) must now be minimized considering a 
and 5 as the variables. The two equations obtained 
by differentiating (JT — U) partially with respect to 
a and 6 are 


al K* C2dze+ rf" C2 dz — p> f*(P? + Q?) dz] + 
b1 — rf" GCG,dz=0 (23) 


a(l—r)fK* GC, dz + bff" C2 dz + rfp" C2 ds — 
pf (P? + Q*) dz] =0 (24) 


where p? = yw?/EI,. 

The frequency equation which is of the second: degree 
in p* is the determinant of the Eqs. (23) and (24). 
After the integrals appearing in these equations have 
been calculated the solution for p? can be expressed in 
the form ; 


p= sai +rjd+(1- im 





(25) 


4L4 gZ 
where 
d = 487r(1 + cos r) — 96 sin r + 873 (26) 
e = —96 sin r + 6(5 — 27?) sin 27 — 367 cos 27 + 


48r (1 + cos rT) — 24r (27) 
f = —33 + 48 cos rt + 247 sin r — 3(5 — 27?) cos 27 — 
187 sin 27 (28) 
g = (1,6807 + 4207? + 2175 + 11177) — (3,360 + 210r‘) X 
sin r + (1,6807 — 14073 + 7775) cos 7 (29) 

and r = OL. 


The: following coefficients of V EI,/yL‘ were cal- 
culated for r = 16: 


0; 3.53 and 14.12 
0.8; 3.68 and 14.0 


II 


abou 


clear 
Fr 
that 
mate 
slot < 
est pi 


the d 
panel 


drag ¢ 
Table 
cates 

valid 

Numtk 
than t 
The 
that is 
in Tat 


lated 
the di 








eight 
for u 


(16) 
(17) 


(18) 


(19) 


(22) 
quan- 
ring a 
ained 
ct to 


3] + 
(23) 


(24) 


legree 
(24). 
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(28) 


\r*) X 
(29) 


7 = 7/2; 3.96 and 14.0 
T= T; 5.5 and 13.6 
7= 27; 9.6 and 11.5 


The exact values for r = O are 3.52 and 14.08. 

As r — « the coefficients approach the common value 
10.4. For r = 1 Eq. (25) gives nearly the same value, 
ground 3.55, for all the above values of r. 


CONCLUSION 


The preceding analysis is assumed to be valid only 
for small values of the rate of twist 6. However, Eq. 
(25) shows that the coefficients depend only on the total 
twist 7 and, hence, large values of 7 can be considered 
provided the corresponding L is large. The effect of 
coupling due to twist is relatively small for values of + 
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as high as 90°. It should be noted that the two funda- 
mental frequencies converge to a common value as the 
twist increases; the smaller increases in value while 
the larger one decreases. 
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about 15 per cent and the average collapse load about 
20 per cent lower than those for guides fastened without 
clearance. 

From the results of these tests it becomes obvious 
that the influence of the clearance of the guides is of 
material importance. Guides with a nonadjustable 
slot and a given tolerance do not seem to give the high- 
est possible loads of the panels; it appears to be neces- 
sary to adjust the guides individually to each panel. 

Another test series, concerning the basic question of 
the detrimental influence of any eccentricity on curved 
panels, has been carried out. It could be shown that, 


Buckling Loads of Curved Thin Aluminum-Alloy Sheets 


(Continued from page 208) 


contrary to the general consideration admitting that 
the slightest eccentricity of the panels is extremely 
harmful to their mechanical resistance under axial 
thrust, a local eccentricity of the panels up to as much 
as several tenths of an inch remains without marked 
influence on the buckling load or even on the collapse 
load. Because of their special character, details about 
these test results can only be published at a later 
date. 
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drag coefficients calculated from Curtiss flight-test data, 
Table 1, as well as other unpublished comparisons, indi- 
cates that the cosine squared wake shape assumption is 
valid for normal wing wakes. At moderate Mach 
Numbers Eqs. (1) and (2) yield a cg, value that is higher 
than the actual drag coefficient. 

The drag coefficients that are calculated by methods 
that include the effects of compressibility are presented 
in Table 2. 

A comparison between the drag coefficients calcu- 
lated with Eq. (3) and those obtained by correcting 
the drag coefficients in Table 1 by the Silverstein- 








Effects of Surface Finish on Wing Profile Drag 


(Continued from page 213) 


Katzoff compressibility factor (F,/Fi) shows that the 
two compare closely. 

The conclusion may be drawn that use of the simpli- 
fied Silverstein-Katzoff integrating equation yields drag 
coefficients that compare favorably with those calcu- 
ated with the longer, exact equations. 
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Airfoil Theory for Rotary Wing Aircraft 


RUFUS 


ISAACS* 


United Aircraft Corporation 


INTRODUCTION 


I THE AUTHOR’S PREVIOUS PAPER, ‘‘Airfoil Theory for 
Flows of Variable Velocity,” hereinafter referred 
to as AT, there was developed an airfoil theory for a 
stream with a greatly varying velocity. In particular, 
the lift was explicitly given for the case where the veloc- 
ity consisted of a sinusoidal variation superposed on a 
constant velocity. 

The chief application in mind of such a theory was 
rotary wing aircraft. In the present paper the theory, 
insofar as this application goes, will be completed. 
Specifically it includes the following addenda to AT: 
The angle of attack as well as the velocity is variable and 
the aerodynamic moment is computed in addition to 
the lift. The lift and moment are given specifically for 
the case where the angle of attack is also of the form 
constant plus sinusoidal variation; the frequency being 
the same as for the velocity variation but the ampli- 
tude and phasé are arbitrary. The result is in the form 
of a complete Fourier series, as some applications call 
for a knowledge of certain higher harmonics. 

To save a great deal of needless duplication the 
chain of ideas in AT will be followed closely; most of the 
content of this paper stresses the modifications neces- 


sary. 
GENERAL THEORY 


The hypotheses of this section are the same as the 
corresponding section of AT except that a is now a func- 
tion of the time. Eq. (1) of AT then becomes 


(1’) 


The reasoning of AT is unchanged until the equating 


v = aw + xa + v7(x, t) 


to zero of the net flow through the plate. Then the 
system (10) is replaced by 
C1 = by — by — 2aw 
Co = by — bd, — ba (10’) 
Cy = Ogt1 — bg-1 (sw = 3, 4,.....) 


and likewise 


o= 2aw + ba a bo + by (12) 
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This leads to the final relation (in the periodic form), 


br(2aw + b&) = T(t) + 


fe ve) Ving 





2b 
— Wet —» 





+1-—- 1] dy 
(14 


The left side will be designated by g(t). Making the 


same transformation as in AT, 
T(t) = Q(W) 
Eq. (14’) assumes the form 


g(t) = Q(W) + SK? Q'(W — r»YIV (2b/d) +1 — 1Jr 
(16) 


The usual problem is: given w(t) and a(t), to find 
T(t). This now admits of a precise mathematical 
formulation as follows: g(t) is a known function. By 
integration W(t) is known and this may be inverted in 
general to give ¢ as a function of W. Substituting its 
value in g(t), there is obtained g as a function of VW, 
Eq. (16’) is then an integral equation with Q as the un- 
known function, whose solution, when W is replaced by 
W(t), solves the problem. 

For later purposes, there is obtained by making the 
usual transformation 


1 ers 
b=-oh aw - »| 





V/a? — 1 


(a — Lee TN 


(16") 


1 + (A/8). 


where a 


THE CASE OF ROTARY WING AIRCRAFT 


Now g(t) will be taken as a periodic function of cir- 
cular frequency w. And so for some comple 
gn(— © <m < =), g(t) may be written 


digne”” 
n 


Then, as in Eq. (17) of AT, Q(W) may 


wt, 


with + 
be written as 
yo. petal (17' 
n 
where w is the mean value of w(t). 
Inserting these series in Eq. (16’) and proceeding 4 
in AT, there results 
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AIRFOIL THEORY FOR ROTARY WING AIRCRAFT 


where A, = R,a, and 


Ro = 1 

Pe (™) = ¥(nk), (n > 0) 
Wo 

R-,= R, 


yis given by AT, Eq. (19). 


wand W are now taken as the same specific functions 
as on AT, page 114. Then applying the same tech- 
nique as in AT there is obtained 


1 ; : 
A, _ - gm ru + c sin rains sed cos * dr 
=a Mm 


Then by well-known expressions for Bessel functions 


Ay = a (o/2t)(g1 aad £1) 
. -> mgni-"Jy-m(na), (n + 0) 


It is now time to take a more explicit expression for 
a(t) such as will be encountered in connection with 
rotary wing aircraft. Thus 


ao(1 + B; sin r + B2 cos 7) 


a= 


Let 
h = 2bway = (2b)(Cz /2) 


where C;, is the mean lift coefficient, i.e., for angle of 


attack ap. 
Then from a short calculation 


fo = woh [1 + (o/2)A:] 

b= -2.= - E + (& + ia( 1 + it) 
22 = Z-2 = —(woha/4)(Bi + Be) 

gn = 0, (|n| > 2) 


It follows further 
Ao/hwo =1+ (o?/2) + o|By — Be (k/4)] 


and using the abbreviations 


J, = J,(nc) 
s* = Jn+i1(no) soca n-1(no) 
forn + 0 
A, = < woh(H, + iH’) 
with 
‘ k 
H, = '/2 Gn (. -—~h~ & ‘) — (26,/ne)J, 
2 n 9 k 
H,’ = _ Bs rm J E -)=— a4 | 
n o 2 


Thus the A, are known for this specific case. 
Then the a, are, and from Eq. (17’) is obtained the 
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circulation, so that the problem is now essentially 
solved. 
CALCULATION OF THE LIFT 


AT can be followed without change until Eq. (25). 
In Eq. (26) the final aw’ needs be replaced by (aw)’. 
As further details consist of a routine application of the 
results in the earlier sections, they are omitted. Letting — 
In be the steady lift for the mean conditions, i.e., 
here 


Lo = raopw?(2b) = hpwy’? 
and L the lift sought, the final result is 
L 


o? k 
) in 1+ 5 +0(- 5.4) + 


[Ke + Bi) + h:]cos +[((¢ - ah (1 + =) + 


oF + 1 | sint + (3 oB2 + h) cos 2r + 


(5 cB + h') sin 27 + >> (1, cos nr + 1,’ sin nr) 

“2 n=3 

where 

ln + ily’ = —2mi-"D { Fala n(no) — Jn-m(ne)] + 
n=1 


4G n(Jn4m(no) + Ja nm(na) }} 
F, + iG, = = [Flnk) + iG(nk)\H y+ Hy) 


H, and H,,’ as above 
tT = oat. 

Computing the elementary lift for comparison 
L 2 2 : 
r= 1+ 5+ 06 + fi(1 +2) cos'r + 

[20 — B, (1 + 2/407)] simr +.... 

it should be noted that here, contrary to AT, the mean 
lift is not necessarily the same in both cases. 


CALCULATION OF THE MOMENT 


A climbing moment will be taken as positive. Fol- 
lowing the procedure used hitherto, the moment 1/ is 
decomposed into M,; + M,, with 


My; = = ow f yy x dx 
— (p/2)(d/dt) f°,(b? — x*)y dx 


The integral in M, has already occurred in the 
lift calculations of AT. From AT, Eg. (25), it 
is 


M, = 


br — 1 = —xb'aw + ib Sa,S,e*"”* 


Thus J;, the integral appearing in M,, presents the 
only new difficulty. Evaluating it: 
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I, = BSttsin® 0(0, t) do 

1/,53 f,*(1 — cos 26)7(6, ¢) sin 6 dé 
1/.b°T — (9/4)b%ce 

1/b°T + (/4)b4a — (2/4)b*(b3 — Dy) 


Using Eq. (16”) and performing reductions by the 
usual techniques, there is obtained 


In = (x /4)b4& + (b?/2)VragT pe”? 


ll 


where 
io = | 
T, = ¢(nk), (n > 0) 
y ae = Fs 

and ? 


fx) =1+i fr*e”xX 


(N= 23)+1Jo 


Now, marshalling of the results leads to 





M =p rb?aw? — * b‘pa — phwlA, X 
8 n +0 


(2 + jents) inkW/b 

4—____—__— Je 
Ra 

In an appendix it will be shown that the expression 
in parentheses is equal to '/,[1 — C(nk)]. 





1946 


Putting 
Mo = (b/2)Lo = prbrapwy? “i (b/2)hpwe? 
there is obtained 
M a w 2 1 a 1 d i" 
aT ae ae, ba Rk? ed t@ v 
Mo aX (=) ? 8 (<) + i dr 2, n? Xx 
(i. - 1H,,') [C(nk) ai 1 Jeit(* ~ e080) 


All that remains is the arrangement as a Fourier series, 
This is a simple matter as it has already been done for 
the infinite series if the bracket is thought of as being 
replaced by C(mk). Then the final result is 


M o2 o? Be 

_~ (1 +2 + oh) +| o.(1 +748) +1| 
— + | 20 +0,(1 + co? +*) +! |sinr + 
| -¢ fs + ps) a | cos 2r + (cB: + ty’) sin 27 + 


(- = +1) cos 3r + (- se + n!) sin 37 + 


>> (tm cos mr + tm’ sin mr) 
m=4 
where tm + itm’ is given by the same expression as 
ln + tm’ except that F(nk) is to be replaced by 
F(nk) — 1. 


Appendix 


In the region R, (J[x] < 0, x + 0) the following 
integrals exist 


E(x) = x fire ™*(1 +X — Vd? + 2A) dA = xar(x) 


g(x) = 1+ ix fore [2(1 + AWA? + QA — QA? — 
4. — 1] dd 


because both brackets = O0(A), (A ©). 
region I[x] <c <0 


Then for any 


t(x) = 1 — Qixa(x) + ix fp %e7* [2A A? + DD _— 
oo 1)+ 1] dx ead 2 —_ 2ixa(x) + 2x a’ (x) 


and so 


bol 
Lemuel 


E(x) + ; xt (x) of (1 2 § xa) + Zax) +; 





From the principles of analytic continuation this latter 
equality must be valid in R, in particular on the posi- 
tive real axis. Here, from AT, Appendix II, Eq. (6), 
1 


a(x) = BL —_- — 
x 


OW ix (2) /.. 
3 a5! hh” (x) 
x x 


Thus, a simple calculation gives 


xl 
a 5° 
~ — 


E(x) + zat (e) =— @ FT, (x), (x > 0) 


Then, using Eqs. (5) and (7) from the same appendix, 
there is readily obtained 


fed ee & iH)® (x) |- 
v(x) 21H) +i) 1 


5 (C(x) =U 








= ben hep he hesk Ek! 





aor srr a 


Ysetonm roasos 


sinh 


Pom 


SUR weno 


moon 





Wo? 


ef (* — conn) 


‘ier series, 
1 done for 
fas being 


+h 
sin; + 
n 27 + 
in 37 + 
»’ Sin mr) 


ession as 
laced by 


his latter 
the posi- 
Eq. (6), 


> 0) 


ippendix, 


(x) - 1] 





INSTITUTE ofthe AERONAUTICAL SCIENCES 


President 
ARTHUR E. RAYMOND 


Vice-Presidents 
P. R. BASSETT 
J. H. KINDELBERGER 
Joun C. LESLIE 
BURDETTE S. WRIGHT 


Director 
S. PauL JOHNSTON 


Executive Vice-President 


BENNETT H. HORCHLER 


E. E. ALDRIN 
PETER ALTMAN 

Rex B. BEISEL 
LAWRENCE D. BELL 
GEoRGE W. BRADY 
W. A. M. BuRDEN 


W. M. ALLEN 
WALTER BEECH 
ERNEST R. BREECH 
RaLpH S. DAMON 


. DonaLp W. DovucLas 


R. H. FLEEt 
Jack FRYE 
RoBertT E. Gross 


Joun C. ADAMS 
GRIFFITH BREWER 
W. A. M. BurpEN 
E. U. Conpon 

L. C. CRAIGIE 


H. H. ARNOLD 
DonaLp W. DoucLas 
HucH L. DryDEN 

W. F. DuRAND 

Srr RICHARD FAIREY 


E. E. ALDRIN 

Hart O. BERG (Deceased) 
W. A. M. BurpEN 

E. L. Corp 

DonaLp W. DouGLas 
SHERMAN M. FarRcHILp 


Joun D. AKERMAN 

Epwin E. ALDRIN 

Harry G. ARMSTRONG 

Kart ARNSTEIN 

J. L. ATwoop 

Pau. S. BAKER 

Ratpeu S. BARNABY 

P. R. BASSET? 

WELLWoop E. BEALL 

Rex B. Bgersev 

G. M. BELLANCA 

Don R. BERLIN 

Henry A. BERLINER 

Wittram BoLiay 

B. C. BouLToNn 

Grorce W. BrRapy 

Lyman J. Briccs 

W. G. BRoMBACHER 

CuarLes P. BuRGEsS 

FrRanK W. CALDWELL 

KENNETH CAMPBELL 

Cuartes H. CHATFIELD 

ALLAN CHILTON 

V. E. CLark 

FRANKLIN R. CoLLBoHM 

Cuarvtes H. Cotvin 

Kart T. Compton 

Joun W. Crow ey, JR. 

Ratpg S. DAMON 

Luis DE FLoREz 

Smira J. De FRANCE 

. P. Den HarTOG y 
ALTER S. DIEHL 

James H. Doo.itTLe 

C.S. Draper 


Past-Presidents 


J. C. HUNSAKER 
CHARLES L. LAWRANCE 
DoNnaLp W. DouGLas 
GLENN L. MARTIN 
CLaRK B. MILLIKAN 
T. P. WricHT 
GEORGE W. LEwIS 
JaMEs H. DOooLittr_e 
FRANK W. CALDWELL 
Hav L. HrBBarRD 
Hucu L. DrypDEN 

R. H. FLEET 
CHARLES H. CoLvIN 


COUNCIL 


LESTER D. GARDNER, Chairman 


Assistant to the President 
GEORGE R. FORMAN 


Secretary 
ROBERT R. DEXTER 


Treasurer 
C. S. JonEs 


Controller 


JosePH J. MAITAN 


CHARLES H. CoLviIn 
JaMEs H. DooLitTrL_e 
WruiaM K. EBei 
SHERMAN M. FaIRCHILD 
DONALD W. FINLAY 


CHARLES FROESCH 
Leroy R. GRUMMAN 
R. Paut HARRINGTON 
W. WALLACE KELLETT 
Orto E. KiRCHNER 


ADVISORY COMMITTEE 


R. M. Hazen 

H. MANSFIELD HORNER 
J. C. HUNSAKER 

Croit HuNTER 

TH. VON KARMAN 

I. M. Lappon 

P. W. LITCHFIELD 


HONORARY 


Ira C. EAKER 

Harry F. GuGGENHEIM 
Henry G. HovuGHTon 
Georce B. PEGRAM 


F. W. MacIn 
ALFRED MARCHEV 
GLENN L. MarTIN 
Tuomas A. MoRGAN 
R. C. Murr 

Joun K. NORTHROP 
W. A. PATTERSON 


MEMBERS 


L. WeLcH PocuE 

J. LAURENCE PRITCHARD 
ARTHUR W. RADFORD 

F. W. REICHELDERFER 


HONORARY FELLOWS 


Str Roy FEpDEN 

J. C. HUNSAKER 
Srr MELVILL JONES 
Georce W. Lewis 
GLENN L. MARTIN 


GEorRGE J. Meap 

Srr FREDERICK HANDLEY PaGE 
D. R. Pye 

Icor I. Srkorsky 

G. I. TayLor 


BENEFACTORS 


Srr RICHARD FAIREY 
R. H. FLEet 


FLORENCE GUGGENHEIM (Deceased) 


Cc. D. Hanscom 
Joun P. V. HEINMULLER 
ALEXANDER KLEMIN 


AMERICAN 


Ivan H. Driccs 
WitiiaM K. Esper 

C. L. EGtvept 
FREDERIC FLADER 
CHARLES FROESCH 

JacK FRYE 

GARLAND FULTON 
Lester D. GARDNER 
RicHarp C. GAZzLey 
MICHAEL E. GLUHAREFF 
MELVIN N. GouGH 
Leroy R. GRUMMAN 
Harovp R. HARRIS 

R. M. Hazen 

ALBERT F. HEGENBERGER 
Epwarp H. HEINEMANN 
S. D. Heron 

Haut L. HrBBarp 

L. S. Hospss 

Pau. E. HovGAARpD 
WILLIAM HovGAARD 
ROBERT INSLEY 
EASTMAN N. JACOBS 
CLARENCE L. JOHNSON 
ALEXANDER KARTVBELI 
Paut H. KeEMMER 

J. H. KInDELBERGER 
Orto E. KrrcHNER 
ARTHUR L. KLEIN 
ALEXANDER KLEMIN 

W. B. KLemMpPerer 
Pavut KOLLSMAN 

Otto C. Kopren 

B. V. Korvin-KRovuKOvsky 


PauL Ko_isMAN 

BELLA C. LANDAUER 

WILLIAM LANDAUER 

GROVER LOENING 

Epmunp C. Lyncu (Deceased) 
F. W. Mactn 


FELLOWS 


S. M. Kraus 

I. M. Lappon 

E. S. LAND 

CHARLES L. LAWRANCE 
Joun G. LEE 

Joun C. LESLIE 
CHarveEs A. LINDBERGH 
WituraM LITTLEWOOD 
GROVER LOENING 

W. RANDOLPH LoveE.ace, II 
Eucene E. LuNpQuist 
R. D. MacCart 

Joun R. MARKHAM 
ERLE MARTIN 

Wa. H. McAvoy 

Cc. J. McCartuy 
CLARK B. MILLIKAN 
ROBERT J. MINSHALL 
Sanrorp A. Moss 
Joserxu S. NEWELL 
ALFRED S. NILES 

Joun K. NortTHrRop 
ARTHUR NuTT 
SHATSWELL OBER 
RoBERT R. OSBORN 

W. Bat_tey OswaLp 
Gerorce A. PAGE, JR. 
F. W. PawLowskI 

F. W. PENNOYER, JR. 
A. A. PRIESTER 
ARTHUR E. RAYMOND 
F. W. RBICHELDERFER 
E.trott G. Rem 

H, J. E. Rem 


JEROME LEDERER 
ROBERT J. MINSHALL 
Ear. D. OsBoRN 
ELMER A. SPERRY, JR. 
ERNEST G. Stout 
Frep E. WEICK 


E. T. Price 

F. W. REICHELDERFER 
LAURANCE S. ROCKEFELLER 
Frep H. Rour 

T. CLauDE RYAN 

Guy W. VAUGHAN 

EDWARD P. WARNER 


H. B. SALLapa 

CARL SPAATZ 

L. C. STEVENS 
NATHAN F. TWINING 


TH. von KARMAN 
Epwarp P. WARNER 
ORVILLE WRIGHT 

T. P. Wricut 


GLENN L. MARTIN 

Grorce J. Mgap 

S. A. Reep (Deceased) 
LAURANCE S. ROCKEFELLER 
Epwarp P. WARNER 


H. C. RiIcHARDSON 
LEONARD E. Root 

C. E. RosENDAHL 
Witiram T. SCHWENDLER 
James M. SHOEMAKER 
Mac SHort 

Icor ALExIs SIKORSKY 
SELDEN B. SPANGLER 
Ev_MER A. Sperry, JR. 
Joun Strack 

WitiraM B. Stout 
Harry A. SUTTON 

C. Fayette TAYLOR 
Epwarp S. TAYLor 
Puiip B. TAYLOR 

Tu. THEODORSEN 

T. E. TIrL_IncHAsT 
S. TIMOSHENKO 
STaRR TRUSCOTT 

L. B. TUCKERMAN 
Rapa H. Upson 

J. G. VINCENT 
RICHARD VON MISES 
FRANK L. WATTENDORF 
P. V. H. Weems 
Frep E. Weick 

R. D. WeEYERBACHER 
Joun B. WHEATLEY 
ANDREW VY. WILLGOOS 
DonacLp H. Woop 
RoBeERT J. Woops 
Joun E. YOUNGER 
ALBERT F. ZAHM 
STEPHEN J. ZAND 














An Appreciation 


TO CORPORATE MEMBERS OF THE 
INSTITUTE OF THE AERONAUTICAL SCIENCES 


Listed below are the leading aeronautical and industrial companies that make possible the publication of the 


Journal of the Aeronautical Sciences. 


They enable the Institute to publish scientific and technical papers in the Journal and to devote its pages ex- 


clusively to editorial material. 


Dues from individual members which would normally be used to publish this scientific and engineering 


monthly are thus released for additional services to the profession and industry. 
The Aeronautical Engineering Review and the Aeronautical Engineering Catalog, together with the Journal, provide 
a comprehensive technical information service for the benefit of the aircraft and related industries. 


ACADEMY OF AERONAUTICS 
AERO INSURANCE UNDERWRITERS 
AEROJET ENGINEERING CORPORATION 
AEROQUIP CORPORATION 
AGAWAM AIRCRAFT PRODUCTS, INC. 
AIRCHOX COMPANY DIVISION, JOYCE AVIATION, INC. 
AIRCRAFT RADIO CORPORATION 
AIREON MANUFACTURING CORPORATION 
AIRESEARCH MANUFACTURING COMPANY 
AIRESEARCH MANUFACTURING COMPANY OF 
ARIZONA, INC. 
AIREX MANUFACTURING COMPANY, INC. 
ALLIS-CHALMERS MANUFACTURING COMPANY 
ALLISON DIVISION, GENERAL MOTORS CORPORATION 
ALUMINUM COMPANY OF AMERICA 
AMERICAN AIRLINES SYSTEM 
AMERICAN OVERSEAS AIRLINES, INC. 
AMERICAN BOSCH CORPORATION 
AMERICAN PHENOLIC CORPORATION 
ASSOCIATED AVIATION UNDERWRITERS 
ASSOCIATED FOUNDRIES & MANUFACTURERS, INC. 
ATLAS SUPPLY COMPANY 
BAKER STEEL & TUBE COMPANY 
BEECH AIRCRAFT CORPORATION 
BELL AIRCRAFT CORPORATION 
BENDIX AVIATION CORPORATION 
BENDIX PRODUCTS DIVISION 
BENDIX RADIO CORPORATION 
ECLIPSE-PIONEER DIVISION 
EXPORT DIVISION 
FRIEZ INSTRUMENT DIVISION 
PACIFIC DIVISION 
SCINTILLA MAGNETO COMPANY 
THE BG CORPORATION 
BORING AIRCRAFT COMPANY 
BOEING SCHOOL OF AERONAUTICS 
BREEZE CORPORATIONS, INC. 
ESSEX TOOL & DIE COMPANY 
FEDERAL LABORATORIES, INC. 
CAL-AERO TECHNICAL INSTITUTE, A DIVISION OF 
GRAND CENTRAL AIRPORT COMPANY 
CARRIER CORPORATION 
CESSNA AIRCRAFT COMPANY 
CHANDLER-EVANS CORPORATION 
CHASE NATIONAL BANK OF THE CITY OF NEW YORK 
THE CLEVELAND GRAPHITE BRONZE COMPANY 
THE CLEVELAND PNEUMATIC TOOL COMPANY 
AUTOMOTIVE-AIRCRAFT DIVISION 
CLIFFORD MANUFACTURING COMPANY 
CONSOLIDATED VULTEE AIRCRAFT CORPORATION 
NASHVILLE DIVISION 
STINSON AIRCRAFT DIVISION 
CONTINENTAL MOTORS CORPORATION 
COX AND STEVENS AIRCRAFT CORPORATION 
CULVER AIRCRAFT CORPORATION 
CURTISS-WRIGHT CORPORATION 
AIRPLANE DIVISION 
DEVELOPMENT DIVISION 
PROPELLER DIVISION 
WRIGHT AERONAUTICAL CORPORATION 
DOAK AIRCRAFT COMPANY, INC. 
DOUGLAS AIRCRAFT COMPANY, INC. 
EL SEGUNDO PLANT 
LONG BEACH PLANT 
THE DOW CHEMICAL COMPANY 
DZUS FASTENER COMPANY, INC. 
EASTERN AIR LINES, INC. 
EATON MANUFACTURING COMPANY 
THOMAS A. EDISON, INCORPORATED, INSTRUMENT 
DIVISION 
EDO AIRCRAFT CORPORATION 
THE ELECTRIC AUTO-LITE COMPANY 





ELECTROL INCORPORATED 
ENGINEERING AND RESEARCH CORPORATION 
ETHYL CORPORATION 
THE FAFNIR BEARING COMPANY 
FAIRCHILD CAMERA & INSTRUMENT CORPORATION 
FAIRCHILD AERIAL SURVEYS, INC. 
FAIRCHILD ENGINE AND AIRPLANE CORPORATION 
DURAMOLD DIVISION 
FAIRCHILD AIRCRAFT DIVISION 
RANGER AIRCRAFT ENGINES DIVISION 
FEDERAL TELEPHONE AND RADIO CORPORATION 
FIRESTONE AIRCRAFT COMPANY 
FLEETWINGS DIVISION, KAISER CARGO, INC 
FLETCHER AVIATION CORPORATION 
G & A AIRCRAFT, INC. 
GENERAL AIRCRAFT EQUIPMENT, INC. 
GENERAL AVIATION EQUIPMENT COMPANY, INC 
GENERAL ELECTRIC COMPANY 
GENERAL INSTRUMENT CORPORATION 
GENERAL MOTORS CORPORATION 
AC SPARK PLUG DIVISION 
AEROPRODUCTS DIVISION 
BUICK MOTOR DIVISION 
CADILLAC MOTOR CAR DIVISION 
CHEVROLET MOTOR DIVISION 
DELCO PRODUCTS DIVISION 
DELCO-REMY DIVISION 
EASTERN AIRCRAFT DIVISION 
FISHER BODY DIVISION 
FRIGIDAIRE DIVISION 
HARRISON RADIATOR DIVISION 
RESEARCH LABORATORIES DIVISION 
ROCHESTER PRODUCTS DIVISION 
THE B. F. GOODRICH COMPANY 
THE GOODYEAR TIRE & RUBBER COMPANY 
THE GRAY MANUFACTURING COMPANY 
GRUMMAN AIRCRAFT ENGINEERING CORPORATION 
GUIBERSON DIESEL ENGINE COMPANY 
W. & L. E. GURLEY 
HANNA ENGINEERING WORKS 
HASKELITE MANUFACTURING CORPORATION 
HAWAIIAN AIRLINES LIMITED 
THE HILLIARD CORPORATION 
THE [INTERNATIONAL NICKEL COMPANY 
JACK & HEINTZ, INC. 
JACOBS AIRCRAFT ENGINE COMPANY 
JOHNS-MANVILLE SALES CORPORATION 
CASEY JONES SCHOOL OF AERONAUTICS, INC. 
KELLETT AIRCRAFT CORPORATION 
KENYON INSTRUMENT COMPANY, INC. 
WALTER KIDDE & COMPANY, INC. 


KOLLSMAN INSTRUMENT DIVISION, SQUARE D COM- 
ANY 


LANGLEY CORPORATION 

LAVELLE AIRCRAFT CORPORATION 

LEAR INCORPORATED 

LIBERTY AIRCRAFT PRODUCTS CORPORATION 
LINK AVIATION DEVICES, INC. 

THE LIQUIDOMETER CORPORATION 

LOCKHEED AIRCRAFT CORPORATION 
LONGINES-WITTNAUER WATCH COMPANY, INC 
MAGNAFLUX CORPORATION 

THE MARQUETTE METAL PRODUCTS COMPANY 
THE GLENN L. MARTIN COMPANY 

THE W. L. MAXSON CORPORATION 

WARREN McARTHUR CORPORATION 
McDONNELL AIRCRAFT CORPORATION 
MENASCO MANUFACTURING COMPANY 
MICROMATIC HONE CORPORATION 
MINNEAPOLIS-HONEYWELL REGULATOR COMPANY 
MOORE DROP FORGING COMPANY 

NATIONAL CITY BANK OF NEW YORK 
NATIONAL CREDIT OFFICE, INC. 


THE NEW YORK AIR BRAKE COMPANY 
NORMA-HOFFMANN BEARINGS CORPORATION 
NORTH AMERICAN AVIATION, INC. 
NORTH AMERICAN AVIATION, INC., OF TEXAS 
NORTHROP AIRCRAFT, INC. 
NORTHWEST AIRLINES, INC. 
OTTO AVIATION CORPORATION 
OWENS-CORNING FIBERGLAS CORPORATION 
PAN AMERICAN WORLD AIRWAYS SYSTEM 
THE PARKER APPLIANCE COMPANY 
PENNSYLVANIA-CENTRAL AIRLINES CORPORATION 
PESCO PRODUCTS COMPANY DIVISION, BORG-WARNER 
CORPORATION 
PHILLIPS PETROLEUM COMPANY 
PIONEER PARACHUTE COMPANY, INC. 
THE PURE OIL COMPANY 
REPUBLIC AVIATION CORPORATION 
J. P. RIDDLE COMPANY 
A. V. ROE CANADA LIMITED 
JOHN A. ROEBLING’S SONS COMPANY 
ROHR AIRCRAFT CORPORATION 
ROOSEVELT FIELD, INC. 
THE RYAN AERONAUTICAL COMPANY 
SCIAKY BROTHERS 
SCOTT AVIATION CORPORATION 
SHELL OIL COMPANY, INC. 
SIMMONDS AEROCESSORIES, INC. 
SKYDYNE, INC. 
SOCONY-VACUUM OIL COMPANY 
SOLAR AIRCRAFT COMPANY 
SPERRY GYROSCOPE COMPANY, INC. 
SQUARE D COMPANY 
STANDARD OIL COMPANY OF CALIFORNIA 
STANDARD OIL COMPANY (INDIANA) 
STANDARD OIL COMPANY OF NEW JERSEY 
SWEDLOW AEROPLASTICS CORPORATION 
TELEFLEX LIMITED 
THE TEXAS COMPANY 
THOMPSON PRODUCTS, INC 
TINNERMAN PRODUCTS, INC. 
TITEFLEX, INC. 
TRANSCONTINENTAL & WESTERN AIR, INC. 
TRIPLETT & BARTON, INC. 
UNION CARBIDE AND CARBON CORPORATION 
BAKELITE CORPORATION 
HAYNES STELLITE COMPANY 
LINDE AIR PRODUCTS COMPANY 
NATIONAL CARBON COMPANY 
UNITED AIRCRAFT CORPORATION 
CHANCE VOUGHT AIRCRAFT DIVISION 
HAMILTON STANDARD PROPELLERS DIVISION 
PRATT & WHITNEY AIRCRAFT DIVISION 
PRATT & WHITNEY AIRCRAFT CORPORATION OF 
MISSOURI 
SIKORSKY AIRCRAFT DIVISION 
UNITED AIR LINES, INC. 
UNITED STATES AVIATION UNDERWRITERS, INC 
UNITED STATES RUBBER COMPANY 
DOMINION RUBBER COMPANY, LTD. 
THE VARIETY AIRCRAFT CORPORATION 
VICKERS, INC. 
VIDAL CORPORATION 
THE WACO AIRCRAFT COMPANY 
WARNER AIRCRAFT CORPORATION 
THE WEATHERHEAD COMPANY 
WESTERN AJR LINES, INC. 
WESTINGHOUSE ELECTRIC CORPORATION 
WESTON ELECTRICAL INSTRUMENT CORPORATION 
WYMAN-GORDON COMPANY 
YOUNG RADIATOR COMPANY 























RATION 
-WARNER 


ISION 


TION OF 





